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FOREWORD 


This  report  describes  work  performed  under  Contracts  AP  33(6l6)- 
7757  and  AP  33(657)-76l7,  "A  Research  Program  for  Understanding 
the  Mechanisms  of  Flame  Inhibition,”  and  AP  33(6l6)-7458,  "Fire- 
Resistant  High  Temperature  Hydraulic  Fluids  .  '  The  former  contract 
was  initiated  under  Project  No.  6075^  "Plight  Vehicle  Hazard 
Protection,"  Task  No.  607505,  "A  Research  Program  for  Understanding 
the  Mechanisms  of  Flame  Inhibition."  The  latter  was  initiated 
under  Project  No.  73^0,  "Non-Metalllc  and  Composite  Materials," 

Task  No.  734008,  "Power  Transmission,  Heat  Transfer  Fluids,  and 
Other  Forms  of  Energy  Transfer  Fluids . " 

The  contracts  were  performed  at  the  Dayton  Laboratory  of  Monsanto 
Research  Corporation.  The  first  two  contracts  were  sponsored  by 
the  Plight  Accessories  Laboratory,  Aeronautical  Systems  Division, 
Wright- Patterson  Air  Force  Base,  Ohio  with  Mr.  Benito  P.  Botteri 
serving  as  project  engineer.  The  third  contract  was  sponsored  by 
the  Directorate  of  Materials  Processes,  Aeronautical  Systems 
Division,  Wrlght-Patterson  Air  Force  Base,  Ohio  with  Mr.  Harold 
Adams  as  project  engineer. 

For  Monsanto,  the  computer  program  was  written,  by  Dr.  G.  H.  Rlngrose 
of  Monsanto  Research  Corporation  and  Dr.  D.  R.’  Miller  and 
Dr.  A.  C.  Pauls  of  Monsanto  Chemical  Company.  Dr.  G.  B.  Skinner 
served  as  project  leader. 

The  authors  are  indebted  to  Dr.  W.  C.  Hammann  of  Monsanto  Chemical 
Company  for  the  hand  calculations  that  led  to  the  adoption  of  the 
general  mathematical  approach.  This  contribution  and  the  many 
helpful  suggestions  made  by  numerous  people  are  acknowledged  with 
gratitude.  Particular  assistance  was  given  by  Mr.  D.  J.  Kaufman, 
also  of  Monsanto  Chemical  Company,  who  was  instrumental  in  trouble 
shooting  and  expediting  the  machine  solutions. 

The  routines  were  originally  programmed  and  run  on  the  IBM  704 
computer  in  the  Research  and  Engineering  Division  of  Monsanto 
Chemical  Company,  St.  Louis.  Later  the  programs  were  adapted  for 
solution  by  an  IBM  7090,  at  the  Aeronautical  Systems  Division, 
Wrlght-Patterson  AFB,  Ohio. 


ABSTRACT 


Two  digital  computer  routines  were  developed  to  process  flame 
speed  data  resulting  from  the  burning  of  compounds  in  air  and 
oxygen,  and  to  correlate  particular  structural  configuration 
with  flame  speed . 

In  both  routines  a  high  degree  of  flexibility  has  been  in¬ 
corporated  to  assure  efficient  utilization  under  several 
forseeable  circumstances. 

The  first  routine,  PSC,  processes  the  raw  experimental  data 
to  obtain  flame  speeds,  equivalence  ratios,  and  the  equivalence 
ratio  at  the  maxlmxim  flame  speed .  This  Information  is  stored 
on  a  master  magnetic  tape  for  subsequent  calculations. 

The  second  routine,  PSR,  permits  selection  of  specific  data 
groups  from  the  master  tape  for  analysis.  A  linear  model 
was  chosen  for  the  correlation. 

This  technical  documentary  report  has  been  reviewed  and  is 
approved . 


flELIAH  TJr'SAVAGE' 

Chief,  Environmental  Branch 
Plight  Accessories  Laboratory 
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FLAME  SPEED  DATA  REDUCTION  AND  CORRELATION  USING  A  DIGITAL  COMPUTER 


I.  INTRODUCTION 


Two  Air  Force-supported  programs  at  this  laboratory  concern  the 
relative  rates  of  combustion  of  a  number  of  chemical  compounds. 
One  program  /AF  33(6i6)-745^  deals  with  the  synthesis  of 
potential  high  temperature  hydraulic  fluids  and  the  other 
/ffF  33 (616) -76127  with  increasing  our  understanding  of  the 
mechanisms  of  flame  inhibition,  the  ultimate  objective  being 
the  development  of  agents  for  extinguishing  propellant  fires. 

Searches  for  effective  fire-resistant  fluids  or  propellant  fire 
extlngulshants  would  be  systematized  and  expedited  if  a  general 
method  for  characterizing,  measuring,  and  predicting  combusti¬ 
bility  were  at  hand.  Unfortunately,  most  flammability  tests 
and  specifications  in  use  are  so  application-oriented  or  other¬ 
wise  restricted  as  to  be  of  only  limited  value  for  general  use. 

A  measure  of  the  Intrinsic  ability  of  a  material  to  support  or 
inhibit  combustion,  as  free  as  possible  from  influences  peculiar 
to  the  testing  procedure,  is  the  first  requisite.  Given  such 
a  measure,  there  is  then  the  possibility  of  relating  combusti¬ 
bility  to  the  molecular  structures  of  materials  tested  and  of 
predicting  therefrom  the  combustibility  of  untested  candidates. 

This  report  describes  computation  routines  designed  to  support 
the  above  approach.  Flame  speed  -  a  state  property  of  a 
combustible  gas  mixture  -  is  the  intrinsic  measure  of  com¬ 
bustibility  selected  for  use  in  characterization,  correlation, 
and  prediction.  Earlier  work,  although  limited  to  pure  hydro¬ 
carbons  burned  in  air,  tends  to  confirm  the  approach  adopted 
here . 

The  routines  are  programmed  in  the  FORTRAN  II  language  for 
IBM  7090  solution. 


Manuscript  released  by  authors  February  1963  for  publication  as 
an  ASD  Technical  Documentary  Report. 
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II.  THEORETICAL  ASPECTS 


Theoretical  flame  speed  is  commonly  defined  as  the  subsonic 
rate  of  perpendicular,  propagation  of  an  Infinite  plane  flame 
front  through  the  quiescent  combustible  gas.  External  forces 
(for  example,  gravitational  or  magnetic  fields)  are  presumed 
absent  or  Insignificant,  and  the  flame-induced  flow  is  assumed 
to  be  laminar  and  one-dimensional . 

When  defined  in  this  Idealized  manner,  flame  speed  is  a  function 
only  of  the  initial  temperature,  pressure,  and  composition  of 
the  combustible  gas.  It  is  thus  an  Intensive  state  property 
of  the  gas.  Just  as  are  density,  enthalpy,  viscosity,  refractive 
index,  and  the  like .  Moreover,  flame  speed  is  a  transport 
property  dependent  upon  the  rates  of  energy,  momentum,  mass, 
and  chemical  exchange  accompanying  combustion. 

Prediction  of  flame  speed  from  gas  temperature,  pressure,  and 
composition  is  a  long-standing  goal  of  combustion  research. 

A  reasonably  general  and  accurate  method  would  be  extremely 
useful  in  solving  a  range  of  practical  and  theoretical  problems.. 
Combustor  and  propulsion  system  design  may  be  cited,  in  addi¬ 
tion  to  the  materials  search  problems  underwriting  the  work  of 
this  report. 

In  principle,  flame  speeds  can  be  predicted  theoretically. 
Simultaneous  solution  of  the  differential  equations  of 
change  is  possible  when  all  Important  physical-chemical 
properties  (e.g.,  heats  of  reaction,  reaction  rate  constants, 
density,  dlffuslvltles,  emissivlty,  etc.)  are  known.  Un¬ 
fortunately,  many  of  these  properties  are  not  yet  known. 
Ccmbustlon  reaction  mechanisms  are  normally  quite  complex  and 
are  consequently  poorly  understood.  Much  remains  to  be  done, 
even  in  the  Identification  of  intermediate  chemical  species. 

While  theory  is  of  little  practical  utility,  empirical  cor¬ 
relations  offer  considerable  promise.  Hibbard  and  Plnkel 
(Ref.  l)  achieved  a  good  correlation  of  maximum  flame  speed 
of  37  hydrocarbons  (mixed  with  air)  versus  the  concentrations 
of  the  various  C-H  bond  types  present.  The  relation  used  was 
of  the  form 


%ax  ‘^J^  J  ^  ®  J  ^  max 


(1) 


The  concentration  (cj)niax  of  f^e  J  type  bond  is  determined  at 
the  fuel-to-oxldlzer  ratio  which  yields  the  maximum  flame 
velocity  Umax  fhe  corresponding  influence  coefficient  for 
this  bond  is  designated  bj.  The  average  per  cent  deviation 


between  the  measured  and  predicted  maximum  flame  speeds  was 
about  2%  for  the  37  hydrocarbons  correlated. 

The  approach  of  Hibbard  and  Plnkel  was  extended  by  Hammann 
and  Blake  (Ref.  2)  In  calculation  of  flame  speed  coefficients 
of  the  various  bond  types  In  fuels  containing  oxygen,  nitrogen, 
sulphur,  boron,  and  silicon  In  addition  to  carbon  and  hydrogen. 
Their  correlations  were  based  upon  data  obtained  experimentally 
by  burning  142  model  compounds  and  upon  data  reported  by  Gibbs 
and  Calcote  (Ref.  3).  Several  modifications  of  the  correlation 
model.  Equation  1,  were  tested  and  Equation  lA  was  found  to 
provide  the  most  acceptable  correlation. 


The  number  of  J-type  contributors  occurring  In  each  molecule  of 
fuel  Is  designated  nj . 

While  the  agreement  between  the  predicted  and  measured  flame 
speeds  was  not  as  good  as  for  hydrocarbons  alone,  the  results 
convincingly  supported  the  general  utility  of  the  linear 
correlation  technique.  The  computer  routines  described  In 
this  report  were  used  by  Hammann  and  Blake. 


'^ax 
^°f ^max 


DEFINITIONS 


The  following  terms  are  used  frequently  throughout  the  report: 

Contributor  A  countable  structural  feature  of  a  fuel 
molecule  (e.g.,  each  hydrogen  bonded  to  a  primary  carbon  In 
H  H 

ethane,  HC-CH,  could  be  classed  as  a  contributor). 

H  H 

Contributor  number  1  An  Identifying  code  number  which 
Is  assigned  to  each  defined  contributor  (e.g.,  the  prlmary-H 
contributor  was  assigned  the  code  nianber  37). 

Contributor  count  nj  (or  n^j)  The  number  of  J-type 
contributors  contained  in  a  molecule  of  fuel  (or  In  a  molecule 
of  the  ith  fuel) (e.g.,  the  contributor  count  of  the  contributor 
prlmary-H  In  the  fuel  ethane  would  be  6). 

Fuel  A  material  that  forms  a  combustible  mixture  when 
mixed  with  an  oxidant.  Included  are  Impurities,  additives, 
and  all  compounds  not  Included  In  the  gaseous  oxidant.  For 
totalling  the  contributor  count  for  mixed  or  composite  fuels. 
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it  is  assumed  that  the  over-all  fuel  has  additive  properties  of 
the  compounds  present.  For  example,  if  is  the  mole  fraction 
of  compound  i  in  the  fuel,  and  n^^  is  the  count  for  contributor 
J,  then 

^  ^l^lj  ^2^2 j  +  •  •  • 

Oxidant  The  portion  of  the  combustible  mixture  which 
Includes  the  oxygen  and  inert  gases  but  excludes  the  fuel. 

Equivalence  Ratio  The  actual  fuel-to-oxygen  ratio 
divided  by  the  fuel-to-oxygen  ratio  stolchionietrically  required 
for  complete  combustion  of  the  fuel  oxidizer.  This  is  a 
measure  of  the  richness  or  leanness  of  the  flame . 

Standard  Error  Synonymous  with  the  statistical  quantity 
"standard  deviation." 

Data  Group  A  unit  of  data  pertaining  to  a  single  fuel- 
oxidant  combination  burned  at  various  equivalence  ratios  but 
otherwise  under  Identical  conditions.  The  components  of  a 
data  group  in  the' order  in  which  they  are  stored  by  the  com¬ 
puter  on  master  tape  6  are: 

(a)  Serial  Number  Assigned  sequentially  by  the  computer 
in  order  of  tape  location. 

(b)  Fuel  Name  Alphabetic  and/or  numeric  representation 
of  the  fuel,  e.g.,  (n-Pentene-2) .  A  maximum  of  12 
characters  (including  blanks)  is  permitted. 

(c)  Fuel  Numbers  Code  numbers  which  are  used  to  classify 
the  types  of  fuel . 

(1)  Fuel  Glass  Number  General  fuel  type  (e.g., 
01-organlc  aliphatic;  02-Qrganic  aromatic, 

00- inorganic ). 

(2)  Fuel  Group  Number  Denotes  subclassification 
of  fuel  type  (e.g.,  06-saturated  cyclic  compound). 

(3)  Fuel  Member  Nvunber  An  arbitrary  code  number, 
used  to  indicate  sequence  of  experimental  analysis 
of  a  particular  group. 

Thus  for  the  fuel  cyclopropane  the  fuel  number  is 
010601. 

(d)  Data  Source  N^i^ber  Designates  the  source  of  the 
data  (e.g.,  Oi-ilterature-,  02-experlmental) . 
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(fe)  Expertoental  Conditions  N\jmber  A  code  to  dls- 
tingulsh  data  taken  under  different  experimental 
conditions,  such  as  temperature,  pressure,  per  cent 
oxygen  In  oxidant,  etc. 

(f)  Flame  speed  at  unity  equivalence  ratio,  Ugtoc 
(cm/sec; . 

(g)  Maxlm\am  flame  speed,  (cm/sec). 

(h)  Fuel  concentration  at  unity  equivalence  ratio, 

(Cf)stoc  (molecules/cc). 

(l)  Fuel  concentration  at  the  conditions  of  maximum  flame 
speed  (Of>)niax  (™ol®cules/cc) . 

(j)  Equivalence  ratio  at  the  maximum  flame  speed,  0inax 

(k)  Niomher  of  different  defined  (i.e.,  coded)  contributors 
in  the  fuel  molecule. 

(l)  List  of  the  contributor  code  numbers  with  their 
respective  co\ints  for  the  fuel  considered. 
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III.  DESCRIPTIVE  OUTLINE  OF  THE  CALCULATIONS 


An  analysis  of  the  computations  required  to  reduce  the  experi¬ 
mental  data  and  perform  the  regression  analyses  led  to  the 
separation  of  the  problem  into  two  sections.  The  first  section 
processes  the  raw  experimental  data  and  calculates  the  flame 
speeds,  concentrations,  and  equivalence  ratios.  These  quanti¬ 
ties  are  then  stored  on  a  master  reel  of  magnetic  tape.  The 
second  section,  namely,  the  flame  speed  regression  section, 
uses  this  reel  of  magnetic  tape  as  its  input  data.  Prom  these 
data,  the  regression  coefficients  of  the  flame  speeds  and  other 
dependent  variables  are  determined. 

Corresponding  to  each  of  these  sections  is  a  Fortran  routine 
that  is  described  in  detail  in  a  later  section.  A  general 
discussion  of  these  routines  follows. 


DATA  REDUCTION  (Routine  FSC) 

The  data  reduction  routine  calculates  the  flame  speed  and 
related  quantities  from  the  raw  data,  tests  these  quantities 
for  inconsistencies,  and  prepares  or  corrects  the  master  reel 
of  magnetic  tape  with  acceptable  data. 

Flame  Speed  The  calculations  are  based  upon  the  conical 
(Bunsen)  burner  method  for  flame  speed  determination.  Measured 
flows  of  the  several  components  of  the  mixture  are  passed  in 
laminar  flow  through  a  burner  tube  at  a  controlled  temperature 
and  pressure.  A  roughly  conical  flame  front  is  stabilized  at 
the  mouth  of  the  tube . 

The  large  density  gradient  in  the  gas  stream  at  the  flame 
front  permits  the  use  of  schlleren  photographic  techniques  to 
record  the  flame  profile.  It  is  assumed  that  the  flame  is 
symmetrical  about  its  vertical  axis .  Figure  1  is  schlleren 
photograph  of  a  typical  flame  cone . 

The  flame  speed  is  determined  by  dividing  the  total  volijmetrlc 
gas  flow  rate  by  the  area  of  the  flame  front.  Since  the  flame 
cone  photographed  is  seldom  a  "right  circular  cone, "  its 
surface  area  must  be  determined  from  measurements  at  several 
intermediate  points .  To  do  so,  the  diameter  is  measured  at 
several  heights  of  the  cone  as  illustrated  in  Figure  2.  The 
surface  area  between  any  two  of  these  diameter  measurements 
is  assumed  to  be  that  of  a  frustrum  of  a  right  circular  cone. 
The  sum  of  the  surface  areas  of  these  frustra  gives  a  fairly 
accurate  value  of  the  true  flame  front  area.  Actually,  the 
local  flame  speed  is  slightly  above  thp  average  at  the  cone  tip 
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and  slightly  below  the  average  at  the  base.  Therefore,  the 
calculated  flame  speed  must  be  considered  an  average  value 
only.  Evidence  has  shown  that  there  is  a  slight  dependence 
of  flame  speed  on  the  total  volumetric  flow  rate  of  the  gases, 
the  burner  part  geometry,  and  the  burner  part  temperature. 
Slight  influences  are  also  made  by  the  method  used  to  define 
the  cone  area  (l.e.,  schlleren  shadow,  or  radiant  photography). 
However,  by  standardizing  on  one  of  the  techniques,  this 
variable  can  be  eliminated. 

Valid  Data  Teats  Once  a  series  of  flame  speeds  has 
been  obtained,  it  is  necessary  to  determine  the  conditions 
for  maximum  flame  speed  (since  the  correlations  proposed  are 
based  on  conditions  at  the  maximum  flame  speed).  This  is 
accomplished  by  fitting  an  empirical  curve  of  flame  speed  ys . 
equivalence  ratio  (u  vs.  0)  to  the  data  and  solving  the  equa¬ 
tion  for  the  maximum.  The  flame  speed  at  stoichiometric 
conditions  (0  =  l)  is  also  determined. 

If  a  maximum  flame  speed  lies  within  tne  data  range  the  data 
are  acceptable  for  further  tests .  If  no  maximum  or  a  mlnimvim 
is  found  within  the  data  range,  the  data  are  not  acceptable. 

Storage  on  Magnetic  Tape  Data  groups  termed  accept¬ 
able  are  stored  on  a  master  reel  of  magnetic  tape  for  use  by 
the  regression  routine. 

Options  The  computational  routine  contains  several 
options  depending  upon  the  form  of  the  input  data.  These 
options  permit  sections  of  the  computational  procedure  to  be 
omitted.  For  example,  a  set  of  flame  speed  vs.  equivalence 
ratio  data  available  from  previous  calculations  could  be 
entered  directly  into  the  sequence,  and  the  original  raw  data 
calculations  would  be  by-passed.  These  options  are  detailed 
in  Section  VI. 


FLAME  SPEED  REGRESSION  (Routine  FSC) 

A  multiple  linear  regression  model  similar  to  that  used  by 
Hibbard  and  Plnkel  (Ref.  l)(but  including  a  constant  term) 
was  selected: 


y  -  >>0  *  (2) 

Independent  Variables  In  equation  (2)  the  independent 
variables  were  chosen  to  be  the  contributor  counts  rather  than 
the  contributor  concentration.  This  simplifies  the  prediction 
and  data  handling,  and  adds  flexibility  by  making  the  "inde¬ 
pendent”  variables  more  independent.  This  is  true  because  the 
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contributor  concentrations  depend  upon  the  mixture  temperature, 
pressure,  equivalence  ratio,  and  oxidant  composition,  as  well 
as  upon  fuel  composition.  A  simple  relation  exists  between 
the  contributor  counts  and  concentration,  as  shown  in  Equation  3! 

Cj  =  molecules/cc  (3) 

where  Cf>  is  the  concentration  of  fuel  in  the  combustible  mix¬ 
ture.  Multiplying  Equation  2  by  Cf  gives  the  equivalent  of  the 
Hibbard -Pinkel  model  (b^  =  O). 

Dependent  Variables  For  the  reasons  mentioned  above, 
the  dependent  variable  y  would  be  set  equal  to  u-ax^^^f  ^jnax 
when  the  coefficients  comparable  to  those  publisned  by  Mobard 
and  Pinkel  are  wanted.  However,  this  is  not  the  only  dependent 
variable  of  possible  Interest,  In  predictions  of  max imimi  flame 
speed,  for  example,  it  is  necessary  to  predict  fuel  concentra¬ 
tion  (or  fuel-oxidant  ratio,  or  equivalence  ratio)  at  maximum 
flame  speed.  This  is  most  readily  done  by  correlating  fuel 
concentration  (etc.)  at  maximum  flame  speed  against  contributor 
counts ,  Correlations  at  stoichiometric  fuel-oxidant  ratio  may 
be  equally  interesting.  To  emphasize  that  such  other  dependent 
variables  can  be  handled  with  equal  facility,  the  indefinite 
dependent  variable  y  is  used. 

Coefficients  With  the  emphasis  so  far  given  to  fuel 
composition,  it  is  appropriate  to  recall  that  calculated  bj's 
are  not  constants.  Values  will  depend  on  several  or  all  of 
the  following  (depending  on  the  definition  of  y): 

1.  initial  mixture  temperature 

2.  initial  mixture  pressure 

3.  some  specification  of  fuel-oxidant  ratio 

4.  oxidant  composition 

5.  contributor  definitions 

6.  dependent  variable  definition 

The  above  list  Includes  the  primary  state  properties.  Since 
experimental  flame  speeds  (distinguished  from  theoretical 
flame  speeds)  are  always  obtained  under  conditions  differing 
from  the  ideal,  values  of  bj  will  be  related  to: 

7.  experimental  conditions  (oxidant  composition, 
temperature,  pressure,  etc . 

8.  source  of  data  (laboratory,  technique,  Journal 
reference,  etc. 

All  listed  factors  should  be  considered  when  sets  of  coeffi¬ 
cients  are  compared. 
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Block  Regressions  Obtaining  a  full  set  of  coefficients 
is  basically  a  stepwise  trial-and-error  process.  The  trlal-and 
error  enters  in  the  definition  of  contributors  (i.e.,  independ¬ 
ent  variables).  Since  it  is  not  known  at  the  outset  which 
structural  features  will  give  the  best  correlations,  it  is 
necessary  to  accommodate  definition  of  new  contributors  as  the 
analysis  progresses.  In  addition,  it  is  desirable  and  most 
efficient  to  calculate  certain  coefficients  using  only  a 
particular  class  of  fuels.  C-H  bond  coefficients,  for  example, 
are  most  readily  and  accurately  determined  from  flame  speed 
data  on  hydrocarbons.  If  accurately  determined,  a  coefficient 
should  be  applicable  unchanged  in  later  regressions. 

The  above  considerations  indicated  that  a  "block  regression" 
procedure  (certain  coefficients  held  constant)  would  be  most 
advantageous .  This  procedure  may  be  expressed  mathematically 
by  Equation  4: 

y  -  -  *>0  +  ■^uVu 

The  bs  coefficients  represent  those  whose  values  are  known, 
and  are  termed  the  "prespecified  coefficients."  The  b©  and 
b^  coefficients  are  those  to  be  determined  by  the  regression 
analysis.  The  constant  term  bo  may  be  either  calculated  or 
prespecified  as  zero. 

Block  Selection  Criteria  All  of  the  data  groups  that 
are  acceptable  for  regression  analysis  are  stored  on  a  single 
master  reel  of  magnetic  tape.  Seven  different  accept-re Ject 
tests  are  available  for  selection  of  these  data.  The  select¬ 
ing  criteria  comprise  (l)  data  group  serial  number,  (2)  fuel 
class  number,  (3;  fuel  class-group  number,  (4)  fuel  member 
number,  (5)  data  source  number,  (6)  experimental  conditions 
number,  and  (7)  fuel  contributor  content. 

For  flexibility,  these  tests  may  be  Ignored,  or  applied  singly 
or  in  any  combination.  The  tests  are  described  in  detail  in 
the  following  section. 
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IV.  MATHEMATICAL  OUTLINE  OF  THE  CALCULATIONS 


FLAME  SPEED  DATA  REDUCTION  CALCULATIONS  (Routine  PSC) 


Functions  Performed  The  functions  performed  by  Routine 
FSC  may  be  classed  In  five  groupings.  The  portions  of  the 
computer  program  used  in  each  function  are  enclosed  In 
parentheses . 


Function  1  (Subroutines  EXPD  1  or  EXPD  3) 


Olven:  (a) 

jb) 

(c) 

jd) 

(e) 

(f) 


Combustion  mixture,  temperature  and 
pressure 

Fuel  stoichiometric  oxygen  demand 
Mole  fraction  of  O2  In  oxidant 
Molar  volume  of  fuel 
Volumetric  flow  rate  of  oxidant 
Flame  envelope  dimensions 


Calculate:  (a) 
jb) 

(c) 

(d) 


Volumetric  flow  rate  of  mixture 
Area  of  flame  envelope 
Flame  speed 
Equivalence  ratio 


(a)  Volumetric  flow  rate  of  the  mixture,  Qm 


Wx  =  q^*/22,klh 

g-moles/sec 

(5) 

Wf  =  Qf/Vj. 

g-moles/sec 

(6) 

wf  and  Wx  are  the  molar  flow  rates  of  the  fuel  and  oxidizer  and 
Qj  and  Qx*  are  the  volumetric  flow  rates  in  cc/sec  (the  latter 
at  S.T.P.).  Vf  is  the  molar  volume  of  the  fuel  and  must  be 
based  on  the  same  fuel  density  as  q^. 

v^  =  62,366(tnj  +  273.15)/Pni  cc/g-mole  (7) 

where  Vjj^  Is  the  molar  volume  of  the  mixture  and  tj^^  (°C)  and 
Pmv  (atm.)  are  the  temperature  and  pressure  of  the  mixture. 

Therefore:  q^^^  =  (w^  (8) 
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(b)  Area  of  flame  envelope,  A 


The  Increment  of  surface  area  between  any  two 

heights  hj^  and  h^_^  is  given  by 


^1-1^ 


+ 


'\-l  -^i 


cm 


where  is  the  measured  diameter  at  the  upper  height  and 

is  that  at  the  lower  height.  The  total  surface  area  of 
the  flame  is  calculated  by  summing  these  Increments  over  the 
whole  height 


(9) 


n 

A  =  ^  AA 
i=l 


_.2 


cm 


(10) 


(c)  Flame  speed,  u 

Once  the  total  flame  area  has  been  determined,  flame 
speed  can  be  easily  calculated: 

u  =  q^A  (11) 


(d)  Equivalence  ratio,  0 

Prom  its  definition  as  the  actual  fuel-to-oxygen 
ratio  divided  by  the  stoichiometric  fuel-to-oxygen  ratio,  the 
equivalence  ratio  is  calculated: 

0  =  (Wf./wQ)/(wj./w^g^)  (12) 

where  WqIs  the  molar  flow  rate  of  oxygen  used  in  the  experi¬ 
ment  and  Wqq^  Is  that  rate  required  for  complete  reaction. 

But  w.  =  w  z  g-moles/sec  (13) 

and  w^g^  =  Wfrg^/2  g-moles/sec  (l4) 

where  z  is  the  mole  fraction  of  oxygen  in  the  oxidant  and 
where  rst  is  ih®  number  of  atoms  of  oxygen  required  to  com¬ 
pletely  oxidize  one  molecule  of  fuel. 

Therefore  0  =  Wj.r  (15) 
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Pimctlon  2  (Subroutine  MAXM) 


Given:  A  series  of  experimental  values  of  flame 
speed  vs.  equivalence  ratio.  (The  flow 
rates  are  the  only  controllable  variables . ) 

Calculate:  (a)  Maxlm;im  flame  speed 

(b)  Flame  speed  at  stoichiometric  condi¬ 
tions 

(c)  Equivalence  ratio  at  maximum  flame 
speed 

(a)  Maximum  flame  speed,  Ujjjqx 


Past  work  has  shown  that  curves  of  u  vs .  0  generally 
are  concave  from  below  with  a  maximum  in  the  vicinity  of  0=1. 
Figure  3  is  an  example  of  such  a  curve . 

Many  of  the  curves  are  fairly  symmetrical  in  the  vicinity  of 
the  peak  and  are  nearly  parabolic  in  shape.  Others  are  fairly 
asymmetrical  and  are  better  approximated  by  a  higher  order 
power  series .  To  find  the  peak  of  a  given  series  of  measured 
values  of  u  vs .  0,  the  data  are  fitted  with  an  equation  of  the 
form 

R  , 

u  =  a  cm/sec  '(}6) 

r=l  ^ 

The  method  of  calculating  the  "a"  coefficients  depends  upon  the 
n\amber  of  points  available,  N.  If  N  =  4,  a  perfect  cubic  fit 
(R  =  4)  is  obtained  by  solution  of  the  simultaneous  equations 

R  T,  T 

If  N  >4,  a  least  square  cubic  fit  (R  =  4)  is  obtained  by  the 
solution  of  the  "normal  regression  equations." 


N 


n=l 


u„0„ 
n  n 


r-1 


R 

q=i 


(a 


N 

^  nfl 


0^^  +  <1-2)  r  =  1  ...R  (18) 


The  normal  regression  equations  are  the  result  of  a  mathematical 
treatment  on  the  original  polynomial.  Equation  17.  An  explana¬ 
tion  of  the  analysis  may  be  found  in  statistics  texts  such  as 
Hoel  (Ref.  6). 

Solution  of  these  equations  gives  the  set  of  "a"  coefficients 
which  minimize  the  residual  sum  of  the  squares: 
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Pig. 


(19) 


I 


Q 


2 


where  represents  a  predicted  value. 
The  standard  deviation  of  u  about  u  is 


S  = 


cm/sec 


(20) 


If  the  calculated  standard  deviation  is  greater  than  the  one 
specified,  the  cubic  fit  is  declared  unsatisfactory  and  a  para¬ 
bolic  fit  using  Equation  l8  with  r=3  is  tried.  If  the  standard 
deviation  is  still  too  high,  the  data  set  is  ignored. 


If  only  three  points  are  available,  a  perfect  parabolic  fit, 
r=3,  is  obtained  in  the  manner  of  Equation  17.  Less  than  three 
points  are  not  considered ' since  they  cannot  define  a  curvewith 
a  maximum.  . 


(b)  Equivalence  ratios 

(i)  Cubic  maximum  Once  the  cubic  fitting 
procedure  has  been  satisfactorily  completed,  it  is  necessary  to 
determine  the  maximum  point  (%ax>*^max)  curve.  Two 

additional  factors  must  be  considered.  The  maximum  velocity 
must  occur  within  the  data  range  considered,  and  a  minimum 
velocity  must  not  occur  within  the  data  range.  If  either  of 
these  conditions  is  violated,  a  parabolic  fit  is  attempted. 


Double  differentiation  of  Equation  17  with  r=4  gives 


u  =  a^,  +  2a^0  +  3ai^0^ 

(21) 

and 

u  =  2a^  +  6ai^0 

(22) 

At  the  maximum 

A 

u  =0 

max 

(23) 

and 

u  <0 

max 

(24) 

Solving  Equations  21  -  24  yields 

^max  =(a3/3a4)(l  -  sj  ) 

(25) 

and 

0  >  +  28^  1-38284/32^ 

(26) 
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The  Inequality  shows  that  the  upper  sign  is  to  be  used  in 
Equation  26  when  a^  is  positive,  and  vice  versa. 

For  a  real  maximum  to  exist,  three  conditions  must  be  fulfilled: 


0  (27) 

33  0  (28) 

(real  roots)  (29) 

(li)  Cubic  minjjmmi  The  equivalence  ratio 
at  the  minimum  point  (I3^j^^,j2>min)  of  the  cubic  equation  must 
also  be  determined: 

<^min  =  -(a3/3a4)(l  +  N/l-Sa^V^)  (30) 


The  upper  sign  is  understood  to  hold  when  a3  is  positive,  and 
vice  versa.  This  equivalence  ratio  must  not  lie  within  the 
data  range. 


(ill)  Parabolic  maximum  When  a  satisfactory 
parabolic  fit  is  obtained,  is  given  by 


^max  “  -32/233 


(31) 


in  which  a3  must  be  negative. 

Function  3  (Executive  routine) 

Given:  (a)  Combustion  mixture  temperature  and 

pressure 

Ib)  Fuel  stoichiometric  oxygen  demand 
c)  Mole  fraction  of  O2  in  oxidant 
d)  Equivalence  ratio  at  maximum  flame  speed 

Calculate;  (a)  Fuel  concentration  at  maximum  flame  speed 
(b)  Fuel  concentration  at  stoichiometric 
conditions 


(a)  Fuel  concentration  at  (cf*)Tnfly  For  any 


value  of  the  equivalence  ratio  0,  the  concentration  is 

given  by 


°f  “  b2,366(t 


(6.0238  X  10^^) pm 


m 


+  273.15)  (1 


(32) 


molecules  cc 
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Therefore  evaluated  at  0  =  0^^^^^  (33) 

(b)  Fuel  concentration  at  0  =  1,  (cp)Bf.r.f>_ 

Similarly, 

(cj.)gtoc  evaluated  at  0  =  1  (3^) 


Function  4  (Executive  Routine) 


1 

1  Given: 

Results  listed  above 

Composition,  source,  conditions  and 
related  input  Information 

1  Action: 

f 

Stores  in  memory  only  acceptable  data 
groups 

1  Function  5  (Subroutine  Tape) 

Given :  ( a ) 

(b) 

Acceptable  groups  stored  in  memory  1 

and  on  master  tape  6 

Changes  to  be  made  in  the  data  already 
on  the  master  tape 

Action: 

Edit  the  existing  master  tape  data  i 

Add  new  data  groups  to  the  master  tape 

In  all  cases,  the  Input  data,  intermediate  values,  and  output 
values  are  printed  out  at  each  step. 


Function  Options  Pour  options  are  provided  to  permit 
portions  of  the  program  to  be  by-passed.  The  options  would 
be  exercised  if  portions  of  the  experimental  data  had  been 
processed  previously  or  if  outside  data  were  to  be  used. 

1.  Complete  calculation  (Functions  1  through  4) 

2.  Partial  calculation  (Functions  2  through  4) 

3.  No  calculation  (Function  4) 

4.  Tape  writing  (Function  5  -  follows  successful 

execution  of  the  other  options) 


FLAME  SPEED  REGRESSION  CALCULATIONS  (Routine  PSR) 

Functions  Performed  Routine  PSR  prepares  the  input 
data  tape  for  the  Esso  Regression  Analysis  subroutine.  Its 
operation  may  be  considered  in  three  stages; 
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Function  1 


I 

t 


Given;  (a)  All  data  groups  currently  stored  on 
master  tape  6 . 

(b)  A  list  of  data  group  accept/reject 
tests 

Action:  Selects  and  stores  up  to  300  data 

groups  passing  the  accept/reject  tests 

Seven  different  tests  are  available  for  accepting  or  rejecting 
data  groups  now  recorded  on  the  master  tape  for  the  regression 
analysis.  These  tests  may  be  employed  In  any  combination  de¬ 
sired,  including  omission  of  any  or  all .  A  limit  of  300 
acceptable  data  groups  has  been  imposed  for  each  regression 
run.  If  more  than  300  could  pass  the  tests,  only  the  first 
300  doing  so  are  accepted  for  regression. 

(a)  Data  group  serial  number  test  Those  data 
groups  whose  serial  numbers  (machine  assigned  in  order  of 
tape  position)  are  given  in  an  input  specification  list  are 
unacceptable . 

(b)  Fuel  class  test  Only  those  data  groups 
whose  fuel  class  ni:imbers  are  specified  in  an  input  list  are 
acceptable . 

(c)  Fuel  group  test  Only  those  data  groups 
whose  fuel  class-group  numbers  are  specified  in  an  input. list 
are  acceptable. 

(d)  Fuel  member  test  Only  those  data  groups 
whose  fuel  class-group-meraber  numbers  are  not  specified  in 
an  input  list  are  acceptable. 

(e)  Data  source  test  Only  those  data  groups 
whose  source  number  corresponds  to  an  input  source  number 
are  acceptable . 

(f)  Experimental  conditions  test  Only  those 
data  groups  whose  .experimental  conditions  number  corresponds 
to  an  input  conditions  niimber  are  acceptable. 

(g)  Contributor  count  tests  Data  groups  may  be 
acceptable  or  unacceptable  depending  on  the  presence  or 
absence  of  given  contributors  in  the  fuel  molecule .  In  addi¬ 
tion,  the  counts  of  any  specified  contributors  may  be  assumed 
zero  whether  they  are  zero  or  not  on  the  master  tape. 
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Function  2 


Given ! 

fa)  All  data  groups  selected  as  above 

(b)  Dependent  variable  numerator  and 
denominator  code  numbers 

(c)  A  list  of  prespecified  coefficients 
with  their  corresponding  contributor 
nxambers 

Calculate: 

Adjusted  and  scaled  dependent  variables 
for  regression 

(a)  Dependent  variable  calculation  The  basic 
regression  problem  Is  to  determine  the  b  and  b  coefficients 
In  the  equation 

y  =  +  ^Vu  (35) 


such  that  the  sum  of  the  squares  of  the  difference  between 
the  actual  and  predicted  values  is  minimized.  Routine  PSR 
provides  the  option  to  define  the  dependent  variable  y  pro¬ 
vided  the  Independent  variables  are  still  contributor  counts. 
Thus.,  If  regression  of  the  maximum  flame  speed  on  contributor 
concentrations  at  maximum  flame  speed  is  desired,  y  would  be 
defined 

y  =  (36) 


The  dependent  variable  y  may  be  any  quantity  expressible  as 
a  simple  ratio  of  any  two  of  the  following  quantities: 

Quantity  Identification  Number 


1.0 

^stoc 

u 

max 

(°f ^stoc 
(°f ^max 
^max 


1 

2 

3 

4 

5 

6 


(b)  Dependent  variable  ad.lustment  If  It  Is 
desired  to  prespecify  (and  therefore  not  calculate)  some  of  the 
b  coefficients  in  a  given  problem,  it  Is  necessary  to  adjust 


I 


1 
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the  dependent  variable  by  subtracting  the  contributions  asso¬ 
ciated  with  the  prespecified  coefficients.  Equation  37  indi¬ 
cates  the  procedure 

y'l-yi-'^Vis  <37) 

(c)  Dependent  variable  scaling  To  avoid  exceed¬ 
ing  the  numerical  limit  of  the  computer  (l03o)  during  the 
regression  calculations,  all  of  the  final  dependent  variables 
are  scaled  (l. e. ,  multiplied  by  a  constant  factor) 

Yi  =  (38) 

so  that  the  largest  Y.  In  a  given  problem  does  not  exceed 
1000 .  ^ 


Function  3 

Given:  ^a)  All  data  obtained  above 

(b)  Input  lists  of  overriding  regression 
control  data  (optional) 

Action:  For  each  valid  regression  problem  store  on 
tape  7: 

ia)  Output  page  heading 
b)  Regression  control  data 
c)  Independent  variable  values  (contributor 
counts ) 

(d)  Adjusted  and  scaled  dependent  variable 
values 

The  regression  analysis  routine,  a  modification  of  SHARE  E-R- 
MPR2,  requires  Input  data  for  the  control  of  calculation  and 
printing.  One  of  the  functions  of  routine  PSR  is  to  supply 
the  needed  data  In  the  proper  form.  Although  this  Is  done 
automatically  using  Internally  stored  values.  Input  control 
data  may  be  given  If  It  is  desired  to  override  any  of  these 
values . 

In  all  cases  the  important  Input  data  and  Intermediate  and 
output  values  are  printed  for  each  function  listed  above . 
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V.  THE  COMPUTER  PROGRAM 


FLAME  SPEED  DATA  REDUCTION  (Routine  FSC) 

Routine  FSC  consists  of  an  executive,  or  master  control  pro¬ 
gram  and  the  eight  subroutines  listed  below. 


1. 

EXPDl 

2. 

EXPD3 

3. 

MAXM 

4. 

TAPE 

5. 

GROUT 

6. 

INPUT 

7. 

VDECOM 

8. 

DECDCP 

Section  IV  contained  a  brief  outline  of  the  functions  of  the 
executive  routine  and  the  first  four  subroutines  listed  above. 
Here,  these  shall  be  discussed  in  more  detail. 

Executive  Routine  The  executive  routine  controls  the 
selection  of  five  of  the  eight  subroutines.  The  first  data 
card  in  each  group  contains  information  regarding  the  option 
to  be  performed,  the  use  of  EXPDl  or  EXPD3  and  the  use  of 
MAXM.  The  executive  routine  reads  this  card  and  sequences 
the  operations  accordingly.  Diagnostic  code  numbers  are 
present  in  each  of  the  subroutines  and  if  any  of  these  are 
exercised  the  executive  routine  will  cause  the  printout  of  a 
diagnostic.  If  subroutine  MAXM  indicates  that  the  data  pro¬ 
cessed  were  satisfactory,  the  executive  routine  will  store 
the  pertinent  values  in  memory  prior  to  their  addition  to  the 
master  tape.  Two  small  calculations,  (cf)g^  and  (cf)inax  are 
also  made.  It  has  been  attempted  to  indicate  each  step  in  the 
actual  program  with  "comment  cards",  and  for  convenience  a 
nomenclature  list  precedes  each  routine.  These  can  be  found 
in  Appendix  A. 

Subroutines  EXPDl  and  EXPD3  Both  subroutines  perform 
the  same  basic  calculation  of  total  cone  area,  flame  speed, 
and  equivalence  ratio .  However,  the  respective  data  outputs 
differ  slightly.  It  was  considered  less  confusing  if  each 
alternative  had  its  own  subroutine.  Which  subroutine  to  use 
is  determined  mainly  by  the  intended  use  of  the  data  after 
processing. 
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If  the  data  to  be  processed  are  to  be  considered  for  addition 
to  the  master  tape^  subroutine  EXPDl  should  be  used.  In  this 
subroutine,  one  experimental  run  number  pertains  to  the  entire 
data  group  and  the  Input  and  output  are  set  up  accordingly. 

If  a  set  of  control  data  are  run  sequentially,  subroutine 
EXPD3  should  be  used.  Each  calculation  Is  listed  sequentially 
in  the  output  to  avoid  confusion  of  the  output  with  that  from 
EXPDl.  Sample  printouts  from  both  subroutines  may  be  found 
in  Appendix  A. 

Subroutines  MAXM  and  GROUT  Subroutine  MAXM  analyzes 
the  data  to  see  If  either  a  valid  cubic  maximum  or  a  parabolic 
maxlmimi  occurs  within  the  data  range.  Subroutine  CROUT  solves 
the  equations  generated  for  the  regression  coefficients.  As 
the  logic  used  In  subroutine  MAXM  Is  more  involved  than  in  the 
others,  a  brief  stepwise  description  is  presented. 

1.  If  the  number  of  experiments  (or  points)  per 
data  group  is  1  or  2,  signal  that  the  data  are  ’’bad"  and  re¬ 
turn  to  the  executive  routine.  If  there  are  3  points,  go  to 
step  6;  otherwise  go  to  step  2. 

2.  Calculate  the  cubic  coefficients  either  from 
the  cubic  polynomial  (4  points)  or  from  the  "normal  regression 
equations"  (more  than  4  points),  as  indicated  in  Section  IV. 

If  the  calculation  is  successful,  go  to  step  3;  if  not,  to 
step  6. 


3.  Check  cubic  maxima  criteria  (coefficients  zero 
or  imaginary  square  root).  If  all  are  satisfied  go  to  step 
4j  if  not,  to  step  6. 

4.  Calculate  ^max  i^min*  Check  for  ^max  Inside 
data  range  and  0ijiin  outside  data  range  of  0.  If  so,  go  to 
step  5;  if  not,  to  step  6. 

✓X 

5.  Calculate  Ug^gg,  u^,  absolute  and  per  cent 

differences  between  u^  and  un,  and  standard  deviation  S.  If  S 
is  less  than  or  equal  to  a  specified  value,  signal  the  data 
"good"  and  return  to  the  executive  routine;  if  greater, 
signal  the  data  "bad"  and  return  to  the  executive  routine. 

6.  Calculate  the  parabolic  coefficients  from  the 
quadratic  polynomial  (3  point)  or  from  the  "normal  regression 
equations"  (more  than  3  points).  If  the  calculation  succeeds, 
go  to  step  7j  If  not,  to  step  9. 

7.  Check  for  negative  coefficient  ao  (curve  con¬ 
cave  downwards).  If  so,  calculate  0ijiax  go^to  step  8; 
if  not,  to  step  9. 
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8.  Check  for  0max  inside  0  data  range.  If  so,  go  to 
step  5;  if  not,  to  step  9- 

9.  Signal  the  data  "bad"  and  return  to  the  executive 

routine . 

Subroutine  TAPE  This  subroutine  controls  all  operations 
that  change  the  contents  of  the  master  tape.  The  first  portion 
of  the  routine  scans  the  tape  and  makes  the  alterations  to  the 
data  groups  that  were  specified  on  an  input  data  card.  The 
second  section  permits  the  alteration  or  addition  of  contributor 
code  numbers  and  names.  The  third  section  adds  to  the  master 
tape  the  new  data  groups  that  were  either  approved  by  MAXM 
or  entered  via  option  3.  The  final  section  is  used  only  for 
the  initial  makeup  of  the  master  tape.  Once  the  tape  contains 
some  data  the  other  portions  of  the  program  may  be  used  without 
a  sequencing  diagnostic. 

Subroutines  INPUT,  VDECOM,  and  DECDCP  These  subroutines 
facilitate  preparation  of  the  input  data  cards  by  permitting  a 
variable  width  format .  The  complete  card  format  is  described 
in  Section  VI. 


FLAME  SPEED  REGRESSION  (Routine  FSR) 

Routine  FSR  prepares  the  input  data  for  the  E-R-‘MPR2  Esso 
Multiple  Regression  subroutine.  It  does  not  alter  the  input 
data  (with  the  exception  of  scaling  the  dependent  variable), 
but  serves  only  for  selecting  the  appropriate  data  for  the 
regression  analysis.  In  addition,  control  data  are  supplied 
for  the  Esso  program. 

Accept/Re.lect  Tests  Seven  tests  are  provided  for 
accepting  or  rejecting  data  groups.  The  first  six  tests  are 
simple  yes/no  choices  and  were  explained  in  Section  IV.  These 
tests  are  listed  below. 


1. 

Data 

group,  serial  number  test 

2. 

Fuel 

class  test 

3. 

Fuel 

group  test 

4. 

Fuel 

member  test 

5. 

Data 

source  test 

6. 

Experimental  conditions  test 

Test  7f  the  contributor  count  test,  is  more  complicated  and 
will  be  explained  here  in  detail. 
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Data  groups  will  be  found  acceptable  or  unacceptable  depending 
upon  the  presence  or  absence  of  particular  contributors  in  the 
fuel  molecule.  In  addition,  the  counts  of  any  contributor  may 
be  set  equal  to  zero.  The  code  nijmbers  and  test  operation  are 
listed  below. 

zero  test  -  The  count  of  the  specified  contributor(s) 
is  set  to  zero 

one  test  -  The  count  of  this  contributor  may  be  either 
zero  or  positive. 

two  test  -  The  count  of  this  contributor  must  be  zero 
or  the  group  will  be  rejected. 

three  test-  The  count  of  this  contributor  must  be  posi¬ 
tive  or  the  group  will  be  rejected. 

four  to  -  Six  tests  are  available  to  select  data  groups 
nine  tests  on  the  basis  of  a  particular  contributor  con¬ 
tent.  The  acceptance  of  a  data  group  may  be 
conditional  upon  it  containing  (+)  and/or 
excluding  (-)  a  certain  number  of  a  list  of 
contributors .  An  example  should  clarify  this . 

Example 

The  data  card  types  are  given  an  alpha¬ 
betic  identification.  The  contributor 
count  tests  are  specified  on  a  "P"  type 
card  and  the  conditional  tests  on  a  "Q" 
type  card . 

(a)  Sample  P-type  Card  1,  0,  15.  1, 
16,  0,  27,  4,  29,  5,  32,  3,  33,  34, 

4,  37,  2,  40,  6,  61,  0. 

Interpretation 

Apply  zero  test  (o)  to  contributors 
number  l-l4,  16-26,  6l  and  up 
Apply  one  test  (l)  to  contributor 
number  15 

Apply  two  test  (2)  to  contributors 
nvunber  37-39 

Apply  three  test  (3)  to  contributor 
n\imber  32 

The  above  tests  are  specific;  the  fol¬ 
lowing  are  conditional: 
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Apply  four  test  (4)  to  contributors 
number  27 j  28,  34-36 
Apply  five  test  (5j  to  contributors 
number  29-31^  33 

Apply  six  test  (6)  to  contributors 
number  40,  4l 

(b)  Corresponding  Q-type  Card 
+3,  -2,  +1 
Interpretation 

Pour  test;  the  data  group  must  con¬ 
tain  (+)  at  least  3  contributors  from 
the  list  27,  28,  34-36. 

Five  test;  the  data  group  must  not 
contain  (-)  at  least  2  contributors 
from  the  list  29-31,  33. 

Six  test;  the  data  group  must  con¬ 
tain  either  contributor  4o  or  4l. 

NOTE;  Conditional  tests  7  to  9  were 
not  used  in  the  example . 

Dependent  Variable  Limitations  Up  to  10  different 
dependent  variables  may  be  regressed  in  a  single  computer 
run  (problem  numbers  1  to  lO).  All  problems  in  the  same 
run  share  a  common  list  of  regression  control  data  and 
operate  on  the  same  set  of  selected  data  groups . 


Regression  Control  Data  The  regression  analysis  sub- 
routlne7  a  modification  of  SHARE  E-R-MPR2,  requires  several 
pieces  of  information  for  the  control  of  calculation  and 
printing.  One  function  of  routine  PSR  is  to  supply  the 
needed  data  in  the  proper  form.  This  is  done  automatically 
using  Internally  stored  values.  Only  when  it  is  desired  to 
override  one  or  more  of  the  Internal  values  is  it  necessary 
to  supply  input  control  data.  Internally  stored  values  are 
given  in  the  following  list. 
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Item 

Dec .  1 
2 
3 

Int.  1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

Values 


Description 


Built-In 

■Value 


Tolerance 

P  value  for  entering  variable 

P  value  for  removing  variable 

Problem  number 

N\imber  of  variables 

Nianber  of  points  (data  groups) 

Weighting  factors  given 
Intermediate  steps  printed 
Raw  sums  of  squares  and  cross  products 
printed 

Averages  to  be  printed 
Residual  siims  of  products  to  be  printed 
Partial  correlation  coefficients  to  be 
printed 

Predicted  values  (Y^'s)  to  be  printed 
Constant  term  (bo)  to  be  assiamed  non- zero 


0.001* 

0.00002* 

0.00001* 

calculated 

calculated 
calculated 
1  (no) 

0  (yes)* 

1  (no)* 

1  (no)* 

1  (no)* 

1  (no)* 

0  (yes)* 

1  (no)* 


marked  with  an  asterisk  may  be  overridden  by  Input  values. 


? 
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VI.  USE  OF  THE  COMPUTER  PROGRAMS 


Both  routines  have  been  written  in  the  standard  FORTRAN  II 
language  compatible  with  most  large-scale  computing  equipment. 

A  maximum  of  four  magnetic  tape  units  are  required  in  addition 
to  any  tape  units  used  for  input-output  operations .  The 
routines,  as  presented  in  Appendices  A  and  B,  are  programmed 
to  be  compatible  with  the  IBM  709O  Monitor  system,  with  Tape  2 
designated  as  the  Input  and  Tape  3  as  the  output. 

Two  additional  routines  have  been  Included  to  avoid  the  storage 
of  the  master  data  tape  over  a  period  of  time.  Routine  PSRDM 
will  prepare  a  master  data  deck,  or  modify  an  existing  data 
deck  from  Information  contained  on  the  master  data  tape. 

Routine  PSRTL  will  prepare  a  tape  from  the  master  data  deck. 
These  routines  are  presented  in  Appendix  C. 


INPUT  DATA  DECK  PREPARATION 


The  programs  accept  the  standard  IBM  type,  80-column  punched 
cards .  However,  only  coliamns  1  to  72  may  contain  information 
to  be  processed  by  the  computer.  It  is  recommended  that 
columns  73-80  be  used  for  some  type  of  identification  both 
for  the  users  standpoint  as  well  as  for  the  computing  center. 

All  decimal  Input  data  and  most  Integer  input  data  are  pro¬ 
cessed  via  subroutines  VDECOM  and  DECDCP  before  being  used  In 
the  computation  portion  of  the  main  routine.  The  use  of  these 
subroutines  simplifies  data  card  preparation  to  the  extent 
that  field  widths  may  be  neglected.  It  Is  only  necessary  to 
provide  a  blank  space  or  a  comma  between  entries  on  the  data 
cards.  One  restriction  Is  imposed  though.  No  entry  may  end 
In  column  72  on  the  data  card. 


FLAME  SPEED  DATA  REDUCTION 


The  Input  data  deck  may  consist  of  several  sets  of  cards 
depending  on  the  nature  of  the  data  and  the  option  applicable. 
If  only  the  calculated  flame  speed  is  of  Interest,  an  unlimit¬ 
ed  number  of  data  sets  may  be  entered.  If,  however,  options 
1,  2,  or  3  are  used  for  the  purpose  of  adding  data  groups  to 
the  master  tape,  an  option  4  must  follow  every  20  sets  or 
less.  This  limitation  is  imposed  because  of  the  Internal 
storage  requirements  of  the  program.  Options  1,  2,  or  3  may 
be  run  in  any  order. 
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To  simplify  the  makeup  of  the  data  decks,  each  card  fulfilling 
a  specific  purpose  has  been  assigned  an  alphabetic  code  letter 
Below  are  listed  the  data  deck  requirements  for  each  of  the 
options  available .  The  actual  makeup  of  the  cards  follows . 


Option  1  Calculation  of  Flame  Speeds  from  Raw  Data 


The  first  four  card  types  set  up  the  calculation  sequence  and 
define  certain  general  conditions .  The  remaining  cards  con¬ 
tain  the  experimental  measurements  for  each  flame  photograph. 
Thus,  for  a  total  of  Ng  photographs,  the  card  order  would  be: 


Card  Order 


Card  Type 


1  A 

2  X 

3  B 

4  C 

5  to  N2  +  4  D 


Certain  entries  on  the  A  card  determine  whether  or  not  a  full 
calculation  or  Just  the  flame  speeds  are  required .  If  the 
latter  Is  the  case,  any  number  of  these  sets  may  be  submitted, 
otherwise  the  maximum  number  before  a  tape  writing  sequence  is 
20.  Ng  Is  limited  to  10. 


Option  2  Calculations  for  Given  Equivalence  Ratios  and 
Flame  Speeds 


The  first  four  card  types  are  similar  to  those  mentioned  above 
The  last  card  contains  the  predetermined  flame  speeds  and 
equivalence  ratios . 

Card  Order  Card  Type 


1 

2 

3 

4 

5 


A 

X 

B 

C 

E 


Option  3  All  Final  Values  Given 

Only  the  Information  on  the  first  three  card  types  Is  required 


Card  Order  Card  Type 

A 
X 
B 


1 

2 

3 
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Option  4  Tape  Writing  and  Editing 

The  first  card  type  remains  the  same.  The  second  card  tsrpe  con¬ 
tains  Information  about  the  remaining  two  card  types. 


Card  Order 

1 

2 

3 

4 


Card  Type 

A 

F 

G 

H 


Card  Content  by  Card  Type  All  numerical  Input  data  for 
routine  FSC  are  processed  by  subroutines  VDECOM  and  DECDCP.  In 
the  listings  belowj  all  Integer  quantities  are  designated  Int. 
and  all  floating  point  quantities  (e.g.  ,002  or  2.E-3  or  2.-3) 
are  designated  Dec.  The  only  restrictions  on  preparation  of 
these  cards  are  that  column  72  must  be  blank  and  each  entry 
must  be  separated  by  at  least  one  blank  or  comma.  Readln  of 
alphanumeric  Input  Is  accomplished  with  an  "A"  format  and  the 
card  col\imn  designations  below  must  be  maintained. 


A  Card 
Entry 


1 

2 

3 

4 

5 

6 

7 

8 


X  Card 


Entry 


Int.  1000,  card  Identification 

Int.  Run  number  (l  to  9999) 

Int.  4,  number  of  Integer  entries 
following 

Int.  Option  number  (l  to  4) 

Int.  EXPD,  subroutine  number  (l  or  3) 
Int.  No,  number  of  Inhibitors  present 
Int.  MAXM,  subroutine  switch  (O-use, 
1-bypass) 

Int.  1,  number  of  decimal  entries 
following 


9 

Dec . 

0 . ,  dummy  entry 

Columns 

1 

Int. 

0000,  card  ldentlfl< 

2 

Int . 

0,  dummy  entry 

3 

Int . 

0,  dummy  entry 

4 

13-24 

Alphabetic,  fuel  name 

5 

25-36 

Alphabetic,  date 

6 

37-48 

Alphabetic,  Inhibitor  name 

7 

49-60 

Alphabetic,  Inhibitor  name 
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A  maximum  of  ten  inhibitor  names  may  be  entered .  These  names  ' 
must  follow  sequentially  on  the  data  cards  and  must  occupy 
twelve  columns  per  name.  Columns  1  to  12  on  each  card,  however, 
must  contain  the  three  Integer  entries  listed  above.  Thus, 

If  No  were  4,  the  fuel  name,  data  and  three  inhibitor  names 
would  go  on  the  first  X  card  while  the  fourth  name  would 
appear  in  Columns  13  to  25  of  the  second  X  card. 


B  Card 

Entry  1 

Int . 

2000,  card  identification 

2 

Int . 

Run  number 

3 

Int . 

7,  number  of  integer  entries 
following 

4 

Int . 

Option  number  (l  to  3) 

5 

Int. 

Fuel  class  number  fo  to  99) 

6 

Int . 

Fuel  group  number  (0  to  99) 

7 

Int . 

Fuel  member  number  fo  to  99) 

8 

Int . 

Data  source  number  (1  to  9999) 

9 

Int . 

Experimental  condition  number 
(1  to  9999) 

10 

Int . 

Ni,  number  of  contributors  listed 
on  following  B  cards  (l  to  100 ) 

11 

Int . 

5,  number  of  decimal  entries 
following 

12 

Dec . 

Flame  speed  at  stoichiometric 
conditions  (cm/sec)* 

13 

Dec . 

Maximum  flame  speed  (cm/sec)* 

l4 

Dec . 

Fuel  concentration  at  stoichiometric 
conditions  (molecules/cc )* 

15 

Dec . 

Fuel  concentration  at  maximum  flame 
speed  (molecules/cc)* 

l6 

Dec . 

Equivalence  ratio  at  maximum  flame 
speed* 

*Enter  0.0  when  these  values  are  not  known  (l.e.,  for  options 
1  and  2) 


B  Card  Continuation 


Entry  1 

Int . 

2001,  card  identification 

2 

Int . 

Run  number 

3 

Int . 

N2,  number  of  integer  entries 
following  (N2  -^10 ) 

4  to  N2+3 

Int . 

Contributor  code  niunbers 

N2+4 

Int. 

N2,  number  of  decimal  entries 
following 

N2+5  to 
2N2+5 

Dec . 

Contributor  count 
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If  more  than  ten  contributors  are  present,  additional  continua¬ 
tion  cards  may  be  used.  For  each  subsequent  continuation  card 
increase  the  card  identification  number  by  one  (l). 

C  Card 

Entry  1  Int.  3000,  card  identification 

2  Int .  Run  number 

3  Int.  2,  number  of  Integer  entries  following 

4  Int.  Option  number  (l  or  2) 

5  Int.  N^,  number  of  experiments  in  data 

group  (EXPDI)  or  number  of  experimental 
runs(EXPD3) (i.e . ,  number  of  D  cards) 

6  Int.  7,  number  of  decimal  entries  following 

7  Dec.  Mixture  temperature  (°C) 

8  Dec .  Mixture  absolute  pressure  (mm  Hg) 

9  Dec.  Stoichiometric  oxygen  ratio  (atoms 

oxygen/molecule  fuel) 

10  Dec.  Mole  fraction  O2  in  oxidant 

11  Dec.  Volume  per  mole  of  fuel  at  given  fuel 

flow  rates  (c c/gram-mo le) 

12  Dec.  Maximum  allowable  standard  deviation 

in  u  vs .  0  curve  for  addition  to 
tape  (assumed  1.0  cm/sec  if  0.0 
entered) 

13  Dec.  Actual  distance  between  teeth  tips 

seen  on  schlieren  photograph  (cm) 
(assumed  to  be  0.2  cm  if  0.0 
entered) 

D  Card 

Entry  1  Int.  4000,  card  identification 

2  Int.  Run  nvmiber 

3  Int.  2,  number  of  integer  entries  following 

4  Int.  1,  option  number 

5  Int.  N4,  number  of  diameter  measurements 

listed  on  continuation  card  (l  to  lOO) 

6  Int.  6,  number  of  decimal  entries  following 

7  Dec.  Fuel  flow  at  mole  volume  of  entry  11 

on  C  card 

8  Dec.  Oxidant  flow  at  STP  (cc/sec) 

9  Dec .  Mole  fraction  of  inhibitor  in  final 

mixture 

10  Dec.  Measured  peak  height  (arbitrary  units) 

11  Dec .  Reference  length  Tern) 

12  Dec.  Reference  length  (arbitrary  units) 
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D  Card  Continuation 


Entry  1  Int . 

2  Int.* 

3  Int. 

4  to  N5  +  3  Int.** 
N5  +  5  to  2N5  +  5 


4001,  card  identification 
Run  number 

Nc,  niomber  of  integer  entries 
following 

Station  height  of  diameter  measure¬ 
ment  (1  =  cone  base) 

Cone  diameter  (arbitrary  units) 


*Por  data  processing  by  subroutine  EXPDl,  the  run  number  on 
each  D  card  must  be  the  same  as  that  on  the  A,  B  and  C 
cards.  For  processing  by  routine  EXPD3,  the  run  number  is 
to  be  incremented  by  one  for  each  subsequent  D  card. 


**The  station  heights  correspond  to  the  "teeth  locations 
referred  to  in  the  C  card  description.  The  distance 
between  twenty  reference  teeth  was  accurately  known.  As 
all  flame  speed  measurements  were  made  on  a  microfilm 
reader,  it  was  found  to  be  more  accurate  and  more  convenient 
to  measure  the  flame  diameters  and  teeth  spacing  in  any 
suitable  \anlts  (say  mm)  and  have  the  computer  rescale  the 
measurements  to  the  desired  units . 


E  Card 

Entry  1  Int . 

2  Int . 

3  Int. 

4  Int. 

5  Int. 


E  Card  Continuation 


5000,  card  identification 
Run  number 

1,  nijmber  of  integer  entries 
following 

2,  option  number 

No,  number  of  flame  speed  - 
^equivalence  ratio  entries 
(1  -  10) 


Entry  1 
2 
3 


Int . 
Int , 
Int. 


4  Int. 

5  Int.* 


6  to  6  +  Dec . 


5001, card  identification 
Run  number 

1,  number  of  integer  entries 
following 
0,  dummy  entry 
N^,  n-umber  of  flame  speed 
equivalence  ratio  entries 


*No  is  entry  5  on  the  C  card  and  must  not  have  a  value 
exceeding  10. 
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P  Card 


Entry 

1 

Int . 

6000,  card  identification 

2 

Int . 

Run  number 

3 

Int . 

5,  number  of  integer  entries  following 

4 

Int . 

4,  option  number 

5 

Int. 

Ng,  number  of  G-card  groupings  to  follow 
(O  to  20) (contributor  count  changes) 

6 

Int . 

Nj,  number  of  entry  pairs  on  H-card 
(6  to  50) (contributor  name  changes) 

7 

Int . 

7,  for  initial  tape  makeup,  otherwise 
zero 

8 

Int . 

0,  for  other  than  last  data  card  set,* 

1,  for  last  set 

9 

Int. 

1,  ntunber  of  decimal  entries  to  follow 

10 

Dec . 

0. ,  dummy  entry 

G  Card 

Entry 

1 

Int . 

7000,  card  identification 

2 

Int . 

Run  number 

3 

Int . 

2,  number  of  integer  entries  following 

4 

Int . 

Data  group  serial  number 

5 

Int. 

Ng,  number  of  entry  pairs  in  contributor 
count  change  list  (O  to  90) 

6 

Int. 

1,  number  of  decimal  entries  following 

7 

Dec . 

0., dummy  entry 

G  Card  Continuation 

Entry 

1 

Int . 

7001,  card  identification 

2 

Int. 

Run  number 

3 

Int.* 

Ng,  number  of  Integer  entries  following 
(0  to  10) 

4  to 

Nq+3 

Int . 

Contributor  number 

Ng 

+4 

Int. 

Ng,  nvunber  of  decimal  entries  following 

Ng+5  to 

4+2Ng 

Dec . 

Contributor  count 

*If  N3>10  additional  continuation  cards  will  be  required.  The 
card  identification  number  must  be  incremented  by  one  on  each 
continuation  card. 
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H  Card 


Columns  1  to  4 
6 

r  8 

13  -  16 
25  -  28 


0000 

0 

0 

contributor  number 
contributor  number 


61 

19 

31 


64 

24 

36 


contributor  number 
alphabetic  contributor  name 
alphabetic  contributor  name 


67  -  72  alphabetic  contributor  name 


*There  are  to  be  five  entry  pairs  per  H  card.  If  N„  from  the 
P  card  is  greater  than  5^  continue  the  listing  on  additional 
cards  at  five  pairs  per  card. 


Fortran  Routine  The  Fortran  II  program  is  listed  in 
Appendix  A  with  sample  printouts.  Duplicate  decks  of  this 
program  are  available . 


FLAME  SPEED  REGRESSION 


The  niomber  and  type  of  cards  comprising  the  data  deck  will 
depend  upon  the  number  of  tests  and  changes  wanted.  To  facili¬ 
tate  the  deck  makeup,  each  specific  data  card  has  been  given  an 
alphabetic  code  letter  from  I  through  U.  Card  types  I,  J,  K, 
and  T  in  the  deck  are  mandatory  with  the  K-type  card  containing 
control  integers  governing  the  number  and  type  of  all  but  one 
of  the  remaining  cards . 


Card  Order 


Card  Type 


1 

2 

3 

Next,  unless  N2 
Next,  unless  N^ 
Next,  unless  N^ 
Next,  unless  Nc 
Next,  unless  Ng 


I 

J 

K 

=  0  L 

=  0  M 

=  0  N 

=  0  6 

=  0  P 
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Card  Order 


Next,  unless  Ny  =  0 

Next,  unless  Ng  =  0 

Next,  unless  Ng  =  0 

Next 

Next,  unless  N^q=  0 


Card  Type 


Q 

R 


*N]^  pairs  of  these  cards  may  be  used  (N^^^IO).  The  T  card 
contains  N^q  and  hence  governs  the  use  of  the  U  type  card. 

Card  Content  by  Card  Type  With  exception  of  the  I,  R  and 
U  type  cards,  all  Input  data  to  routine  PSR  are  Integers.  A 
fixed  format  is  used  to  accomplish  this  with  four  spaces  allotted 
to  each  entry  (l8l4).  The  type  U  cards  are  decomposed  in  sub¬ 
routines  ‘VDECOM  and  DECDCP  with  the  same  limitations  as  described 
for  routine  ESC.  When  card  column  designations  are  given,  they 
must  be  adhered  to. 


I  Card  -  Subtitle  Information  for  Routine  FSR  Printout 


Colijmn  1  Alpha  Blank 

2-72  Alpha  Subtitle 

J  Card  -  Title  for  Regression  Program  Printout 

Column  1  Alpha  One  (l) 

Column  2-72  Alpha  Title 


K  Card  -  Routine  FSR  Control  Card 


Entry  1 
2 

3 

4 

5 

6 
7 


Int .  Regression  run  number  (l  to  9999) 
Int.  Nt,  number  of  dependent  variable 
choices  in  this  run 

Int.  Np,  number  of  entries  in  input  list 
(l  card)  of  unacceptable  data  group 
serial  numbers  (O  to  300 ) 

Int.  No,  number  of  entries  in  input  list 
(M  card)  of  acceptable  fuel  class 
numbers  ( )  to  13 ) 

Int.  N4,  number  of  entries  in  input  list 
(N  card)  of  acceptable  fuel  class 
group  numbers  (O  to  4o) 

Int.  Nc,  number  of  entries  in  input  list 
(0  card)  of  unacceptable  fuel  class- 
group-member  nvunbers  (O  to  80) 

Int.  Acceptable  data  source  number  (O  to 
999»  0  means  all  sources  acceptable). 
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Entr; 


L  Card 
Entry 

M  Card 
Entry 

N  Card 
Entry 

0  Card 
Entry 

P  Card 
Entry 

Q  Card 
Entry 


8 

Int . 

Acceptable  experimental  conditions 
number  (O  to  999 t  0  means  all 
conditions  acceptable) 

9 

Int . 

Nc,  number  of  entry  pairs  In  Input 
list  (P  card)  of  contributor  count 
tests  (0  to  200) 

10 

Int. 

number  of  entries  In  Input  list 
((i  card)  of  conditional  count  test 
criteria  (O  to  6) 

11 

Int . 

N3,  one  If  overriding  decimal  re¬ 
gression  control  data  (R  card)  to  be 
read;  otherwise  zero 

12 

Int . 

Nq,  number  of  entry  pairs  In  Input 
list  (S  card)  of  overriding  Integer 
regression  control  data  (O  to  7) 

13 

Int . 

Highest  contributor  number  In  use 
(1  to  200,  higher  ones  Ignored) 

14 

Int . 

1,  for  the  last  run,  otherwise  zero 

to  Ng 

Int . 

Unacceptable  data  group  serial  numbers 
(N2  ^^300) 

to 

Int . 

Acceptable  fuel  class  numbers  (N^^20) 

to  2N2^ 

Int . 

Acceptable  fuel  class-group  numbers 
(N2|  £.40}  paired. 

to  3N^ 

Int . 

Unacceptable  fuel  class-group-members 
(N^^8o)  In  trios. 

:e  "Contributor 

Count  Tests") 

3  2N5-I 

Int .. 

Contributor  n\imber  N5  pairs 

4  2Ng 

Int. 

Contributor  count  test  (Ng'£.200) 

to  Ny 

Int . 

Conditional  count  test  criteria. 

First  entry  goes  with  test  4,  second 
with  5,  etc.  (Nj^S) 
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R  Card  (see  "Regression  Control  Data") 

Columns  1  to  10  Dec .  Tolerance 

11  to  20  Dec .  P  value  for  entering  variable 

21  to  30  Dec.  P  value  for  removing  variable 

^  Card  (see  "Regression  Control  Data") 

Entry  1,3  2Nq-1  Int.  Integer  regression  control  Item 

^  number  N  ^  7 

2,4  2Ng  Int.  Overriding  Integer  regression  ”9 

^  control  data 

T  Card 

Entry  1  Int .  Dependent  variable  nimierator 

Identification  number  (l  to  6) 

2  Int .  Dependent  variable  denominator 

Identification  number  (l  to  6) 

3  Int.  N,q,  niiiTiber  of  entry  pairs  In  Input 

list  (U  card)  of  prespecified  re¬ 
gression  coefficients  for  this 
problem  (O  to  50 ) 

U  Card 

Entry  1  Int.  1,  card  number 

2  Int .  Run  number 

3  Int.*  N-j^n,  number  of  Integer  and  decimal 

pairs  to  follow 

4  to  3+N-i  -1  Int .  Contributor  number 

+4  Int. 

Nii+5  to  4+2N^^  Dec. .  Prespecified  coefficient 

*If  N]^q>10  additional  U  type  cards  will  be  required  until 
all  of  the  contributor  numbers  and  coefficients  are  entered. 

The  card  number  must  be  Incremented  by  one  for  each  additional 
card . 

Fortran  Routine  The  Fortran  II  program  for  routine  PSR 
Is  listed  In  Appendix  B.  Duplicate  decks  of  this  program  are 
available . 


MASTER  DATA  DECK  LOADING  AND  MODIFICATION 

Routine  FSRTL  The  only  data  cards  required  by  the 
routine  are  those  of  the  master  data  deck  Itself.  The  card 
Image  of  this  deck  will  be  read  onto  logical  tape  6  which 
will  then  be  given  an  end  of  file  and  rewound. 
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Routine  FSRDM  If  It  is  desired  to  update  the  master 
data  deck,  routine  FSRDM  will  cause  any  or  all  of  the  data 
groups  on  the  master  tape  to  be  punched.  The  master  tape  must 
be  logical  tape  6 .  All  control  data  for  the  routine  are  Integers 
and  a  fixed  format  Is  used  which  provides  a  slx-coliimn  field 
width  (1216). 

First  control  card  for  Routine  FSRDM 


Entry 

1 

Int . 

Number  of  data  groups  in  the  master 
data  card  deck 

2 

Int . 

N^,  number  of  data  groups  which  have 
been  altered  since  last  deck  makeup 

3 

Int . 

Punch  contributor  name  list? 

( 1-yes,  2-no) 

h 

Int . 

Punch  dependent  variable  name  list? 

( 1-yes,  2-no) 

If  none  of  the  data  groups  have  been  modified,  =  0  and  the 
routine  will  punch  only  those  data  groups  added  since  the  last 
makeup.  If  x  0  an  additional  control  card  is  required. 

Second  control  card  for  Routine  FSRDM 

Entry  1  to  Int.  Serial  number  of  data  groups  that 

have  been  modified  since  last  deck 
makeup . 

The  Fortran  II  programs  for  routines  PSRTL  and  FSRDM  are  listed 
In  Appendix  C.  Duplicate  decks  of  these^programs  are  available. 
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VII.  NOMENCLATURE 


a  Power  series  fitted  coefficient  p 

A  Area  of  flame  schlieren  cone,  cni'^ 

b  Regression  coefficient 

c  Number  concentration,  cm“^ 

D  Diameter  of  flame  schlieren  cone,  cm 

E  Dependent  variable  scaling  exponent 

h  Height  above  burner  port,  cm 

n  Contributor  count,  no. /molecule 

N  A  number 

p  Absolute  pressure,  mm  Hg 

q  Volumetric  flow  rate,  cm^/sec . 

Q  Residual  sum  of  squared  differences 

r  Atoms  of  oxygen  per  molecule  of  fuel 

R  Maximum  value  of  r  subscript 

S  Standard  deviation 

t  Temperature,  °C 

u  Flame  speed,  cm/sec. 

V  Molar  volume,  cm3/g-mole 

w  Molar  flow  rate,  g-moles/sec . 

X  Independent  variable 

y  Dependent  variable 

y  Adjusted  and  scaled  dependent  variable 

2  Mole  fraction 

A  Increment  prefix 

Tt  3.14159  ••• 

0  Equivalence  ratio 


Subscripts 


f  Of  fuel 

1  Data  group  index,  block  regression 

J  Contributor  index  (or  number) 

m  Of  combustible  mixture 

min  At  minimum  of  cubic  fitting  equation 
max  At  maxlmimi  flame  speed 

n  Height  index 

o  Of  oxygen,  or  zero  index 

p  At  peak  of  flame  cone 

q  Dummy  r  index 

r  Term  index,  u-0  curve 

s  Specified  (or  prespecified) 

stoc  Stoichiometric 

u  Unspecified  (to  be  calculated  by  regression) 

X  In  oxidant 
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<  * 


Superscripts 

Calculated  from  a  fitted  equation 
At  STP  (0°C,  760  mm  Hg  abs.) 
Adjusted 

First  derivative 
Second  derivative 
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APPENDIX  A 

Fortran  Program  and  Sample  Printouts 
for  Routine  PSC 


HOUTINC  FSCK 


"‘ibia  Ffa'io  tiHi  ‘  speeo  calculations  '  moo 'i 

eitECUIlVC  «OufiNYT6A“FLVHe‘'$PBeD  CALCULA"tIONS  HRC  -  OAVTON 


NONCNCLAfURI; 

“  Al'i  rrANO  Aim  -  FUEL  NAME 

_ -  NUMAER  OF  STRUCTURAL  CONTRIMUTORS  FEE  MOLECULE  _ 

CNX  -  FUEL  CONCENTRATION  AT  NAAINUM  FLAME  VELOCITY 

_ _ - _ 

CSfOC  -  FUEL  CONCENTRATION  AT  STOICHIOMETRIC 

_ £fWO.LT.10NS.„|NOLE.C.UlES  ./_Cjt.J.  _ 

OU  AND  TE  -  FIRST  AND  SECONO  HALVES  OF  DATE 

.  I.  R.  AMft  N  -  COUNTING  INTEGERS . . . . 

ILAdl  -  FUEL  CLASS  NUMBER 


OU  ANO  TE 
It  K._ANO_K_ 
ILAdl 

-llBdl 

llCdl 

_UJ.I.t _ 

UdI 

-JAtU _ 

lAOO 

IGTL 

..lEXFO _ 

iNHlBlUtll  I 

IHMlQgtJdl 

IPN 


•  FUEL  MEMBER  NUMBER 

«  DATA  SOURLc  .NUMBER _ 

-  EXPERIMENTAL  CONDITIONS  NUMBER 

~^I^Blft'OF  DA'tA^CRuW  acceptable  FOR  WRITING 
_=jmQ«LA01lTJIPLJJ«TXG.EJLiH^^ _ 

-  OPTION  CONTROL  INTEGER 

_r_EXPD  .  SUDROUTlNE.SELeCJ.IflN_IKT.E.ft£B _ 

_=_lNHlft.|.lO|l.NAME5..  _ _ _ _ 

NUMBER  OF  EXPERIMENTS 

-  SPECIES  CONTRiaiitnft  CODE  number 

•>  SIGNAL  INTEGER  FOR  ACCEPTABLE  DATA 
_=_RUN.  NUMBER  OF  .DATA  CROUP _ 

-  RUN  NUMBER  CHECK 

.-‘^_£<H»I VALENCE  ;Ul|ns  _ 

>  NUMBLA  OF  iNHIttilURS  USED 
_«  AIQMS  OP  DUVUCM  in  CftMPLgfFLV  OXIPllE  ONE  " 
MOLECULE  OF  FUEL 

-r_£QUlVALENCE  RATIO  AT  MAXINUH.^LAHE-yELOCliy _ ..... 

-  ABSOLUTE  PRESSURE  OF  GAS  MIXTURE  INH.  MERCURVI 

_=_MAXIMUM  ALLOHABLE  STANOARD..OE.VIALU»-ttE._U..VS _ 

0  CURVE  FOR  ADDITION  TO  TAPE 
»  EQUIVALENCE  RATIU _ 

•  FLAME  SPEED  (CM. /SEC) 

--^Cr.WAL..OiSTANCe  BEfWEENJTEETH  TIPS_SEEN  ON _ 

SCHLIEREN  PHOTOGRAPH  (CM. I 
-JtJlUT.URE.  JEMPERATURE  , ..  IDE&i.Jt» _ _ 

•  MAXIMUM  FLAME  SPEED  (CM./SEC) 

(CM. /SEC) 

..^uJVOLUHE  PER  HOLE  OF  FUEL. _ (CC*/GftAit:NOLli _ 

•  VOLUME  PER  MOLE  OF  MIXTURE  ICC •/GRAM-MOLE  I 

»  MOLE  FRACTION  OXYGEN  IN.OXIDANt _ . 


COMMON  INT»  DEC*  1C«  JU»  J2t  JS»  N1N»  NiX 

OIMEMSION  INTdOl*  OH(LO)  . . . . 

DIMENSION  A142O)*A2(20)»Btll0O»2O)»DAI6»2OltllAI2O>»llSltOI» 

I _ ttCiifil#  12(20)  t.I.>.(  20 ) » I A  ( ^0 )  >V i. t.  100»20)tO(lBUUdO) 

I  »INMlBldO«20)«  lNHtB2d0»20) 


FSCX  7 

_FSCX _ B 

FSCX  9 

_FSCX _ 10 

FSCX  11 

FSCX _ 12 

FSCX  13 

...FSCX  14 

FSCX  IS 

_FSCX _ U. 

FSCX  IT 

_FSCX  .  _  18 
FSCX  19 

_FSCX  20 

FSCX  21 

_FSCX _  22. 

FSCX  23 

«_FSCX_..  24 
FSCX  25 

FSCX  .  ,26 
FSCX  27 

_FSCX _  28 

FSCX  29 


FSCX  35 

...FSCX  36 

FSCX  37 

FSCX  38 

FSCX  39 


U  ANO  UIN) 
_ TEETH _ 


DIMENSION 

_ I _ 

2 

_ .NIN..».| _ 

to  lAOO  •  0 
30  CALL  INPUt 
NEX  •  NEX 


_ A„CARO..AMJ.<a!ll!WlU_- 

31  NRUN.  ■  Jll  _  _ _ _ 

ICTL  •  INTd) 

__  _IIXP0  »  INT(21 _ 

NUMINH  ■  1NTI3I 

_ ^MAX(JN  Ji_UlTJiU _ 

IF  (NRUN)  32»210»32 
__32_l F  ( NUMINH )  320»320i33 
320  NUMINH  «!' 

OPTION  SELECTION 


33  IF  dCIL  •  4)  4A.tS0,30 

_40  IF  dCII  )  .twJo.vu  _  _ 

SO  I  >  lAOO  ♦  1 


_ FSCX _ 40 

FSCX  41 

...FSCX  ..  42 
FSCX  43 

—  FSCX  .  44 

FSCX  45 

.^FSCX _ 46. 

FSCX  47 

...FSCX  .  48 

FSCX  49 

..  FSCX  $0 

FSCX  SI 

—  FSCX  ..  .52 

FSCX  $3 

FSCX  54 

FSCX  55 

FSCX  56 

FSCX  '  57 

_  FSCX  .  58 

FSCX  59 

_  FSCX  6C 

■FSCX  61 

_ FSCX _ 62 

FSCX  63 

__  FSCX  64 

FSCX  65 

_  .  FSCX  66 

FSCX  67 

_ FSCX _ 6B 

FSCX  69 

_  FSCX  70 

FSCX  *  71 
_  FSCX_  72 
FSCX  “■  73 

_ FSCX _ 74 

FSCX  75" 

FSCX  _  76 
rscx  "77 
F5CX  70 

FSCX  ■  ‘  79 
•  FSCX _ 80 

FSCX  »r 


52  READ  INPUT  TAPE  2*  1002 

_ 9  »  Ald),A?()),DU.Tr  .(INIlinM  J,n,|NH|82(J|l 

INUNINH) 

_  (;0  TO  30  __  -  _ 

c 

c  „  6  CARO  ASSIONNfNTl _ _ _ 

C 

.  .  Jll  NUNCON IQ^_ 2.«)0 _ 

IF  INUMCON)  523»  522t  523 


ICANT  •  INIdl 

_IUd»  •,INT.|2| _ ^ _ 

Iieilt  •  INT<3) 

„.nClllj!  INT.(4| _ 

I2dt  •  INTIS) 

_ 1 3  L II  »  1NTL6J _ 

Kl  •  INTI?) 

00  5220  R  •  2(6 _ 

5220  OAlKfll  •  OECU'll 

CO  TO  30  _ 

521  IFf  NRUNT  -  211)  300t  5210*  >00 
_5210.L tifit  ■  ia»  ttUMCON 
IW  •  UNIT  -  9 

.  IFILINIT  -  Rll  325«Ji29»„52A _ 


FSCX  100 
FSCX  101 
FSCX  102 
FSCX  103 
^FSCX  104 
FSCX  105 
.  FSCX  106 
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ROUTINE  ESCl 


"flame  SREtO  ULCUtATlW  liifCuflVI  fkdUAil 


Kl 


524  LIMIT 
529  M  •  0 

_ DO  52A_K_«  ICO.  LiHIT _ 

M  •  M  ♦  I 

_  JURtU  -  IMTIMI _ 

524  AHKit)  •  OCCINI 

IF<  LIMIT  -  Rl»  _  9*1#. 

55  t4Ui  #  Kl 


C 

.C_ 


MRlNTOUr  OF  INAUT  DATA 

UAITE  QUTFUT  tTfFT*  1010 
H  •NAUN,OU»T€ 

UR1T£  OUTFUT  TAME  5«  1011 
_HA J  TE  QUTFUT  TAFt  5.  1012 


.  Aim.  A2(n  •I1NH|«1IJ,1)«  rNHIB2U«lltJ«U  MUNINH) 

NR|.T£..0U.T.eof_TAME_>j_JLOU _  _ 

iiAiii.  iiem,  ucm.  izm.  uin 


9  •  iiAiiii  nem.  ucm 

WRITE  OUTFUT  TAME  5j_l.0l51 _ 

9  .(bAlk.n.K-E.A) 


WRITE  OUTPUT  TAPE  5.  1015 
.9  ,  IJllK..I.UOllJ(jIJu«teli|C.U„ 


._.RUK..NUMOER_tONSI4TE»lCV  CHECK  _ 
IF  (NftUN  •  NRIMTI  200.55.200' 


C  CONTROL  INTEGER  CHECK 

.C _ _ _ 

55  IF  ilCTL  •  tCANTl  200.40.200 

c  ““option '^e’ction 


70  GO  TO  30 

C  CARb*A5SlCNHE‘Nfs~ 


701  NRUNI  •  Jll 

ICANT  ■  INTm 


IPN 

-TH 


INTI2I 

•.0ECm_ 


0ECt2> 
..OFSR.  ■.0ECI3I„ 
VOX  •  0EC<4) 
VMFO_  >_OEC<5l 

sros  •  0ECI4I 
TEETH  •  UECITI. 


.c_.. 

c 

c 


74  WRITE  OUTPUT  TAPE  5*1017 

WRITS  OUTPUT  TAPE_5#__1A11 _ 

9  •  IM.PN.OFSR.VOXfVOFO 


RUN  NUKOER  ClINSUTENCV  CHECK 
IF  INRUN  -  NRUNTr2d6*Ti»‘200  ■ 
_  CONTROL  INTEGER  CHECK _ 


_71,|F  IlCTL.-  ICAN.TI_200L75t.200 _  _ 

75  VOM  •  62366.0*iTH  ♦  27>.15)/PN 

CSTOC  •  6.0250e25/lVOM*ll.O  ♦  0P5R/t2.0»YDXm 

_ EXPERJNINTAL  0ATA_P_R0CE5SIN6 _ _ 


_ GO  .TC..I80,90».»  _LCJL. _  _ 

80  GO  TO  182,95.85).  lEXPO 
82  CALL  EXPOlINKUN.lCTi.OFSR.VOXt.VAFO  .VBM. 
1  NHUNI  ) 

J?fl-T5_ft4_ 


IFH,.0,._Mj_Ka5LB_L 


85  CALL  SXP05tNRUN.lCTL.OFSR.VOX.VOFO  .VBN.  IPN,  0.  Ut  K6bi»  tItfH 

_  .  1  NRUNT  I _ _ _ _ 

84  IF  IKOOE  -5)  81.100.  81 

90  GO  TC  50  _  _  _ _ 


C  CARD  ASSIGNMENTS 


..901  NUMCCN  ■  1C  -  5000 

IF  I  NUMCON  I  V04,  V02. 

_ S02_NRUNT  ■  Jll  _  _ 

ICANT  «  INTIl) 

DO  905  N  »  1«  IPN 


905  OIN)  •  OECIN) 

.  .  GO. TO  50. _ _  _ 

904  00  905  N  •  1.  IPN 
_909  UIN)  ■  OECIN) 


RUN  NUMBER  CONSISTENCY  CHECK, 


IF  INRUN.-  NHUNT)  2O0x9U2OO  _ 

CONTROL  intscer_j;hec.k _ 


91  IF  IlCTL  -  ICANT)  200,100.200 
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gWIM}  •  |WK  *  HXf  •  VRM  /(I'.Q  -  fNl|Mp 

EXROl  . 

..  153 

C 

EXROl 

154 

r  CAlCIHATinN  Of  EOUIVALENCE  RATIO 

EXROl 

155 

c 

EXROl 

156 

.EXROl  . 

-  157 

IRK  alPK*! 

EXROl 

158 

.  EXROl 

159 

C  .  CALCULAttON  OF  CflUC.  AMA  . . .  _ 

1 

.  EXROl  160 

C 

EKRDl 

161 

AR  •  0.0 

EXROl 

162 

00  40  K«l.  IKK 

EXROl 

U3 

OILS  •  S0RTPUH(K4ll  »  WIKIl442_4  llOUg  -  0JK41M/2.014A2I _ 

EXROl.. 

-  164 

40  AR  •  AR  4  lOIK)  4  OlK4l))#OELS 

EXROl 

165 

AINI  «  1.S70T964AR 

EXROl 

166 

C 

EXROl 

167 

C  CALCULATION  OF  FLANE  SREEO 

EXROl 

168 

.C 

EXROl 

169 

UINI  »  ONINI/AIN) 

EXROl 

170 

100  CONTINUE 

EXROl 

171 

IFtL4|RN  -56)  101.101.102 

EXROl 

172 

102  WRITE  OUTRUT  TARE  3.  1018 

EXROl 

173 

L  -  0 

EXROl 

174 

C 

EXROl 

175 

C  RRINTOUT  OF  FLAME  SREED 

EXROl 

176 

C 

EXROl 

177 

101  WRITE  OUTRUT  TARE  3.  1016 

EXROl 

178 

00  110  N«l.iRN 

EXROl 

179 

110  WHITE  OUIRUT  TARE  3.  lOlT 

EXR01_ 

180 

9  .N.  OFOINI.  OXOIN).  FNilN).  ONINI 

1  XINIi  UINIt  OINI  " 

EXROl 

181 

NHUNT  •  IRUNT 

EXROl 

..182 

120  RETURN 

EXROl 

183 

C 

EXROl 

184 

C  DIAGNOSTIC  -  UNACCERTADIE  CARO  ORDER 

EXROl 

185 

C 

EXROl 

ll>6 

200  NEK  -  2 

EXROl 

187 

WRITE  OUTRUT  TARE  3.  3000.  1C.  Jll 

EXROl 

188 

KOOE  •  1 

EXROl 

189 

GO  TO  120 

EXROl 

190 

C 

EXROl 

191 

C  FORMAT  STATEMENTS 

EXROl 

192 

C 

EXROl 

193 

EXROl_ 

...194 

9  IIHO.  lOX.  22NfLAM€  FRONT  OINENSIUNS  /  1H0» 

4H  SET ■  ) 

EXROl 

195 

10121  FORMAT 

.EXROl 

.  196 

9  IIHO,  14,  2X.  19HS1ATI0N  IMEASUREO) 

.4X.  1019  ) 

EXROl 

197 

1013  FORMAT 

EXROl 

198 

9  IlH  .  6X.  lYMOlAMETER  IMEASUREO).  AX.  lOF 

9.3  ) 

EXROl 

199 

-EXROl  . 

200 

9  (IH  »  6X.  I4HHE1GNT  ICM.)  .  9X.  lOR 

9.5  ) 

EXROl 

201 

EXROl 

202 

9  IlH  *  6X.  14H01A)I^TER  ICM.)  .  9X.  lOF 

9.5  )-  • 

EXROl 

203 

IXROl_ 

204 

9  IlN  ,  6X.  12HSCALE  FACTOR  3X.  F  10.5  ) 

• 

EXROl 

2U5 

1016  FORMAT  ...  ,  .  - 

...  EXROl 

.  206 

9  IIHU.  4H  SET,.  7X.  9MFUEL  FLOW  »  4X»  12HOXIOANT  FLOW  .  3X. 

EXROl 

207 

1  13HM0LE  FRACTION  5X..  )1HV.OLUM1.FLOW  ■  7A, 

9HC0«MLA&CA-a _ _  - 

EXROl 

208 

2  5X.  IIHFLAME  SREEO  .  9X.  1 INEQUIVALENCE  / 

EXROl 

209 

-1  IN  .  liX.  9HICC./SEC).  6X<  9HICC./9ECr<  AX... 

AMiMUlBttn*  . 

EXROl. 

.  210 

«  •«  .  m<CC./SfClf  tX,  FNISO.  CN.I  1  4X,  9HICN./SCCI  ■ 

CXROl 

211 

9  9X  .  5HRAf|0  )..  .  .  .  . 

EXROl 

212 

49 


EX^O)  SUMOUTINi  NODI  SENT  6,  1962  HRC  OAVTON 


NONENCLAfUNE 


NAUNT 

„Q  AND  Q(N)_ 
QFSR 


-  CUMULATIVE  CONE  AREA  •Z/NI 

_»  TOTAL  C<»*^AAEA  ISO.  CM^) _ 

•  CONE  OlAMETEft  AT  ABOVE  ELEVATION  ICM.) 

^  JlEASilREO  OIANE  I  £g_A.T_  AI0V6^I.E  VAT  ION  _ 

-  INCREMENT  SLANT  HEIGHT  ICM.) 

«  REFERENCE  LENGTH  NEA5URE0  IN_VNITS_Og^  OAtIL 

•  REEERENCE  LENGTH  MEASURED  IN  CN« 

-  MOLE  fRACTI 0N_0e_J  NH 1 0 IT  DR  IN  Ml XTURJ _ 

•  ELEVATION  ABOVE  BURNER  TON  (CN.) 

_=_llEAjtOR€0  NEAX  .HEJGHT  .  _ 

>  ONTION  CONTROL  INTEGER  CHECK 

-  MgASURED  ELE IQN  ABOVE  BURNER  TON _ 

-  TOTAL  NUMBER  OF  CONE  MEASUREMENTS 

>  NUMBER  OF  EXPERIMENTS  _ 

-  NUMBER  OF  OIAHETEK  MEASUREMENTS  NER  RUN 

-  COUNTING  INTEGERS  „  _ 

•  SIGNAL  INTEGER  FOR  ACCEPTABLE  DATA 

»  RUN  NUMB ER  OF  DATA  CROUP _ 

•  RUN  NimBER  CHECK 

_-_EOUiVALENCE  RATIO— _ 

-  ATOMS  OF  OXYGEN  TO  COMPLETELY  0XI0I2E  ONE 

_ M.QLftUiE  OF.  .F«6L. _ 

-  FUEL  FLOU  ICC. /SEC  AT  VBFO  MOLAL  VOLUME  I 


TEMPERATURE  AND  PRESSURE  I 
_=-fMtOANT  FLO>H.CC./SEC.Ar,.0_PEfi.  C«  76Q  WUJI 

-  SCALE  FACTOR  ICM. /UNITS  OF  DAIKU 

-  ACTUAL  DISTANCE . BETVEEN  TEETH  TINS  SEEN  ON 

SCHLtEREN  PHUTOCRAPH  ICM. I 
^\PLAM£_SPEED  ICM./SEO _ 

-  VOLUME  PER  MOLE  OF  FUEL  ICC./CRAH-MOLEI 
-^VOLUME.  PER-MOLE  OF  MIXTURE...  ICC./GaAIbJlOCEJ- 

-  FUEL  FLOW  IGRAM-MOLES/SEC.) 

_=_0/l  tOANT. .  FLORL.J  CRAMrMOLES/.SiC,. » _ 

*  MOLE  FRACTION  OXYGEN  IN  OXlOENT 


C  VBFO 

C _ V8M  . 

C  WF 

C‘_^ _ VX.  _ 

C  YOX 


SUBROUTINE  EXPOSINHUN* 1CTL*OFSR»YOX» VBFO»VBN» IPN»OtU*ROD£| TEETH* 

_ L-  .MRUN1__I _ _ _ 

COMMON  INT*  DEC*  IC»  Jll*  J2*  JS,  NIN*  NEX 

_ OIMENS.IOW-.IMHUI>  OECU.GJ _ _ 

OIMENStON  AILOOl,  OTIOOI*  DAIIOOI*  HilOO),  INIIOOI*  FMIIlSt 


WRITE  Qurpur  TAPE  )*  1010 

KOOC  •  S  “  *  ~ 

00  100  N«1*1PN _ 

*  10  CALL  INPUT 

_ NEX.  .HEX  -  - 

|F(  NEX  •  41  200*  It*  200 

1 1  NUMCON  IC.*  4000  . 

1F(  NUMCON  1  IS*  12*  1} 

""C  0~CARO  ASSIGNMENTS 

_ _ _ _ 

12  NRUNT  •  Jll 

..  .  ICANT  •..INTIII.  _ _ 

IRK  -  INTIEI 

_  UFOIN^  ■  OECtll 

OxblNl  •  OECUl 

_ FMlINl..-  OECI.?! _ .  - 

HPKA  -  OECIAI 

6LCM  .  •  0ECI5I _ 

ELA  -  OECIAI 

_ CO  T0_10 _ 

U  IFI  NRUNT  -  Jill  200*  141*  200 

_ 14ULIMIT^  Jl05LllMM£aH _ 

IGO  •  LIMIT  >  V 

_ .IFI  LIMIT  ~  IMK  1  142.  142*  14J 

140  LIMIT  >  IRK 


DO  143  K  •  IGO*  LIMIT 

_  M  •  ♦  I  _ 

IHIKl  •  INTIMI 

_ 143  OAIKI  ■  OECIMl _ 

IFI  LINff  -  IRK  1  10*  144* 
144  IMUM  *NRUM  ♦  N  ^1 
IFITEEINI14»I5«14 

_ 19  TEETH  •0.2 _ 

•  14  SF  •  ELCM  /ELA 

_  IPX  •  IRK  ♦  I _ • 

*C  "  ’ 


.  .  1FUPK.-I.0ftll9*l9,30 _ _ _ 

to  00  20  K«  I  •  IRK 

_  HIKI  •  TEETH  •FLOATFUHIRI-IJ _ 

20  OIKI  •  0A|R1*SF 

CALCULATION  OF  PEAR  HEIGHT  IN^. 

iillPKl  AAiPKAVrE 

01  IPX  I  -^O^O _ 

JPKl  •  I 

PRINTOUT  OF  FLAME  FRONT  OlMENSToNS 


C  TEST  FOR  MORE  THMI  TEN  MEASUREMENTS 

IFIIFR-WI  24*24*2'i  ^ 

_ 2  U0O..2  }.Jf  _ 

IF(L-4012M*2I0*2>1 

.291  MRITE  OUTPUT  TAPE.  _ 


51 


SUBftduTiNE  EXBD) 


ICONTINUIOI 


CHMBiNeNlAL  FLANC  sBC^o  datIT  AlAudfldN 


„ PRINTOUT  POt  CROUBS  Of  T^N  WeASURgMEWTI 


.8.10  UHITC  0UTRU1  TAP£_J,_JLftjL2li,.  .  ,  . 

H  IRUN.  IlHIK).  K-  JPKlt  JBK  I 

WRITE  0U1RUT  TARE  3.  lOlIt  . . . 

4  (OAlKt.K*  JRKlt  JPK  ) 

_  WRITE  OUTPUT  TAPE  3.  lOiAt 


(HIKIf  K>  4PKi(  JPK  I 
WRITE  OUTPUT  TAPE  3,  lOlSi 
*t  (OIKI*  K«  JPKl*  JPK  ) 

JPKt  •4PK  ♦  1 

23  L«  I*  » 

24  IFIL»AWI2S.2S.241 


241  WRITE  OUTPUT  TAPE  3* 

_  t»0 _ _ 


1020 


C_  PRINTOUT  PjpR  LESS  THAW  TEN  NEASUREWENTS 
C  HITE  OUTPUT  T  PE  3.  1012. 


9  iKUNt  HPKA*  (IHUJ*  K*  JPKIt  IRK 

_  WRITE  OUTPUT  TAPE_^3.*  _ _ _ 

9  lOAIKlt  K>  JPKl*  IRK  I 

_ WHITE  0UTPi>T_fAP£_>*_l0l4, _  _  _  _ 

9  IHIKI*  K*  JPKl*  tPK  1 

WRITE  OUTPUT  TAPE  3.  101S> 


lOIRl*  K«  JPKl*  IPK  I 
WRITE..0U.TP.0.T.  rAPe_J,_JL0.l.5l#._  _ 

9  ELCH*  ELA 

f  X*X<f. _ _  _ _ _ 


IP  (IRUfi-  NRtmT)  30«3l|30 
30  KOOE  -  I 

IRUNT  •  NftUNT  _  _ 

GO  TO  100 


CONTROL  INTEGER  CHECK 
IP  IICTL  ‘-‘iCANTT“30'*‘35*30” 


CALCULATION  OP  OXIOANT  AND  PUEi  PLOW 


WX  •  0X0(N)/22414.0 
WP  •  OFOINl/VefO. 


c 
c 

.c. _ 

c 
e 
c 


OH(N)  •  (WP  ♦  WX)  •  VBM  yil.O  -  PNl(N)) 


CALCULATION  OP  EUUIVALEITCE  RATIO 


OIN)  •  UP*0PSR/(2*0*WX*V0X) 


CALCULATION  OP  CONE  AREA 


AR  •  0.0 _  ■ 

00  40  R  •  I*  IRK 

DELS  •  SOHIPttW(K»n  -  WU))»t2  ♦  HOIK)  -  D(K*i))/l*QlPP21_ 


40  AR  •  AR  ♦  (OIK)  ♦  0(R«1))*DELS 
AIN)  •  l..ST0796*AR  . 


_ CALCULATION  OF  PLANE  SPEED 


C 

-C_ 

C 


_  U(N)  •  QM(N)/A(.N) _ 

100  CONTINUE 

IP(L«IPN  >96)  101*I01.*102 
102  WRITE  OUTPUT  TAPE  3*  lOlB 

_L-‘i _ 


PRINTOUl  OF  FLAME  SPEED 


101  WRITE  OUTPUT  TAPC_At. 
00  lib  N-i*lPN 
J|  RUN  ■  NRUN  -1  ♦N 


1016 


110  WRitE  OUTPUT  TAPE  3*  1017* 

_ 9  IRim 

NRUNT  •  IRUNT 

120  HETUHN  _ 

C  DIAGNOSTIC  -  UNACCEPTABLE  CARP  ORDER 


200  NEK  ■  2  _ 

WRITE  OUTPUT  TAPE  3* 

_  KOOE  •  l_  _ 

GO  TO  120 


3000*  1C*  Jll 


PORNAT  STATEMENTS 


1010  PORNAT 

9  (IHO*  lOX* 
1012  FORMAT 


22HPLATIE  FRONT  DIMENSIONS  /IHOi  4H  RUM  ) 


9  IIHO.  14.  2X.  2RHPEAK  HEIGHT  IMEASURED  UNITS)  FB,2  / 
1  IN  »  6X*  I9HSTATI(M  (NEASUREO)*  4X*  1019) 


10121  FORMAT 

9  (IHO*  14* 

1013  FORMAT _ 

9  (IH  * 

1014  FORMAT 
9  (IH  • 

1019  FORMAT 
9  (IH  • 
10151  FURMAT 


19M$TATlbN  (MEASURED) •  4X(  1019) 
“l9H0IAMiTCR"(MEAS0RE0l*  4X,"l0F  9.3  )~ 


6X* 


14HHEIGHT  (CM.) 
~I4M0VaMETER‘|CM. I  ■ 


9X(  lOF  9.9  ) 


9X(  lOF  9.9  ) 


9  (iH  *  6X*  26HACTUAL  LENGTH  OF  REFERENCE  3X*F10.B*3I*  3HCM* 

_  . 1 . JU«^28!!ifilU«Ui£JLJL^^  . .  . 

1016  FURMAI 

9  (IHO*  4H  RUN*  TK*  VHFUEL  FLOW  *  4X*  12H0X10ANT  FLOII.X  IKl 


EXPOS 

lOT 

EXP03 

lOB 

EXPD3 

109 

EXPOS'" 

“uo 

EXPOS 

111 

EXPOS 

112 

EXPOS 

US 

EXPOS 

114 

EXPOS 

US 

SXP03 

"UG 

EXPOS 

117 

EXPOS 

iia 

EXPOS 

119 

EXPOS 

120 

EXPOS 

121 

EXPOS 

122 

EXPOS 

123 

EXPUi 

124 

EXPOS 

12S 

EXPOS 

126 

EXPOS 

127 

EXPOS 

12H 

EXPOS 

129 

EXPOS 

ISO 

EXPOS 

131 

EXPOS 

•132 

EXPOS 

ISS 

EXPOS 

134 

EXPOS 

139 

EXPOS 

136 

EXPOS 

137 

EXPOS 

13B 

EXPOS 

139 

EXPOS 

140 

EXPOS 

141 

EXPOS 

142 

EXPOS 

143 

EXPOS 

144 

EXPOS 

149 

EXPOS 

146 

EXPOS 

147 

EXPOS 

146 

EXPOS 

149 

EXPOS 

150 

_ IXPOS._ 

191 

EXPOS 

152 

..  .  EXPOS 

ISS 

EXPOS 

194 

EXPOS 

159 

EXPOS 

196 

_  EXPOS 

.157 

EXPOS 

ISB 

EXPOS 

159 

EXPOS 

160 

EXPOS 

161 

EXPOS 

162 

EXPOS" 

163 

EXPOS 

164 

EXPOS 

165 

EXPOS 

166 

EXPOS 

167 

EXPOS 

168 

EXPOS 

"169 

EXPOS 

IFO 

EXPOS 

171 

EXPOS 

172 

EXPOS 

179 

EXPOS 

174 

EXPOS  179 

EXPOS 

176 

EXPOS 

177 

EXPOS 

17a 

EXPOS 

179 

EXPOS 

IBO 

EXPOS 

"iBl 

EXPOS 

182 

EXPOS 

183 

EXPOS 

1B4 

EXPOS 

IBS 

EXPOS 

186 

EXPOS  197 

EXPOS 

168 

EXPOS 

189 

EXPOS 

190 

EXPOS 

191 

EXPOS 

192 

EXPOS 

"19S 

EXPOS 

•  194 

EXPOS 

199 

EXPOS 

196 

EXPOS 

■  197 

EXPOS 

196 

EXPOS  199 

EXPOS 

200 

EXPOS 

201 

EXPOS 

202 

EXPOS 

2US 

EXPOS 

204 

EXPOS 

■20> 

EXPOS 

206 

EXPOS 

207 

EXPOS 

208 

EXPOS 

209 

.  EXPOS 

«210 

EXPOS 

2U 

EXPOS 

212 

52 


SUftROUTlNC  fXPOS 


ICONTINUfO) 


I  IIHUOII  FMCflON  •  SR«  IlHVOLUNE .  FLOH  ,  7I»  . 

I  SXi  IIHELAMC  SMCO  •  SXi  UHEOUIVALCNCE  / 

_ JL  IM  ,  UX,  Wiy./iEO^  4X.  9HICC./SEC>  .  4X,  jHjWlfiTgU _ 

A  EX  t  9Hl7C«/iECIt  WUO.  CN.»  »  AX*  «HlCiU/SEC)  » 

_ 5;  9X  ♦....SHXA.aO..  .1 _  _  ...  _ 

1017  FORMAT 

9  IlH  •  lA,  AX.FtO.9«AX.FIO«i.AX.FlO.5«AX.FI0.»i6XtF10j»iAXrFJfliA 
1  AX(FLO.S) 

.JJUM.  £0AMA.T _ _ _ 

9  UHll 

_10M  FORMAT _ _  _ _ _ _ 

9  UHU  lOX*  ElMELAME  FRONT  DIMENSIONS  /IHO  •  AH  MM  I 
SOOO  FORMAK  .iHO*.  2SX«.  AlNOAlA.CARO  OROCR.  INCORUCT.^  CASPJMHiilL. 

1  tA»  SX«  I2HRUN  NUNOER  If  I 


EXFO)  21) 
fXFO)  21A 

EXFO) _ 21S 

“■EXFD)  2U 
EXFO)  217 
~EXFO>  21i 
EXFO)  219 
~'EXF0)  220 

EXF03 _ 221 

'~EXFO)  222 
EXFO)  22) 
EXFO)  22A 
EXFO)  .  22) 
EXFO)  226 
EXFO) _ 227 


LLT  ft. JT.M mu 414.1,1  l.rjlL«L 


CHAXM  SUIAQUTINC  MAXN  PON  ROUTINE  20|lt  FEB  2B  1942  'hRC~  OAVTON  ' 
C  NOMENCLATURE  ~ 


c 

A(R) 

c  . 

ARG 

c 

aij.x) 

C 

c 

IPN 

J£, 

c 

K9 

,c 

..  RODE 

c 

MAXQN 

c 

c 

NCF 

c 

0(JU 

c 

OM 

c 

OINF 

c 

c  .  _ 

OL 

c 

OMX 

j; _ 

STDS 

C 

c_.  .. 

, .  SUM 

c 

UIN) 

.c _ 

_ UOIN) 

-  COEFFICIENTS  OF  FITTED  COUATION 

^*JMIA.Vf«ENT  OF  QUADRATIC  SQUARE  ROQT _ 

-  SUN  OF  SQUARES  AND  CROSS  PRODUCTS  OF  THE ~ 

_ UlNfKI  MATRIX 

-  NUMBER  OF  DATA  POINTS 


N9  -  SIGNAL  FLAG  FOR  SEQUENCE  OF  CALCULATIONS 

.  MODE . . ::_SICNAL  INTEGER  FOR  ACCEPTABLE  DATA  _ 

MAXQN  -  SNITCH  TO  BYPASS  MAXM  SUBROUTINE 

_  ^  (ANY  NON-ZERO  VALUE  BYPASSES  SUBROUTINE) 

NCF  -  NUMBER  OF  COEFFICIENTS 

^0  (Nl  -  EOUIVALEMC  E_RA  T 1 0 _ 

CIM  -  LARGEST  EQUIVALENCE  RATIO  INTET 

..  .OI.NF _ -  EOUIVALFNCE  RATIO  AT  INFLECTION  POlNl.  D<LCU81C 

CURVE 

.  OL _ ::_SNALLEST  EQUIVALENCE  RATIO  IN  SET _ 

OMX  -  EQUIVALENCE  RATIO  AT  NAXiMUM  FLAME  VELOCITY 

_ -  MAXIMUM  ALLQNABLE  STANpARO  DEVIATION  OF  U  VS 

0  CURVE  FOR  ADDITION  TO  TAPE 

.  SOM  .  _ _ =_SUH  OF  THE  SQUARES  OF  THE  OE)tIAULQ*lS _ 

UIN)  -  FLAME  SPEED  (CM. /SEC) 

_OOIN) _ ^_r_DEV)ATION  BETHEEN  MEASURED  AND  PREDICTED _ 

FLAME  SPEED  (CM. /SEC) 

PRCDlCT‘E0“'F?AHrTplEb~*^  SDRED  Kp 

..  UHX _ r-MAXlMUM  .FLAME  SPEED  (CM. /SEC.)  _ _ 

UP(N)  -  PREDICTED  FLAME  SPEED  (CM. /SEC) 

..USTO _ =_STAN0AkD  DEVIATION  OF  FLAME  SPE|Ql  |CN./:S60 _ 

USTOC  -  FLAME  SPEED  AT  STOICHIHETRIC  CONDITIONS 

_ U.I1.  /.SiCi . . . 

X(N»R)  -  INDEPENDENT  VARIABLES 

SUBKOUTJNE.  NAXMI..|PN,..0*  U,  OMXi  UNX*  USTOC.  RODE.  STDS.  _ _ 

1  lEXPO.  MAXON.  NRUN  ) 

CQNHON  INT._O.ECj_JICa.  illlt.  J2|  JB»  NIN.  NER.  _  _  _ 

DIMENSION  INTIIO),  DECIIO) 


MAXM  1 

MAXM  2 

MAXM  3 

_MAXM _ 4 

MAXM  '  "•» 

MAXM  6 

MAXM  7 

MAXM  8 

MAXM  V 

MAXM  _ 10 

MAXM  11 

MAXM  12 

MAXM  13 

MAXM  14 

MAXM  IS 

MAXM _ 16 

MAXM  17 

MAXM  10 

MAXM  19 

MAXM  20 

MAXM  21 

_MAXM _  22 

MAXM  23 

.MAXM.  24 

MAXM  2b 

MAXM  26 

MAXM  27 


_ KOOE  -  0.  ....  _ _ _ _ 

IF(MAXON)  SOO.  20.  5D0 
-..20.N«fT£  OU.TPUT  .tAP.E 

IF  (IPN  •  2)  500*500,2S 

_ 2  JLOi-  A  JftUJ _ 

OH  •  OL 

00  30  N*1.1PN  _ _ _ 

OL  •  HINIFfOL.OlHI) 

UH  •  HAXIFIOH.OINJJ _ 


*IN.*.l)_f.  UO  .  .. 
X(N,2)  •  0(N) 

_ F < N.J)_»  0LNJf2_ 

X(N.4)  •  01N)«*3 

-  30  MIN. 5.)_«_0,(JI ) _ 

C 

.C. _ _  TEST  eOlLHORfeJUJi 

C 

_ IF  ilPN  -^)  ISO. 

C 


_ SOLUtlON  OF  /OyR_lAumQNS_IH_Fjlu.LU||RN0WlS^ 


.•.4L,WR| TE  OUTP.UT„.TAPt >t-Ji)lQ. _ 

WRITE  OUTPUT  TAPE  3.  1011. 

_ .9_. .  I(XIK.LI.L»1.3  I.K»1.4) _ 

GO  TO  200 

ATTEMPT  TO  FIT  FOUR  OR  MORE  DATA  POINTS  TO  A  FOUR  CONST 
_  ..  _E«UAUUN _ 

„9Q_00_>.0_J«U4 _ 

00  60  K*l«5 

CUj HI  ■  ftJ _ 

00  60  N«I.1PN 

'  generation  oTI^uir  regress ion'Touamons 

60”6( j.R)  TTu.'x)  ♦aiN.J)«X(N|K) 

SOLUTION  OF  four  REGRESSION  EQUATIONS  IN  FOUR  UNRHONNS 
'CAU  CROUf  'iVTcTa) 

_ IFI  .SENSE.  LIGHL_ai_4i._jyi _ 


61  WRITE  OUTPUT  TAPE  3._IU10._.  _ 

WRITE  OUTPUT  TAPE  3.  1011. 

9  IIG<RtL)4l«l«9..).».R£Uf^J _ 

GO  TO  200 

WRAtTuNS  UN'I^UR  or  more  OATT  POINrS 


_MAXH _ 64 

HAXH  65 

_MAXM,._.  66 
HAXH  67 

_MAXM _ 6B 

MAXM  69 

_MAXM  _ _ 70 

MAXM  71 

_MAXH _ 72 

MAXM  73 

_MAXM _ 74_ 

MAXM  75 

_MAXM  _  76 

MAXM  77 

_MAXM  _  78 

MAXM  79 

.MAXM _ ^BO 

MAXM  ai" 

_MAXM _ 82 

MAXM  as 

MAXM  64 

MAXM  BS 

.MAXM _ ^86 

MAXM  87 

..MAXM...  68 

MAXM  B9 


_■ _ MAXM 

-  90 

NAXM 

91 

_ ^MAXM 

—92 

MAXH 

93 

RAXM  .. 

.  94 

MAXM 

95 

MAXM 

96 

MAXH 
-  HARM _ 

97 

—9a 

MAXM 

99 

MAXM 

100 

MAXM 

101 

.  .  MAXM 

102 

MAXM 

1U3 

_ MAXM  .. 

.104 

MAXM 

lOS 

MAXM 

106 

54 


su«iiouiiNe  'naxm 


NCr  •  4  .  _ 

K4  •  I 

.1  f  OlViOt  CMeCK  Ti« 
UMX  ■  0.0 
OMX. 


HAXM  107 
MAXM  lOt 

_MAXM _ 109 

110 


QlNf  •  -AOl/iS.O^AIAll 

IF  OlVI.OE  _ 

K9  •  2 

_ARf.  3L^*AAll#A14iyAia»»2_ 


IF  OlVlOC  CHECK  310*  tO 


TEST  FOK  SOUAAE  KOQT  OF  NE6AT1VE  NUHOEA 


IF  URU)  3t0«)10ta5 

J _ 


0 1  Wf •  M  .0  ♦_tl  GMF  <  SORTfURtfl^  Aiyj_L^ 


c 

C  TEST  Jto  it£  IF,  CtHJIVALtilCC  RATIO  AT  HAKlHUli  FLAWI  jPBBP  IS  MITHIN 

C  OAIA  RANGE 

IF  IGNX  -  0H»  90,‘90*i00' 

40  IF  IQL  -  ONKt  9S.9S.300 _ -  _  .  . 

9S  K9  «  4 


C 

>.c._ 

c 


QHIN  »  PIMF  *11.0  .  SICtlFIR0tTFtARGI>AI3>n_ 


.  -  .r.E.S.r  ..FJ).lJHlNillUIU»01Hl-JM-OAIA-M^ _ 


_  1  F.rON  l.»l_=JJHULOA*.Ufl*.Ui>_ 
100  IFiOL  •  OHIN)  300«110*110 

_UA-R.9_-_J! _ 


C 

C 

C 

C 

Jt- 


ENO  or  OFERATION  ON  CUOIC  FIT 


ATfEHFr  fO  FIT  THREE  DATA  FOINTS  TO  A  THREE  CONSTANT  EQUATION 


NAXM 

_NAXK_ 

NAXM 

_HAXH_ 

^'maxh 

NAXN_ 
NAXH 
_HAXM  , 

"naxm  ' 

MAXN 

MAXM 

_HAXM 

"'haxh” 

_HAXH 
MAXH  ' 
NAXM 

'naxm  “ 

_NAXM_ 
NAXM 
_NAXN  „ 
MAXM 
_NAXM_ 
MAXH 
.NAXM_ 
MAXH 
-HAXM... 

NAXM 

.MAXH 

MAXH 

_MAXM_ 

NAXM 

.NAXM 

MAXH 

MAXM 

MAXM 

_MAXH_ 


.111 

112 

_117 

114 

1ft 

116 

117 

lie 

119 

I2U 

121 

122 

123 

124 
.125 

126 

-127. 

126 

129 

130 

131 

132 

-133. 

134 

13$ 

136 

137 
136 

.139 

14U 

141 

142 

143 

144 
.14$ 


MAXM 

147 

CALL  CROUr  I3*X»AI 

MAXH 

146 

HAXM 

149 

C 

NAXM 

ISO 

MAXH 

-iSl 

C 

MAXH 

1S2 

MAXH 

IS3 

HKirs  ourFUf  tare  3#  1013# 

MAXH 

164 

9  ((XfK«LbL*l«.4).d(9Xal) _ 

MAXM 

ISS 

GU  TO  soil 

HAXM 

1S6 

.  MAXH  . 

.157 

C  ATTEMMr  TO  FIT  FOUR  OR  MORE  DATA  FOINTS  TO  A 

TNREE  CONSTANT 

MAXM 

1^6 

c  eOUATION  . . 

MAXM 

159 

C  _  ‘  - 

..MAXM 

160 

200  00  210  N*t*iFN 

MAXM 

161 

210  X(Ni4}  •  U(N) 

MAXM 

162 

00  220  J*l»3 

MAXM 

163 

OU  220  K«1.4 

MAXM 

164 

GU»K1  •  0.0 

HAXM 

16S 

00  220  N«1.IFN 

MAXH 

166 

C 

HAXM 

167 

C  GENERATION  OF  THREE  REGRESSION  EQUATIONS 

HAXM 

166 

C 

MAXM 

169 

220  GIJ.K)  •  GU.KI  ♦  XIN.JIMIN.X) 

MAXH 

170 

C 

HAXM 

171 

C  .  SOLUTION  OF  THREE  EQUATIONS  IN  THREE  UNKNOWNS 

MAXM 

172 

C 

MAXM 

173 

CALL  CR.OUT  (3.G.A> 

MAXH 

174 

C 

MAXM 

17S 

C _  TEiX  FOR  UNSUCCESSFUL  SOLUTION 

MAXM 

176 

C 

MAXM 

177 

IFISENSE  LIGHT  3  1  221,  230 

MAXM 

176 

C 

MAXM' 

179 

C  FRINTOUT  OF  UNSUCCESSFUL  FIT 

MAXM 

160 

C 

MAXM 

161 

.  MAXM 

162 

WRITE  OUTPUT  TAPE  3»  1013» 

NAXM 

163 

9  IIG(K.L>.L»1.4t«R«tf3) 

* 

MAXM' 

164 

GO  TO  SOU 

MAXH 

las 

c 

MAXM 

166 

C  OrEMTIONS  ON  IHII»  0*  NONE  O.TA  .OINIS 

MAXM 

167 

HAXM 

166 

230  NCF  •  3 

• 

MAXM 

169 

MAXM 

190 

AI4I  •  0.0 

MAXM 

191 

HAXM 

192 

OHX  •  0.0 

NAXM 

193 

MAXM 

194 

240  K9  •  6 

MAXM 

19S 

MAXM 

196 

C  OEfCIININ.rlON  OF  EOUIV.UNCE  MTIO  tT  NMINUN 

FLANE  SPEED 

HAXM 

19/ 

MAXM 

^196 

ONX  •  -AI2l/(2.0»A(3n 

MAXM 

199 

NAXM 

2UU 

C  TEST  TO  SEE  IF  EQUIVALENCE  RATIO  AT  NAXINUN  PLANE  SPEED  IS  NITNIN 

MAXM 

201 

NAXH 

202 

c 

NAXM 

203 

MAXM 

.204. 

230  IF  tOL  *  ONXI  2A0t2A0t300 

MAXM 

20S 

c  - - 

MAXM 

206 

C  END  OF  OFfiRATION  ON  PARA30L1C  FIT 

MAXM 

207 

MAXM 

206 

C  OFtlltriON  ON  OAT.  AFTfA  SUCCESSFUL  FIT 

MAXM 

209 

HAXM 

210 

260  K9  «  T 

MAXM 

211 

C  .  ...  — - - 

—  .  -  ....  , 

55 


XCOatUMlDI 


MAXIMUM  fLANE  SMEED  OETERNINATION 


CALCULA.TtON  OF  MAX.I.NUM  FLAME  SFEED  _ _ 

iX*IAI3»  ♦  OMX»AIAlli 


e 

c _ _  CAtCML.AltQ.N.JJf-fXAMl_$J«.eo.AI  AN  EOUIVALEHCi  MAIIO  OF  ONE 

JIQ  USTOC  ■All)  »  AI2I  ♦  At>i  ♦  AIAI  _  _ 

SUN  »  0.0 
00  a20  N«1.1FM 


UPIN)  •  Ad)  ♦  0(N)»1AI2)  ♦  0(N)«IAI3)  ♦  0<N)*A{4))) 


COMPARISON  OF  PREOICTEO  ¥S«  MEASURED  EOUIVALENCE  RATIOS 
UOIN)  ■  U(N)  -  'upTnI  . .  “  ‘  ~ 


CALCULAdON  OF  THE  PERCENT  ERROR  Of  THE  OEVUTSON 


UOP(N)  •  lUU.O*UOtN)/UP(N) 


SUMMATION  OF  THE  SQUARE  OF  THE  DEVIATIONS 


.  MAXM  21) 
MAXM  2IA 

MAXM _ 215 

~’nAXH  216 
_MAXM  „217 
MAXM  2U 

_MAXM _ 21V 

MAXM  220 

MAXM _ 221 

MAXM  222 
_MAXM_  22i 
MAXM  224 

MAXM _ 225 

MAXM  226 

_MAXM _ 227 

“mAXM  ‘22« 

_HAXM _ 22V 

MAXM  2)U 
_MAXM  _  231 
MAXM  "232 


c 

)20  SUN  -  SUM  *  U0IN)**2 

MAXM 

234 

If  (IPN  -  NCF)  330.330.140 

'i)o  usro  «  0.0 

MAXM 

236 

CO  TO  3S0 

c 

MAXM 

238 

t  CALCULATION  OF  STANDARD  DEVIATION  OF  FLAME  SPEED 

MAXM 

240 

.  J.40.UST0  -  SORTFISUH/FLOATFIIPN-NCF)) 

MAXM 

350  WRITE  OUTPUT  TAPE  3|  1014 

MAXM 

242 

..WRITE  OUTPUT  .TAPE  3t  1015 

WRITE  OUTPUT  TAPE  3,  1016* 

MAXM' 

244 

4  (A(K).K>i.4).USTD.0MK.UMX 

IFdEXPO  -  1)  351»352*)S) 

MAXM 

351  WRITE  OUTPUT  TAPE  3.  10171 

MAXM 

24r 

00  354  N«l»  IPN 

MAXM 

248 

IHUN  •  NRUN  «  N  -1 

354  WRITE  OUTPUT  TAPE  3*  1010 

MAXM 

250 

«  .  IRUN.  OINl.  UIN).  UPINI.  UDfN).  UOPIN) 

CO  TO  353 

NAXK 

252 

352  WRITE  OUTPUT  TAPE  3.  1017 

00  3521  N  •  1»  IPN 

MAXM 

254 

3521  WRITE  OUTPUT  TAPE  3.  lOlA 

MAXM 

255 

4  •  N  •  01N)»  UINIt  UPIN)t  UDtN)»  UDPIN) 

MAXM 

256 

353  IF(I14-4I  360.360.370 

MAXM 

258 

:  ERROR  PRINT  •  CUdlC.  P1T..0F  CURVE  J.>N&ATI.SFACXOIU _ 

„  .  MAXM 

259 

HAXK 

260 

360  WRITE  OUTPUT  TAPE  3*  .1014 

MAXM 

261 

00  TO  (361(362»363(364),  K4 

MAXM 

262 

:  OIAONOSTIC  •  PARABOLIC  FIT  INDICATED  »  A3  SMALL  OR  TEPO 

MAXM 

264 

1  .  ...  , , 

MAXM 

265 

361  WHITE  OUTPUT  TAPE  3.  1020 

MAXM 

266 

CO  TO  200 

MAXM 

267 

c 

MAXM 

268 

C  ERRGK  PRINT  •  NO  MAXIMUM 

MAXM 

264 

.c  . 

MAXM 

270 

362  WRITE  OUlPUT  TAPE  3*  1021 

MAXM 

271 

00  10  200 

MAXM 

272 

c 

MAXM 

273 

C  ERROR  PRINT  -  MAXIMUM  OUTSIDE  DATA  RANGE 

MAXM 

274 

MAXM 

275 

363  WRITE  OUTPUT  TAPE  3.  1022 

CO  TO  200 

MAXM 

277 

c 

MAXM 

.278 

C  ERROR  PRINT  -  MINIMUM  POINT  IN  DATA  RANGE 

MAXM 

279 

c 

MAXM 

280 

364  WRITE  OUTPUT  TAPE  3t  1023 

MAXM 

281 

<iO  TO  200 

MAXM 

282 

370  IF  (K4  •  7)  3Tl»3«0tS00 

MAXM 

283 

C 

MAXM 

284 

C  ERROK  PRINT  -  PARABOLIC  FIT  OF  CURVE  UNSATISFACTORY 

MAXM 

285 

c 

MAXM  286 

371  WRITE  OUTPUT  TAPE  3.  1024 

MAXM 

287 

IF  (K4  -  61  372.373.500 

372  WRITE  OUlPUT  TAPE  3»  1021 

MAXM 

289 

CO  TU  500 

MAXM 

373  WRITE  OUTPUT  TAPE  3.  1022 

MAXM 

291 

GU  TO  500 

380  IF  (S1US)  38l»381f)02 

MAXM 

243 

3R1  STOS  >  1.0 

MAXM 

C 

MAXM 

295 

C  TEST  TO  SEE  IF  STANDARD  DEVIATION  MITHIM  SPECIFIED  RANGE 

MAXM 

296 

C 

MAXM 

297 

382  IP  lUSrn  •  SIDSI  400.400.383 

C 

MAXM 

249 

C  ERROR  PRINT  -  STANDARD  DEVIATION  TOO  HIGH 

C 

MAXM 

301 

3R3  WRITE  OUTPUT  TAPE  3.  1021. 

9  STOS 

MAXM 

303 

CO  TO  SOU 

MAXM 

400  KOOE  •  3 

MAXM 

305 

500  RETURN 

c 

MAXM 

307 

C  FORMAT  STATEMENTS 

c 

MAXM 

309 

1001  FORMAT 

4  IlHl) 

MAXM 

ill 

Kill  FORMAT  „  „  .  . . 

.  MAXM 

312 

9  IIH  «  lOXt  El2,5t  lOXt  EU.StlOXi  E12.5|I0X|  El>.5tlOXiEl2.! 

I 

MAXM 

313 

IUI2  FORMAT 

MAXM 

314 

V  UHU.  10X»  73HPAKA80LIC  FIT  OF  FLAME  SPEED  CURVE  FAILED  -  MATRIX 

MAIN 

315 

1  OF  IHE  COEFFICIENTS  IS  / 

.  ,.naxm  . 

316 

2  IHOr  lOX*  )2MEguiV«  RATIO  t  8X»  I6HEQUIV.  RATIO  ••2f  4Xi 

MAXM 

317 

3  IIHFLAME  SPEED  yy)  ...  .  _  .. 

MAXM 

iie 

56 


SUMOUTINf  MAXH' 


.ICQIlUMMfO)  . 


MAKtWiM  fikM  SmO  OtTCMINATION 


lOU  rOftMAT 

9  UH  •  lOX. 


E12*St  lOI.  E12.S,10X*  EU.SiLOKf  E12.9) 


lOX*  22HFLAIIE  SEEEO  CUtVE  DAT*  > 


9  (1M0» 

_ LaU  EOAHAT _ _ _ _ 

9  (IH  •  AXi  4<4X*  I^E12«4I«  5X»  0AFB.4  i  •X|F9.S>  9X|F9.4  ) 

._.iai7  FOAMAT _  _ 

9  (IMO.  IQX« 

_ 1. _ -UiU, 

2  IH  21Xi 

_ »  . 

lOlTt  FORMAT 

9  (IH0«..10X, 

I  tlX* 

i  >M _ 


SHOATA  SET  •tOXdlHEOUlVALENCE  lOXf  AHMEASUAEO 
9MFA60llirE0  IPX.  9H0EWIATI0M  UK.  TWFEACENT 


SHAATIO  12X<  1  inflame' STEED  IXtUMFLAME  S«Et~ 

9MOlV|AflW:  /) _  _  _  _ _ 


9H  _ ilMH  4lOX,llHEOUlVALENC.E_. 

9HFAEUICTE0  lOX*  9H0EVUTI0N 
SHAAflQ  i2K.ilMFLAM6 


10X|  BHMEASIMEO. _ 

ilXi  THFERCCNT  / 


3  2AX«  9HOEVlAriON  /> 

_.iftia.A0HMAj _ _ _ _ _ _ 

9  (IH  •  tOX*  lA*  10X,Fll.4»10X|f9.Atl0X|F9.4tt0XfF9.4|l0X|F9»4> 

_  U19  FO«MAT„ 


9  (lHO«  tOX»  ASHCUaiC  FIT  OF  CURVE  UNSATISFACTORY 
_102Q  FORMAT 


I 


9  (IH#«  4YX»  23HFARABOLIC  FIT  INDICATED  i 

tU2l  FORMAT  _  .  _  .  _ 

9  I1H«*  4TX«  tAHNO  MAXIMUM  FOUND  ) 

iOTT  FORMAT _ _ _ 

9  (IM**  47X*  26HNAX1NUN  OUTSIDE  DATA  RANGE  I 

-JLliiULQWAjt. 


9  <tH«»  47Xt  27HNtNimiM  MINT  IN  DATA  MANGE  ) 

._102A  FORMAT. _  _ _  _ _ 

9  IIHCI  6X*  39HFARABOi.lC  FIT  OF  CURVE  UNSATISFACTORY  >  I 

.1025  FORMAT_.. . . . .  . . 

9  MHO  lOX*  SlHSTMiOARO  DEVIATION  GREATER  THAN  F6.3  I 

_liaO_.F.ORMAI _ 

9  (IHO,  lOX*  49NCUBIC  FIT  OF  FLAME  SPEED  CURVE  FAILED  -  MATRIX  OF 


ITHE  COEFFICIENTS  IS. 

2  INO.  tOX  12HE0UIV.  RATIO  •  BXtlAHEOUlV.  RATIO  ••E  t  6Xi 

.  i  _ l4HE0Ul.V.«u-RATIQ.,««3.  VXiIlHFLAME  SPEED  /  ) _ 

1019  FORMAT 

i  I7HEUUI VALENCE  RATIO  ,  4X»  TMNAXIMUN  / 

.  .  2-lH  2HAI.U4J(*._2HA2.i..l4Xr  .2HA3.  »  14Xt..2HA4...«^Xi _ 

3  9HOEVlATtON  »  3X»  UHAI  MAX  FLAME  SPEED  »  3X« 

_UMFUME_SPEED._y.  .1 .  . . 


NAXM 

NAXM 

MAXM 

“maxm“ 

MAXM 

MAXM 

MAXM 

MAXM 

MAXM 

~MAXM 

MAXM 

MAXM 

MAXM 

NAXM 

MAXM 

“maxm" 

_MAXM 

NAXM 

MAXM 

MAXM 

_MAXM 

MAXM 

_MAXM 

MAXM 

MAXM 

NAXM 

MAXM 

MAXM 

.MAXM 

MAXM 

MAXM 

MAXM 

_MAXM. 

MAXM 
.  MAXM 
MAXM 
MAXM 
MAXM 
_MAXM 


END 


MAXM 

MAXM 

NAXM 

NAXM 

NAXM 


319 

320 

321 

322 

323 

324 
.325 

326 

327 
32B 

329 

330 

331 

332 

333 

334 

335 

336 
33/ 
334 

339 

340 
.  341 

342 

343 

344 

345 

346 

347 
34B 
349 
35U 

—351 

352 

353 

354 

359 
356 

—357 
35B 
.  359 

360 

361 

362 
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SUBROUTINE  TABS 


CTAPE 

C 

C 


SUBROUTINE  TAPE  FOR  ROUTINE  2016  MARCH  S»  1462  NRC  -  DAYTON 
NOMENCLATURE  . 

n  -FUEL  MANE  ’  ^ 

_ 0F.5.TRyc.TWR*A-tQW.»l-»UTiRJ.Jl!L.NQWf^^^ 

IN  REVISED  LIST 

_ :.JC.Q»t1JlBUT.0R  NAME.iJFROM  MRftSJ _ _ 

•  NUMBER  OF  STRUCTURAL  CONTRIBUTORS  PER  MOLECULE 
»  BLANK  lUSCB  IN  CONTRIBUTOR  NAME  LIST  I 


c 

AIU)  and" 

c  _ 

ACCLILiKX 

c 

c . 

.  ACNLIR).  . 

c 

tSl(Kil) 

c 

BLK 

c 

U(KI 

X 

ONMIK) 

c 

UU  ANO 

c  . 

FGCJL(J) 

C 

.« _ 

FCCL(K) 

TAPE 

TAPE 

TAPE 

TAPE 

tape' 

TAPE 

TAPE 

_ TAPE 

TAPE 

_ TAPE_ 

TAPE 

TAPE 

TAPE 


C._ _ 

c 

c_ 

c 

X _ 

c 

c  .  . 
c 

c 

£, _ 

C 

C _ 

c 

c _ 

c 

c 
c 
c 
c 
c 

c_ 


oNLIJ) 

J»  J.  ,Kj_L _ 

IIAIU 

iiBin _ 

ncd) 

I2UI  .  _ 
13(1) 

1A(  I) _ 

lAOO 

ICAHT 

ICTL 

IFNC  _ 

IFNG 

IFNM _ 

INOS 

J.N.£,C _ 


-  DECIMAL  DATA  STORAGE 

_z_0£p6N0ENI  VARIABLE  NAME  .LIST _  _ 

•  FIRST  AND  SECOND  HALVES  OF  DATE 

-  NUMRFR  OF  STRUCTURAL  C0NTRI80T0RS._R6fi_M0LECULE _ TAPE, 

FUR  CONTRIBUTOR  NUMBER  LFCCCNiKI  TAPE 

-  NUMBER  OF  STRliCTURAL  CONTRIBUTORS  PER  MOLECULE _ TAPE  ^ 

FOR  CO'NIRlBUrOR  NUMBER  LPGCCNU)  TAPE 

-  FIRST  AND  second  HALVES  OF  FUEL  .HRHC _ TAPE  , 

-  CONTRIBUTOR  NAME  (FROM  TAPE!  TAPE 

-  INDICIES  _  _ tape., 

-  FUEL  CLASS  NUMBER  TAPE 

-  FUEL  GROUP  NgMLEJL _  TAPE 

•  FUEL  MEMBER  NUMBER  TAPE 

-  DATA  SOURCE  NUMBER  _ TAPE., 

-  EXPERIMENTAL  CONDITIONS  NUMBER  TAPE 

-  NUMBER  OF  STRUCTURAL  CONTRIBUTORS  CONSIDERED _  .TAPE  . 

TAPE 

_TAPE_ 
TAPE 
,TAPfc 
TAPE 
TAPE 


•  NUMBER  OF  DATA  GRUUPS  ACCEPTABLE  FOR  WRITING 


-  OPTION  CONTROL  INTEGER 
«  FUEL  CLASS  NUMBER 


o  FUEL  GROUP  NUMBER 
-  FUEL  MEMBER  NUMBER 


.  IPACC.. 
...  IPCN._, 


-  DATA  SOURCE  CODE  FOR  GROUP  SERIAL  NUMBER  ISFN 

serial' NUMBER  ISFN 

^..NUMBER  OF  contributor  COUNT.  CHANGESJF.OLLOHINCL. 
(HEHAINOER  OF  -G-  CARO) 

-  NWBER  OF  CONTRIBUTOR  NAME  CHANGESJF.ULLOW.ING _ 

IStX  PER  -H>  CARO) 

■«  NUMBgR  OF  PAIRS  IN  CONTRIBUTOR  COUNT  LlSt__ 


_ |PI  . 

IPIN 
ISAN- 


FOLLOWING  ON  >APE 
_:uJIUMB£R  OF  CAROS  FOLLOWING  WITH  GQRTIUB.UT.OR.  . 

COUNT  CHANGES  l-G*  CAROS) 

_2JWMB£R_0F  .CROUPS  ON  TAPE  2  .B£FQBE_RD.DJ  UONS_ 


-  NUMBER  OF  GROUPS  ON  TAPE 
JftAIA_GROUB_  SERI  AL_.NUMU.ER_ 


AFTER  ADDITIONS 
IFROM  lAPE  2i 


trP  -  CODE  INTEGER  FUR  INITIAL  TAPE  PREPARATION 

jl(K,n _ .<t_SF£Cl£S  CONTRIBUTOR  CODE  NUMBER _ _ _ 

JGNLtK)  •  CONTRIBUTOR  CODE  NUMBER 

XI  . . -.NUMBER  OF  STRUCTURAL  COM.TR I BUTORS, CONSIDERED _ 

LACNILtR)  '  CONTRIBUTOR  CUDE  NUMBER  (FROM  CAROS) 

.  _.LFCCCNWl _ S_CQNTR1BU.I0R  .COO.E_NUMOElL_|Fi(OMJA«a _ 

LFMH(L)  -  DATA  GROUP  SERIAL  NUMBER  IFKOM  CARDS) 


TAPE 
TAPE 
TAPE* 
,  TAPE 
TAPE 

_ TAPE 

TAPE 

_ ^TAPE 

TAPE 

TAPE 

TAPE 

.TAPE 

TAPE 

TAPE 

TAPE 
..  .tape 
TAPE 


1 

2 

3 

_  4 

■5 

6 

7 

_  B 
9 

^10 

11 

-  12 
13 
U 
13 
—  IB 

17 

18 

19 

20 
21 

-22 

23 

24 
23 
26 
27 

—28 

29 

JO 

31 

32 


33 

36 

37 


41 

42 

43 

44 
43 

_A6 

47 

4B 

49 


LIPACCIL) 


-  NUMBER  OF  CONTRIBUTOR  COUNT  CHANGES.  FQR  DATA 
GROUP  SERIAL  NUMBER  LFMNIl) 

-  COUNTING  INTEGER  FOR  LIMES  PRINTED  PER  PACE _ 

-  r'un  number  of  data  group 

RUN  NUMBER  CHECK 


NRUN 

nrunt  _  _  _ 

subroutine  tape'  lNRUN,ICTL»AltA2t  IIA,llB»IlCiI2,I3»0A»i4»jt,BI» 

1  lAOOi^fiUf  TE*  LAST)  . . 

COMMON  INT,  UEC*  1C»  All*  J2,  J3»  N1N»  NEX 

DIMENSION  INTIIO).  OECIIO)  . . 

biMENSIUN  Al(2O),A2l2O)tBltl00»20)»OAI6» 

1  , 1^(20) tJt(100.20)«ACCL(20, 90) »AGNLi50)»D(6)iDNNI6)*FGCLI40l».  „ 

2  FGCJL(2UO)iGNL(200)*JGNLI30),LACN120»90)»LFNNI20)*LFCCCNI40)» 

3  LIPACC(20),  I1A<20)«  11B(20}*  I1CI20)  .  , 

H  «  4 

_ WHITE  OUTPUT  lAPE  3,  >000, _ 


OU.TE 

__  WRITE  OUTPUT  TAPE  3» 
9  NKUN 

REWIND  6 
10  CALL  INPUT 
NEX  «  NEX 


IFI  NEX  -  6)  600»  11»  3} 
F  CARD  ASSIGNMENTS 
Jll”  ■ 


11  NMUNT 
_ ^ICANJL 


IN^Td) 


IPFM  •  |Nri2) 

IPCN  -  INflJ)  _ 

ITP  »  INT(4) 

LAST  ■_|Nr(3)_ 


IF  (NRUN.  -  NRUNT)  .2.9.«.21|29 
OPTION  CONTROL  .INTEGER  CHECK 


21  IP  I  ICTL  -  ICANT)  29.22.29_ 


TESTS  FOR  executable  CONOUJONS.  . 


22  !F  (ITP  ♦  .7)  25,  26,25.. 

23  IF  (IPFM)  26,26,32  ■ 


TEST  FOR  CROUPS  READY  TO  BE  ADDED 


26  IF  (IPCN) 

27  IF  iuno) 


27,27,40 

30,30,40 


OI.A.ONOSTXC_ 


-  INITIAL  TAPE__PJIEPAIimD)T- 


26  WRITE  OUTPUT  TAPE  3,  3002 


TAPt 

30 

TAPE 

31 

.52 

TAPE 

33 

TAPE 

54 

TAPE 

35 

TAPE 

56 

“TAPE 

57 

TAPE 

36 

“  TAPE 

59 

TAPE 

60 

TAPE 

61 

TAPE 

62 

“TAPE 

"63 

TAPE 

64 

TAPE 

63 

TAPE 

66 

TAPE 

67 

TAPE 

6B 

“TAPE 

69 

TAPE 

70 

'TAPE  " 

71 

Tape 

72 

‘TAPE 

73 

TAPE 

74 

TAPE 

fi 

TAPE 

76 

"TAPE 

77 

TAPE 

7B 

TAPE  " 

79 

TAPE 

60 

TAPE 

61 

TAPE 

82 

TAPE 

63 

TAPE 

64 

TAPE 

63 

TAPE 

86 

TAPE 

87 

TAPE 

68 

TAPE 

89 

TAPE 

90 

TAPE 

91 

TAPE 

92 

TAPE 

93 

TAPE 

94 

TAPE 

95 

TAPE 

96 

TAPE 

97 

TAPE 

98 

TAPE  ‘ 

99 

TAPE 

too 

TAPE 

101 

TAPE 

102 

TAPE 

1U3 

TAPE 

104 

^TAPC  ■ 

lUS 

TAPL 

106 
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_ iHHQW  finj  ~  DATA  OUT  Of  DROta 

UNITE  OUfWiT  _ _  .  _ 

9  NHUN»M(UNr«tCTL«lCANr 

UO  TO  AOO  _ _  _ 


._30  UNITE  .  OUTPUT.  TAM  >.  1004  . 

QO  TO  495 
32  I  ■  0 
»21  L  «  L  ♦  1 

-  <;o  TO  10 _ 


_.33  IFI  Hit  ~  Tl  600«  330#  400 _ 

330  NUMCON  •  1C  «  7000 

_ IFI  NUUCOWI  332.  331,  332  _ 

331  4.FNNIL)  «  INTIl) 

__  .IFACC  •  IMTI2I _ 

60  TO  10 

_ 532  IFI  NNUNt  -  Jill  AOO,  333.  OOP 

333  LIMIT  •  lO«NUMCdN 

_ IIA  «JJHIL.t  y _ 

tfUIMtf  •  IFACCI  335»  335.  334 
_334.. LIMIT  •,  IFACC.  _ _ 

335  M2>  0 

„00  336  A  •_I60  .  LIMIT _ 

M2«  M2«  1 

_ LACH.I  L.KI  .  1MTIM21 _ 

336  ACCLIL.Rl  •  UECCN2) 

IFC  LIMIT.-  IFACC),10.,_33*_10.__ 
35  LIFACCIL)  «  IPACC 

IFI  L  -  IPFM  »  32j|t-A0xJkfl _ 


40  NEAO  IMPUT  TAP.t..6*_iO0# _ 

NEAO  INPUT  TAPE  6«  2001*  IFI 


NEAOV  SCRATCH  TAPE  0 


COMPUTATIONS  FOR  ALTERATIONS  AND  ADDITIONS  TO  DATA  OROUPS  ON  TAPE 
00  90  l•i»IP|  ”  ' 


NEAO  EACH  6H0UP  STORED  ON  TAPE  6 

READ  INPUT  |AP€~^"2002i  ^ 

9  l$FN>FNl»FN2»IFNCtlFN6»IFNN|INOStlNEC»IOIKI»K*2»AI» 

_ X _ _  -  IPPi;C.iLP6CCNIRI,F0CLlRl.K»l>IPFGCI 

,  TES.T..F,OR  PRQPEIL.0AT,A.5RpijlL.SfAIAL.  NUMRER  SEQUENCE  ON  TAPE _ 

IF  II  -  UFN)  4U45*41 _ _  _ _ _ 

41  REUINO  6 
WEUINO  ■ 

TEST  TO  LIMIT  THE  NUMOER  OF  LINES  PRINTED  PER  PAPE  OP  OUTPUT 

IFIM  -52)  42,43,43__  _ ......  .  _ _ 

43  UNITE  OUTPUT  TAPE  3»  3020 


42  UNITE  OUTPUT  TAPE  3*  3005* 

_  9  ISFN.I _ : _ 

60  TO  500 

"C  OYPASS  OF  TESTS  IF  NO  CAROS  To  Of  READ 

.C _  ... _ ; _ 

45  IF  IIPFM)  OOtOO.AA 

’c  COMPU r A riONS  for  addition  of  infornation  to  tape  i 

C  -  ■  - 

46  DO  47  L-UIPFM 

"C  “  IE'S!  POfTALURAT  IONS  TO  OROUP  I S'FN~6N~tAPE  iY  DATA  PRON  -O-  CARO" 

C _ _ _ _ 

IF  It  -  LFNNUI)  4T»50«4T 

4T  CONTINUE - 

60  TO  «0 

.50  00  9l..3>l«200  _ _ _ _ 

C 

C  lERO  COUNT.  .Ll>l_St.Ti»&«C _ 

C 

51  F6CJLIJ)  »  0.0 _ 

DO  52  R-1»IPF0C 

C„ _ _  .  _ _ : _ 

C  STORE  TAPE  VALUES  OF  THE  NUMOER  OF  STRUCTURAL  CONTRIOUTORS  PER 

C  .  _  .  MOLCCULU«LQ&_CAi)i.WWJIW.TPR.WTMIN„TI«JLPJ^JJ.tl^^ 

C 

_ J  ..•.JljTfiCIliLtKl _ 

52  F6CJL13)  •  F6CLU) 

C  AOOITION  OF  NEU  VALUES'OP' THE  NUM  OP  STRUCfURAV'l^fRrOUTldm^ 

C  PER  NULECUi^JEPlUfMOSE  C.ONTNIRUTORS  tlST.SO._Q(LiNJLrDr.  CAOO .. 

C 

_ I PACC  ■  LiPACCtD  -  - ^ - 

00  53  Ml*  IPACC 

J  •  LACN(l*R|  .  . . . 


TAPE  107 
'  TAPE  too 

_ TAPE  109 

TAPE  Ub' 
_  TAPE  111 

'tape  112 

_  TAPE  113 

TAPE  114 

TAPE _ 115 

TAPE  116 

_  TAPE  ITT 
TAPE  lie 

TAPE  IIV 

TAPE  120 

_ ^TAPE _ 121 

TAPE  122 

_  TAPE  123 

TAPE  124 

_  TAPe_  125 

'tape  126 

TAPE _ 127 

TAPE  125 

_ iTAPE  129 

TAPE  130 

_ _ TAPE„  lil 

TAPE  132 

_ TAPE _ 133 

TAPE  134 

_  TAPE  _  135 

TAPE  136 

TAPE  137 

TAPE  136 

TAPE _ 139 

TAPE  140 

_ TAPE  141 

TAPE  142 

_  TAPE  _143 

"tape  144 

_ TAPE _ 145 

-  TAPE  146 

_  TAPE  147 

TAPE  148 

_ TAPE  149 

TAPE  150 

_ TAPE _ LSI 

TAPE  152 

_ TAPt  153 

TAPI  ■  154 
TAPE  155 

TAPI  156 

_ TAPI _  157 

TAPE  158 

.. .  TAPE  159 

TAPE  160 

TAPE  161 

TAPE _ 162 

T'Ai^E  163 

_ TAPE  164 

TAPE  165 

_ TAPE  166 

TAPE  '  167 

_ TAPE _ ^166 

""■^TAPE  169 
TAPE  170 

"TAPE  “*ITI 

_ _ TAPE  _  172 

TAPE  173 

TAPE  174 

TAPE  175 

_ TAPE  _  176 

TAPE  177 

_  TAPE  _178 

“tape"  179 
TAPE  180 

TAPE  181 

_ TAPE _ 162 

TAPE  163 

_ ^TAPE  _ 164 

TAPE  165 

_ TAPE _ 166 

TAPE  167* 

_ TAPE  _  188 

TAPE  169 

_  TAPE  _  190 

TAPE  191 

TAPE _ 192 

TAPE  193 

_ TAPE  194 

TAPE  195 

_  TAPE  _  196 

TAPE  197 

_ TAPE _ 196 

TAPE  199 

TAPE  _  200 
TAPE  201 

_ TAPE  202 

TAPE  203 

_ TAPt _ 204 

TAPE  205 

_ .  TAPE  206 

TAPE  207 

. TAPE  208 
TAPE  209 

_ TAPE.^  210 

TAPI  211 

TAPE  212 
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SUBROUTINE  TARE 


ICQNTINUEOI 


MASTER  LiaMRV  TARE  NOOIFICATICM 


9)  RGCJLUk 
K  -  0 


ACCL.UiKI 


C 

c... . 

c 


LIS't*  CONTRIBUTOR  CODE  NUMBERS  AND  NUMBER  OR  CONTRIBUTORS  RIR 
MOLECULE..  IQ-AC  JJRJTMO  P*i  TARE.,9  _ _ 


no  55  J«lt200  _ 

IR  (RGCJLtJI)  94,55*54 
94  K  •  K_*  L. 


LRGCCNiR)  •  J 
FGCLtKI  •  FCC^LIv'l 
95  CONTINUE 
IRFGC  •  K 


TEST  TO  LIMIT  THE  NUMBER^ OR  LINES  RRIMTEO  RER  RAfiB  OR  QUTRUT 


_  IFIM  ♦.  1 
97  URITE  OUTPUT  TARE  3« 
M«0  . 

54  H  ■  M  «  lPFGC/4  *6 


PRINTOUT  OF  ALTERED  LISTS 


WRITE  OUTPUT  TARE  3«  3006 

WHITE  0OTPUT_TAP£_3t_30OTt _ 

9  ISFNt  UF6CCNIKI|FGCL{K)»K»lilRF6C) 


TARE 
TARE 
TAPE 
TAPE" 
TAPE 
■  TAPE 
TAPE 
TARE 
TAPE_ 
"tape 

TAPE 

TAPE 

TAPE 

TAPE 

TAPE 

"tape 

TAPE- 

TAPE 

TAPE. 

TAPE 

__TAPe 

TAPE 

_TAPE 

TAPE 

TAPE 

TAPE 

TAPE 


TAPE 
TAPE 
TAPE 

ISPM.Fliil«FM3.IPMC»IPNG»iFNM.INDSilMEC.IDIKI  _ TAPE 


C  INTERIM  STORAGE  ON  TARE  B 

*’e0  WRITE  OUTPUT  TARE“b»~20”o2'“ 

_  9  .  _ _  .  . . 

1  IPFGCt I LFGCCNI K ) (FCCL I K ) »K« 1 . I PFGC I 

_5.O_CL0j!UJLH.U.I _ — 

C 

C  END  OF.  ALTERATIONS  AND  ADDITIONS  TO  DATA  GROUPUMLlARE _ u.. _ 

C 

C  ALTERATIONS  OR.  ADDITIONS  TO  tONIRIBMT.OR  NAME.  L.l^ _ 

C 

C  READ  CONTRIBUTOR  NA¥eS  AND  OEPENOEnT^VAiTiABLE  NAME  LIST  PROM  TAR^ 

READ  INPUT  .T.AP.E.6»_j20D3Li _  ...  _ - 

9  (CNL(J)*J«1»200) 

_  READ  INPUT  TAP^  6^2003.*  ...  _  _ _ 

9  (bNHlKl|K-t*6) 

_ REWINft_.R__ _  ,  ■  ,  ..  .. 


aRW|Nn  8 

TEST  FOR  AIT£RAT10N~0R  ADDITION" 
IF  IIPCNI  I0U*10U*95' 


H  CARO  READING 


99  READ  INPUT  TAPE  2»  lOOt* 

9  UGNL|Kl|AGNL<K)tX*l«lPCN) 


■  TEST  TO  LIMIT  THE  NUMBER  OF  LINES"PirnitED"PER  RAGE  OR  OUTPUt 


1F(M  »  lPCN/7  '  52)  97*9B»98 

99  WRITE  OUTPUT  TAPE  3#  351.0 _ _ 

M>0 

97  M«M  ♦  lPCN/7  *6 


.  PRINTOUT  .OF  .AI.T5.RATJPN.S  .QR  .AOD.IT10NS_JO  CONTR|BUTPH  NAMg 

WRITE  OUTPUT  TAPE  .3*.  3008  _  .  _ _ _ 

WRITE  OUTPUT  TAPE  3f  3009» 

IJGNHK).AGNHR).K»1^IPCNI _ 


00  96  K«t»IPCN 
.  J  .•.JCNLAKI _ 


C_  AOO  ALTEREO  OR  HEM  NAMES  TO  LIST 
C 

94_CNL<JI  »  AGNLIR) _ 


.  REVISE  GROUP  COUNT 


_ LOO_IPIN.-  IPl  ♦  lAOO 


PREPARE  CORRECTED  TARE  4 


.WRITE  .OU.TRU.1_WRe_.4i_2ftSL0_ 


WRITE  OUTPUT  TAPE  6»  2001f 

9 _ _IPIN _ 

C~'  (A)  CORRECTED  GROUPS  PREVIOUSLY  ON  TAPE 

C  _  _ 

00  UO  I  •  ItIPt 
READ  _iNPS»f  TAPE  8t  1092« 


TAPE 
TAPE_ 
TAPE 
-TAPE 
TAPE 
TAPE 
TAPE 
_TAPE 
TAPE 
TAPE 
TAPE 
TAPE 
TAPE 
_TAPE_ 
TAPE 
TAPE 
TAPE 
TAPE 
TAPE 
_TAPE 
TAPE  ■ 
TAPE 
TAPE 
TAPE 
TAPE 

""Yape" 

TAPS 
"  TAPI 
_TAP| 
TAPE 
__TAP| 
TAPE 
_TAPE 
TAPE 
TAPE 
TAPE 
_rAP£ 
TAPE 
__TAPE 
TAPE 
_TAPt 
TAPE 
_tAPE_ 
TAPE 
_TAPE 
TAPE 
_1APE 
TAPE 
—JAPE 
TAPE 
_  T*R£ 
TAPE 
_  TAPE 
TAPE 
_ TAPE 


ISFN«FNl»FN2»iF'NCtrFNG|lFNN|IND$ilNEC|IDIR>tMI»B)» 

I  _  IPFGC,IL.FGCCN.I5.U.FG.g.L.W)tR-l.ilRe«.> _ 

ilO  WRITE  OUTPUT  TAPE  4t  2002. 

9  : _ lSFM.fMI>fM2.lfN(^lJ:lC^.iFNM^N0iltlNEC.iPIRI»K*2iBlt 

1  |PFGC*ILFGCCN|K)tF6CLIRI|K>lf IPFGC) 

IF  IIADOiaiSaiStUI  .  _ 


TAPE 

_ TAPE 

TAPE 

_ TAPE 

TAPE 

_ ^TAPE 

TAPE 
_  TAPE 
TAPE 


TEST  rOJl.IM.IT  T_ME_NUMe£R  OF .L  INES.  PRINTED  PER  PAGE  OF  OOTRUt _ TAPE 

TAPE 

IFIM  -  961  ill.H4.il4 _ ^lAPE. 


114  WRITE  OUTPUT  TAPE  3t  3020 
M"0  .  _ 

113  M«H*2 


DIAGNOSTIC 


•  NO  NEN  DATA  GROUPS  ADDED  TO  TARE 


WRITE  OUTPUT  TARE  3f  3011 
GO  TO  14U  _ 


TAPE 
TAPE 
TAPE 
Tape 
TAPE 
_ TAPE 

TAPE 

TAPE 


213 

214 

215 
'214 

217 

2tS 

219 

220 
221 

’  222' 

223 

224 

225 

226 

227 

228 

229 

230 

231 
'  232 
-233 

234 

235 

236 

237 
23tt 

_239 

240 

241 

242 

243 

244 

245 
244 
247 
24B 

249 

250 

-251 

252 

253 

254 

255 

256 

257 
256 

259 

260 
261 
262 

263 

264 

265 

266 
267 

_266 

269 

270 

271 

272 

273 
_274 

275 

276 

277 
..  278 

279 

281 

..262 

263 

_264 

285 

_286 

267 

288 

289 

290 

291 
—292 

293 

294 

295 

296 

297 
—298 

299 

300 

301 

302 

303 
—304 

305 

306 

307 

308 

309 
—  310 

311 

312 

313 

314 

315 

316 

317 
316 


60 


SUBROUTINE  TARE 


-  ADDITION  DF  NEW  CAOURS  TO  TARE 


C  TEST  TO  LIHIT  THE  |H>NKRjPR  ilNES  RRINTED  Ri^^MHLfiE-fiUtW - 

C  . 

IIT  IR(H  -  S2I  11B»11B|114 -  - - - 

116  MRITE  OUTRU1  TARE  S*  S020 

_  N-0  - -  - 

118 

C  RAINTOUT  OF  NEM  MDURS  ADOEO  TO  TARE 

URtTE  OUTRUT  TARE  S«  SOIE 

__L2.Q_oo_UftJiLiJUfifi - - - 

ISFN  •  IRl  *  I 

R1  ■  141 U. _ _ _ 

1F(M  •  SR)  12y*l2Ba26 

128  WRITE  OUTRUT  TARE_2*_3020 _ _ 

WRITE  OUTRUT  TARE  )t  3012 

_ M-P _ _ _ 

129  H>M*I 

_  WRITE  OUTRUT  TARE_3,..?D1.3i. _ _ — — 

9  AI(n*A211)t  ISFN  •llAlt)tllBII)*IICn)tl2inil3nit 

_  I  _ DA  1 2j.l  lj..-OA  1 4  t.l.it,DA  L3.»  J.U.  P*.l « t  U  I  P*  t  tliJ - 

130  WRITE  OUTRUT  TARE  6*  2002«  _ 

_ 9 _ lSfN.AUI),A21ll.llAni.llEUjilLCjni  2!  tUlIb- 

1  lDA(R«l)»K*2*6)*Rl*4JIIR*l)fBllRin«K*TiKII 

“c"’  TEST  TO  UNIT  THfiliUHOER  OF  LIMES  RRINTED  RER  RAOf  OF  OUTRUT 


tF(M«l4IL)/6  •  SO)  131,132|132 
H-O  ‘ 

131  _  . 

WRITE  OUTRUT  TARE  3#  3014 

00  13$  l•l».IAOD _ 

ISFN  •  IRI  ♦  I 

_ 11  «  l^tU - 


_|  P I H*  I4I.I  )/6_.-_»4JLJLi2UM».UA. _ 

134  WHITE  OUTPUT  TARE  3*  3020 

_ itRl.TX-flliUJiUT  T4RE  it  >014 

M>6 

133  MH4|4II)/6  ♦2.. _ _ _ 

139  WRITE  OUTRUT  TARE  3t  3007* 

9  ISFN, _ I  iKRAtJIlJ — 

■  c  _  icj _ :Lj>j»xufliLa€-tJaam»^^ 

"c 

140  WRITE  OUTRUT  TARE  .6, ..2003,. .  . . . . 

9  CGNLtJ)»J*l*200)  .  .  _ 

WRITE  OUTPUT  TARE  6*  2003* 

_ y. _ _ 

ENO  FILE  6 

.  REWIND  6  ..  . . . . . 

REWIND  8 

CO  TO  SOO  . . . . . 


ORERATIONS  FJW.JNJXIAt„NAKf_MR_ .OF.„TAR|  2  _ 

19S  IF  lIPCM).  490,490,.2l0  .  _ _ 


.JIIO  READ  INPUT  TARE  2.  2000 _ 

C 

C  N  CARO  READING _ 


IJCNllR)*AGNllR)*R«ltlRCN) 
OERENOENT  VARIABLE  NANlPriTT'CAHO'' REA^N^ 


READ  INPUT  TARE  2*  1002* 

1 _ I0NN«R).R»U6) 

DO  219  3«1,200 


C  BLANK  OUT  ENTIRE  CONTRIBUTOR  NAME  FIELD 

'219  CNLU)'>  BLK  .  . . 

_ 00  220  __K«UIRCW _ ^ _ 

J  •  JGNLIK) 

C  STORE  CROUPS  FR0N"-N^~CARD  IN  CmTRIBUtOR~NAiil  RIBLO 

220  CNLUI  -  AGNtfK)  ’ 

_ IP  MAOQI  490.49Q.23Q _ 

C 

C  .  PRINTOUT  OF  INIT.IAL_0AJJ1_EP.R_WFE . . 

e 

230  WRITE  OUTRUT  TARE_3,_2O00: _  _ _ 

WRITE  OUTPUT  TARE  3*  3010 

_ WRITE  OUTPUT  TARE  9.  3009. _ ■ 

9  l2GlltlKI*AGNLU)*RM*IRCNl 

WRITE  OUTRUT  TARE  3,  .3019,  _ _ 

9  IONMIR)«K«t*6),BlK 

C  RRERARf  TARE  6 

WRlTE'OUiFUrTARC  •*  2000 

WRITE  OUTRUT  TARE  6,  2001,  _  ..  .  _ 


TARE  319 
TARE  320 
TARE_  321 
TARE  322 

_  TARE  323 

TARE  324 

_ TARE  32$ 

TARE  326 
TARE  327 
*“TARE  328 
TARE  329 
TARE •  330 

TARE  331 
TARE  332 

TARE _ 333 

“^TARE  334 
TARE  33$ 
TAPE  336 
_  _TAR6  337 
TAPE  338 

_ ^TARE _ ^339 

TARE  340 
TARE  __  341 
TARE  342 

_ :.TAR£  343 

TARE  344 

__^TARE _ ;.345 

TAPE  346 
_  TARE  347 
TAPE  348 
TAPE  349 
TAPE  390 

_ TAPE _ iSl 

TAPE  3S2 

_ TARE  393 

TAPE  354 
__  TARE  355 
TARE  396 

TARfc _ 357 

TAPE  3S8 
„  TARE  359 

TARE  360 
TARE  361 
TAPE  362 

_ TARE _ 363. 

TARE  364 
TARE  36$ 
tAPE  366 
TAPE  367 
TAPE  368 

_ Tape  .  369 

TAPE  370 
;  TAPE  371 

.  TARE  372 

TARE  373 
TARE_  374 
TAPE  ‘379 

_  TAPE  376 

’TAPE  377 

_  T3PE  378 

TAPE  379 

_ TAPE  380 

TAPE  381 
_  TAPE  _  382 

TAPE  383 
TARE  _  3B4 
TAPE  389 

__TARE _ 386 

TAPE  387 

_ TARE  388 

TAPE  389 
_  TAPE  390 

TAPE  “'391 

_ TAPE _ ^392 

TAPE  393' 

_ TAPE  394 

TAPE  399 
tape  _  396 
tape'  *397 
TAPE  398 
TARE  399" 
_  TARE  400 

TAPE  401 
TARE  402 
TAPE  403 

_ TAPE _ 404 

TAPE  405“ 

_ TAPE  406 

TAPE  407 

_ TAPE  408 

TARb  409 

Tape _ 4to 

TAPE  411 
_  Tape  412 

TARE  413 
TARE  414 
TAPE  419 

_ TAPE _ 416 

TARE  417 
TARE  4U 
TARE  419 
TARE  420 
TARL  421 

_ TARB  _  422 

TARE  423 
TARb  424 
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SUBROUTINE  TARE 


.(CONTINUED) 


nas7er  limabv  tare  wfc-'^MCUIAN 


H  lAOD 

00  240  I>*)•1A00 

.Kl  •  IKIl- 


240  WkITE  OUTRUT  TAR£  b,  ZOOfi 

9  1  iLAi<iu*2ii)«iM(n»iie(nii  _ 

1  IOA|K*tUM2«6f*K(*lJllKfll|BllKil)|i(«liKI) 

WRITE  OUTRUT  TAR^  6»  200?,  _ 

9  (6NL(J)«J«lf2bO) 

_  WRITE  OUTRUT  TAPE  6.  2003. 


(ONNlK),K*UA) 

END  Fite  6  _ 

WRITE  OUt^T  TARE  ?•  3016 
REWIND  6  _ 

ENO  OF  INITlAi  RRERARATION  OF  TARE  6 


..  GO  JO  50JL. 


C 

C  _  _  ERROR  RRIMT__ 

C 

400J<R1TE  OUTRUT  TARE  3.  3017. 
9  IRCNtlAOO 


>  NO  OATA  FOR  INITIAL  TARE  IHieRABATION 


495  ..REWIND  6 
500  RETURN 
_^00  NEX  ■  2 


WRITE  OUTRUT  TARE  3» 
CQ  TQ  SOO _ 


4000f  fC*  Jll 


FORMAT  .ATJHe.N«I.T.i__ 


C 

C-. 

C 

tool  .FORMAT _ _ _ ^ _ 

9  I t2X  •I4,2XtA6»14,2X»AAtl4»2X»A6.l4»2XiA4»14*2X*A6) 

1002  FORMAT  1  12X,  7A6  I _ 


2000  FORMATIIHI,1IX60H 

2001  FORMAT  . 

9  IlHO.  4.141.  _ _ 

2002  FORMAT 

9  IlHQ  l4,2A4.5l4,IR5et2.4.l4/ilM  lA.ElE.A. 14.6 12.4.16. £12^^14. 


IEI2.4,14,E12,4,14,612«4)) 
2003  FURHAT 


UNO  Vl4X«A4)/(lH  4X,A4,4X,A4>4X,A4,4X,A4»6X|A4»4X|A6,4X|A6|6X» 
1A4»4X*A4.)) _ 

3000  FORMAT 

ROUTINE  2014  MODIFICATION  I  ’  ’  OAfE  2A6  ) 

3001  FORMAT  _ _ _ -  - -  - 

9  nHa*IOXt  AttHTARE  2  WRITING  INFORMATION  *  RUN  14  ) 

3002  FORMAT  _ _ _  .  >  _ -  _ 

9  (IH0»15X»  24H  INITIAl  TARE  RRERARATION  ) 

..  30.03  .  FORMAT.  .. 


9  (  IHUt  14X,30M0ATA  OUT  OF  ORDER  •  RUN  NUMBER  t4» 

1  14,X,  _  ,20HMlSRtAtE.D  RUN  NUMBER. 

2  /IMO  14X,22h6rTION  CONTROL  INTEGER  .I4» _ 

3  22X,2lHCONTROL  INTEGER  CHECK  14) 

_3QP4  FORMAT _ 

9  (iHOVlOX, 

.3005  FORMAT  ..  . 

9  UHOtlOX* 
lER  •  14  , 

3004  FORMAT 

_ 9j.lH0,.lfi.Xj, 

I  IN  ,|4H 


20HN0  ACTION  CALLED  FOR  ) 

57MGRbu>  SERIAL  NUMBER  INCONSISTBNCV  ^ 
5X»  liWLOCATION  ■  14.)  ,  .  .  _ 


34HAtTERATIONS  TO  DATA  GROUF^ON  TAF^j)^^^^ 


OATA  CROUP  .4X,llHCONrRIBUTOR(4X*llHCONTRIBUTOR»4X»liH 
2CONrKI0UTOR,6Xf llMCONTRIBUTOR ,4X»llHC0NTRIBUT0R,4Xt UHCONTRIBUTOR  _ 
3/  IH  ,I17H  SERIAL  ’  CODE  COUNT  /  CODE  COUNT  /  CODE 

4  COUNT  /  COOE  COUNT  t  CODE  COUNT  /  CODE _ COONT_/_ 

5/  IN  ,tlttN  NUMBER  NUMBER  MOLECULE  NUMBER  MOLECULE  NUMBE 

6R _ MOLECULE  NUMBER  MOLECULE  NUMBER  MOLECULE  NUMBER  MOLECULE  /) 


3007  FORMAT 

_  9  UHU.  IlOjJJt,  6n4,2X.F7#3>2X)  /(IM  ,  17X,  I4,2X,F7i3|2X»  I4,.«,FT.3„ 

1,2X,I4,2X,F7V3,2X,16,2X,F7,3,2X,I4,2X,F7.3,2X,14,2X,F7.3  )) 

3006  FORMAT  _  _  _ _ 

9  (iHOflOX,  32HCHANGES  IN  CONTRIBUTOR  NAME  LIST  ) 

3009  FORMAT _ 

9  I  IMO,  t4H  CONTK|BUTOR«4X,UHCONTRIBUrOR,6X,llHCONTRIBUTOR|6X|llH 

lCOWTRIBUTOk.fcK.HHCONTRlBUT0R.6X,llHC0NTRIBUT0R,6  _ 

2/  IK  ,  5H  CODE  l3X*4HCd0E  13X,4HC00E  *3X,4HC0DE  13X,4HC0DE  13X| 

_ 3  _  4HCU0E  13X,4HC0DE  /  _  _ 

4  IH  ,ll7HNUMeER  NAME  NUMBER  NAME  NUMBER  NAME  NUMBE 
9R  NAME  NUMBER  NAME  NUMBER  NAME  NUMBER  NAME  / 


4/llH  •  t4,2X,A6t3X,Uf2XtA4»3X,16,2X,A4-^x7lAf2X>A4>3X»14»2XiA6,3X 

7,I4,2X,A4,3X,I4»^X»A4.J  I  ...  _ ! _ 

3010  FORMAT 

9  UHO.lQXt  29MINITIAL  CONTRIBUTOR  NAME  LIST  ) _ _  __ 

soil  FORMAT 

_  9  ( IHO, IPX.  32HN0  NEW  DATA  CROUPS  ADOEO  TO  TAPE  )  _ 

3012  F'dRMAT 

9  (IHOflOX^  25HDATA  CflOUPS  ADDED  .TO  TAPE  ./ _ _ 

I  IH  •  lOH  FUEL  NAME  IvX, I2HC0DE  NUMBERS  22X| 14HST0ICH1ONETRIC  IIX 
2,13HMAXIMUM  SPEED  4X»  IIHEOUIVALENCE  /IH  4.1>Rt.L0.iM$J.SJAI<..CLASX-C„ 
3R0UP  MEMBER  OATA  EXPERIMENT  FLAME  FUEL  FLAME 

_ 4 _ f.Ufct _ BRritt.RT  l- 


5  IH  ,'39X,T7HS0URCE  CONDITIONS  SPEED  CONCENTRATION  SPEED 
4  CONCENTRATION  MAXIMUM  /|H  •57Xt61H(CN../SEC)  ..LMQL.6.CUjk|U/CC)„(.C.Ni. 
7/SEC)  IHULECULES/CC)  FLAME  SPEED  /) 

3013  FORMAT  ........  -  .  _ 

9  IlH  ■  2A4«  17*214*  317*  4X*  F9,4»  2Xt IPE I2.4»2X(  0PF9.4*2X»lPEi2« 
_ 24,  3X  *  _0PFB,.4_) _ 

3014  FORMAT 


9  nHO*lOX*  32HA001TI0NS  TO  OATA  CROUPS  ON  TAPE...  _ _ 

1  IH  flAH  DATA  GROUP  *4Xt  I  IHCONIRIBUTORiAX*  UHC0NTRIBUT0R*4X*  I  IH 
2CONrRIOUfnR»4X*l|HCONTRiBUTOK*4X,llMCONrRlBUTOR|6X|llNCON.rR.IBUTOR.  . 
J/  IH  *117M  SERIAL  CODE  COUNT  /  COOE  COUNT  /  COOE 

4  COUNT.  /  COOfc  __CPUWT  ^^_COP.E C_O.UNT,/_  _C»JL_QflH«.Ly _ 

'5/  IH  (llBN  NUMBER  NUMBER  MOLECULE  NUMBER  MOLECULE  NUMBE 

4R  MOLECULE  NUMBER  MOLECULE  NUMBER  MOLECULE  NUMBE.R.  MOLECULE  /). 


TAPE 

423 

TAPE 

424 

tAPE 

427 

TAPE 

42B 

TAPE 

429 

TAPE 

430 

TAPE 

431 

TAPE 

432 

TAPE  433 

TAPE 

434 

TAPE 

433 

TAPE 

434 

TAPE 

437 

TAPE 

436 

TAPE 

439 

TAPE 

440 

TAPE 

441 

TAPE 

442 

TAPE 

443 

TAPE 

444 

TAPE 

445 

TAPE 

444 

TAPE 

447 

TAPE 

448 

TAPE 

449 

TAPE 

450 

TAPE 

451 

TAPE 

452 

TAPE 

453 

TAPE 

454 

TAPE 

455 

TAPE 

454 

TAPE 

457 

TAPE 

456 

TAPE 

459 

TAPE 

440 

TAPE 

441 

TAPE 

442 

TAPE 

443 

TAPE 

444 

TAPE 

465 

TAPE 

466 

TAPE 

447 

TAPE 

448 

TAPE  469 

TAPE 

470 

TAPE 

.  A7l 

TAPE 

472 

TAPE 

473 

TAPE 

474 

TAPE 

475 

TAPE 

474 

TAPk 

477 

TAPE 

47B 

TAPE 

479 

TAPE 

4B0 

TAPE 

4BI 

TAPE 

482 

TAPE 

483 

TAPE 

4B4 

TAPE 

485 

TAPE 

484 

“tape  4B7 

TAPE 

488 

tape 

489 

TAPE 

490 

TAPE 

491 

TAPE  492 

TAPE 

493 

TAPE 

494 

TAPE 

495 

TAPE 

496 

TAPE 

497 

TAPE 

498 

TAPE 

499 

TAPE 

500 

TAPE 

301 

TAPE 

502 

TAPE 

503 

Tape 

504 

TAPE 

505 

TAPE 

506 

TAPE 

507 

TAPE 

508 

TAPE 

509 

TAPE 

510. 

TAPE 

511 

TAPE 

512 

TAPE 

513 

TAPE 

514 

TAPE 

515 

tAPE 

514 

TAPE 

517 

TAPE 

518 

TAPE 

519 

TAPE 

520 

TAPE 

521 

.TAPE 

522 

TAPE 

523 

TAPE 

524 

TAPE 

525 

TAPE 

324 

TAPE 

327 

-TAPE 

52B 

TAPE 

529 

TAPE 

530 
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tUtROUTINt  VPiCON 


iMfacmm  tmmuw  — ^nmilfinM 


SUMOUTINt  VOiCONiN  , 


RAMAViR 


OlNfNSION  tOU*Wy|Mt»  DUmvIMI 
—fltmmiflii 

I  KOurRTifai*  Aoumini 

..jCOMMON  IOUKMV*  irLUl»UlllMUS*  fMCPTi  lffft>r*T 
1  KINRUT*  NUNOCR*  RAUAV*  t*  Li  NCX.  HI 

.  iOUIVAtCMCt  I  0UIMV».-10UtWV  I 


VOfCON 

-VOiCON... 


•  lOfCMtlPfCRflt 


RtMlNft  S . 
MUMX^  •  N 
MIX  *1 
Ml  •  I 


UlMRUTf  RINRUTIt  UOUTMt  MMTRTI 


VOCCOM 
jtULsm- 
mcoH  s 

.voecoN _ 4 

VDCCON  T 

-VOKON _ • 

VOfCOM  9 

WMgnM  tn 
VOCCOM  11 

-VOECOM _ 12 

VOECON  n 

-VOECOM _ 14 

VOECOM  19 
^0£COH _ 14. 


60  T0(  1*10«  SO).  KMSt 
I  CONTlMtlE. 

"C  '  SET  UO  CHAUCfERS  SOU  LATeTTut"' 

_• _ OLUS  *.204040404040 _ 

I  AMINUS  •  404040404040 

Jl _ ftgff  -  waaaamaaaaj. 


VOCCOM  IT 

-VOECOM _ 14 

VOCCOM  19 

-VOECOM _ 20 

VOCCOM  21 


1 

COMMA  •  794040404040 

VDCCON 

29 

.4 

4 

4LANK  •  404040404040 

VOCCOM 

29 

„  C 

MCAO  ALMMAMUMCMIC  CNAM^TAlt 

9 

HEAD  IN4UT  TAME  2.1000*  UINOUTU). J*1.T2I 

VOCCOM 

27 

l«t 

60  Toao.  10>  SOI*  K9ASS 

VOCCOM 

29 

.  C  .  . 

OCCOMMOSITION  OP  19fCfiC4& 

10 

00  21  M  •  Ml*  MUMDCP 

VDCCON 

91 

.. 

NL  •  N 

VOECOM 

92 

c 

SCAACH  POD.  START  OP  NUN4E4 

VOCCOM 

99 

1A1 

_IP  IRINPUrill  .  IMtAMMI  lAA.  II.  IA» 

102 

IP  UlNPUTIt)  -  ICOMMA)  12  *  11*  12 

VDCCON 

99 

-  11 

la  !♦! 

VOCCOM 

IPI  t  -72ll0l*t01*  9 

VOECON 

97 

.C_. 

SELECT  INTEGERS 

VDCCON 

12 

LI  a  L 

VOECON 

99 

-M  a  O  _  _ 

DO  20  J*l*4 

VOECOM 

41 

KQUTPTILI  a  KtMPUTIt) 

IP(  RtNPUTtI)  -  ININUS)  122*  120*  122 

VOECON 

49 

122 

IPi  KIMPUTftI  >  IPLUlt  129.  120.  129 

120 

Mai 

VUECON 

49 

tit 

L  a  L  *  1 

!•  )♦! 

VDCCON 

47 

IPtRINPUTIII  •I4LANR»I21  .19  .121 

in 

iPUINPUTIII  MCONMA)  20  *19  *  20 

VOECON 

49 

c 

RIGHT  ADJUST  IN  PiCLO.....  . 

VOCCOM 

19 

IP(  J-61  14*21*14 

VOECOM 

91 

KOO  •  J  •  M 

WAMflii 

ft 

L2  a  LI  ♦  4  -  K 

VOECON 

94 

WASTMt 

19 

KOMTPTC  L2  )•  ROUTPTt  L9  ) 

VOCCOM 

94 

tfDPrriM 

MOO  a  4-J  ♦  N 

VDCCSH 

94 

00  14  K  •  KGO*  too 

VDCCON 

40 

w APR ns 

14 

KOUTPTI  L4I  a  0 

VOECOM 

42 

20 

CONTINUE 

VOECOM 

44 

..  _ 21 

vtoPr.nM 

lENO  a  G*  NUINKP 

VOCCOM 

44 

HRIfC  OUTPUT  TAPE  9*1000*  UOUTPTIJI*  J«|*  lEHD) 

VOECON 

44 

C 

READ  INTC6CR  LIST 

VOECON 

TO 

RCHINO  9 

VOECON 

72 

WABPnM 

*9 

90 

CALL  OCCXP 

VOECON 

74 

GO  T0{90*9I*  NCR 

yoccoM 

74 

91 

PtJI  a  PARRAVUI 

VOCCOM 

74 

1400 

P0RNATIT2AI) 

VOECOM 

40 

_ 1001 

VOPCOM 

END 

VOiCOH 

42 

65 


suMouTim  oecocf 


SUIftOUTINC  fttCOCf 

_ miNirti  irmimi  ifi 

1  KINmr*  NUMKN  Mmv«  f*  L>  NCIi  Nl 

_ MtitwAtmf  I  tntmmw  i _ _ _ 

COUIVALENCC  lAOUT^T*  KOUTETt 

_ DIMENSION  AOUTETCTEIt.  KOUTfT(7t)«.  EAEDAV.IlOX-*.JllJiMIU 

DIMENSION  IDUMMVOOIf  DUMNVOOI 

_  NEX  •  N6X  _ _ _ _ 

M  •  L 

ftn  TO  _ _ 

e  DECOMPOSITION  OP  DECIMAL  AND  IXPONENTIAL  NUNtEM 

..C _  LIMIT  OECOMPOStTlON  lQ  A  NUHiftl _ 

M  IPINUMOCP  >*)S0S*S0>»S02 

._SDa..lOEC  •  6  _ _ 

lENO  •  I2*1DEC 

_ fiO  .Til^OA . . 

SOS  I DEC  •  NUMOCP 

..lENO  -  12*  NUMDCP _ — 

SOA  lP(t  -721  SltSltSOS 

„.50S  NEX  -2  _ _ — 

60  TO  300 

Ai  ■  I  _ _  I...,. 

00  100  N*  Nl«  IKC 

_ ..Nl  ■  N  _ _ 

C  SEARCH  POP  START  OP  NUMRER 

_ SIO  IPlKlNPUrill-  JDLANR1-..52«.^1«_12 _ 

92  IPIRlNPOTItl-  ICOMMAI  94.  S3t  94 

9*  1«  I  ♦  t  -  - 

IPI  t  -  72)  910#  910#  SOS 

_ 94  Ml  •  M  _ _ 

C  STORE  NUMOERS  UP  TO  DECIMAL  POINT 

.941  IPIKINPUTII.)  -  IDECPTI-SS#.  ASt-SSj _ 

99  KOUrPTIN)  •  RINPUTII) 

_ 1-_L_*J _  I 

H«  M  ♦! 

...  .60  TO  941  ..  . _ _ : _ 

C  TEST  FOR  END  OP  NUNOER  OR  EXPONENTIAL 

_..60.  IFIRINPUTIL)  .a_IE _ J-Al#  .70«AJ _ 

41  IFIKINPUTin  -  IPLUS  I  42#  70»42 

43  IFIKlNPuni)  -  ICONMA  )  44#  40#44 

.  .44.  lFIKlNPUT(l).-.14LANR-)-4S»...40#4S _ 

C  STORE  DECIMAL  POINT  AND  NUMOERS 

_ 49  ROUTPTIMI  -  RINPUTII) _ _ 

1>  !♦! 

_ _ -  - 

60  TO  40 

C..  COMPLETC  EXPONENTIAL.^mD..TMMMItt)  D  I  WiATIOMt _ 

70  LOO  •  Ml  ♦  7 

00  71  M#  LOO  - - 

71  KOUTPTiJU  •  0 

_ M  •  Ml  KM _  ,  ,,  I  ■■ 

C _  .  STORE  €  -IN  LOCATIQN-0. _ 

KOUTPTIM!  •  IS 


72  !•  1*1 

C _  .  TEST  P.OR  .SIGN  OP- EXPONENT. _ 

73  IPIKlNPUTin-ININUS)  74#  74*74 

.  _74.IP(RlNPUTfI)-IPLUS  )  79*.  74.* 79 - 

79  KOUTPT(N«tl  •  IPLUS 

_ GO.  JO  J.t _  ,  ,,  II. 

C  STORE  SIGN 

_ 74  XOUTPTIM+I)  •  RINPUTII)  _ _ 

l«l«l 

C TEST  FOR  END  OP  EXPONENT _ ^ _ 

77  IPIRINPUTIl^l)  -  lOLANK)  74*  7«*  74 

74  K0UTPT(N«2)  •  O 

_  KOUTPTtN^S)  .•  RINPUTII) _ _ 

GO  TO  792 

_ 7.91  ROUTPTIM+21  •  RINPUTII) _  . 

K0UTPTCM*3)  •  KINPUTIUl) 


M  •  Ml  ♦  12 

_  GO  TO  100  .  .  - 

C  COMPLETE  DECIMAL  FIELD 

_ 60  LOO  •  Nl  ♦  11  - - ; - 

00  61  Jl  •  N  •  LOO 

_■  RL  unuTPTiJli  •  Jk - 

M  •  Ml  ♦  12 

_  ._100  CONTINUE-  .  _ _ 

C  MRITE  ALPHANUMERIC  CHARACTERS 

.WRITE  OUTPUT  TAPE.  9»JOOO*  lADUlF.TIJI*  JaU  lEMDI - 

IPINUMOCP  -  lOEC  )  201#  201*  200 

aofl  M  «  i_  - 

Nl  •  Nl  «  I 

_  lENO  •  12*INUM0CP  -..A) - 

lOEC  •  NUMOCP 

GO  TO  904  .  _ 

201  REWIND  9 

REM  INPUT  TAPE  9*1D10*IPARRAYIJ)*J»1*  NUNOCPI 

REWIND. 9  - - 

300  L  •  M 

RETURN  .  ..  - - 

1000  F0RN4TIT2AI) 
lam  FfMIMATIAPll.RI 
END 


_DEC0Ct _ U. 

DECOCP  IT 

-DECOCP _ 16 

DECOCP  19 

-DECOCP _ 20 

DECOCP  21 


DECOCP 

39 

JIECDCP _ 

-.40- 

DECOCP 

41 

DECDCP  .. 

—  42 

DECDCP 

43 

DECOCP... 

44 

DECDCP 

49 

.DECDCP 

44- 

DECDCP 

47 

DECDCP  . 

.44 

DECDCP 

49 

DECDCP  .. 

90 

DECDCP  44 

-DECOCP _ 49 

DECDCP  70 

-DECOCP _ 71 

DECOCP  72 


DECDCP  74 

-DECDCP _ 79 

DECOCP  74 

-DECDCP _ 77 

DECDCP  76 

DECDCP  40 

-DECOCP _ 61 

DECDCP  62 

-DECOCP _ 63 

DECDCP  64 


DECOCP  64 

-DECOCP _ 67 

DECOCP  66 

-DECOCP _ 49 

DECDCP  90 

-DECOCP _ 91. 

OECOCP  92 

.  DECOCP  93 

DECOCP  94 

.  DECOCP  .  -  99 
DECOCP  94 

-DECOCP _ ^97- 

OECOCP  96 


EXPliRIHlNUL  CAU 


riXItRE  lE^'MRXTljilE  ^  ‘2S.0  OEG.  C.  RIxrURE  PRESSURE  7A0.0  PM.  PCRCURV 

ATCPS  CP^XVCEA  1C  CCPPLUEtH  CXIOUE  CPE  PCLECUIE  OF~FUrir  UCO 

. .  PCIE  FRXCTICK  Cx'vCEa' lA  (»rCP'Nt  ‘  O.ElCb 

vCLUPc  M~Hcii^CP  TleI  ^V.IaVae  04  cc.  per  crxp-moie*'  . . 

fcxpe  frcni  cipensions 

'•u'i"  ■  ^ 


fll  PEAK  HEICFIT  (PEPSUREC'lNITSI  Al.OO 

ilATlQiL  iPtXlUAFOl _ t  t _ i. 

CIAPETER  IPEASIREOI  4t.CCC  SC.CCC  2I.CCO 

_ mCPJ  ...  ICPil _ ft, _ CjufCjtflC_ft.4C«ft 

CIAPETER  tCP.I  O.AlSSt  O.ASCas  0.)1S20 

^  ACTt*LJ.EKC.TK  £E.  REFERENCE— ._jU0ftfi££_-£M« _ 


U.OOO  A. SCO 

C.60Q0C  fl.iCOCO  Q.44RSR 
C.EtOlJ  0.047SA  0. 

PEASURED  .UNGTH-CFJtEEEtEMCft. 


TIA  FFAK  HtIChT  li>fASLRgC  LMIIXI  A2.CQ _ 

STXTICA  (PEASCREOI  I  2  2  4.9 

CIAPETER  JPEPSUREftJ _ 4C,KCL..,_J0..Xftfi__2.l>.ftft0!  .il.50.0 _ A^lftfiL _ 

FEICFT  (CP. I  0.  0.2eCC0  0.4CCCO  C.AOOOO  O.ROOCO  0.99091 

CIAPETER  (CK.I _  -  O..EAT.EE.  X.4W2ft..  ft.9J920,  .C.20261_-0,09T9.6._0[,^ _ 

ACTUAL  LEN6TF  CF  REFERENCE  4.CCCCC  CP.  PCASLRED  lENSTh  CF  REPCKENCE  2AA.90C00  UNITS 


T9T  PEAK  NEtCFT  IPEASUREO  LNlTSt  49.00 

„  .  SIAllOA—UEASUREDJ _ 1 _ 2 _ Jl _ i _ S _ 

CIAPETER  IKEASUREOI  42.0CC  3t.C00  22.CC0  15.000  1.000 

_ PAIOFT  tCP.I _ _ C.2X.taC_C..4CCC0L  .0,40000-.  -D.80fift0_fl,3ft2M _ 

CIAPETER  ICP.l  0.43SU  0.44E01  0.932TO  C.224S4  0.1209R  0. 

AgfLAL  LtNCtP  OP  RFFFRgiigf  A.ttttt  CP. _ PEASLRfD  LENGTH  Cf  KEFERENCE _ 2A4.S0C0C  UIHU 

TIE  PEAK  HEIGHT  IPFASURfC  UHITSI  73. CO _ 

STATION  (PEAStREOI  1  2  9  4  9  A 

ClAPEIER-lPEASLBCftl _ _„42,CCfl _ 32.9ftC _ 24.000 _ 19.000  .,12.500 _ S.SOO _ 

FEICFT  (CP.I  0.  0.2CCCC  0.4CCCO  C. 60000  0.80000  l.OOOOC  1.09981 

eiAPgfXR  lew.  I _ fl.AHH  C.4R444  C^4i9A  C.2T1 14  0.iT324  0.OR2aA__O. _ 

ACTUAL  LENGTF  CF  REFERENCE  4.0000C  CP.  PEASUREC  t.tNCTh  OF  REFERENCE  265.90C00  UNITS 

n9~PEAK  HEIGHT  (PCASUREO  LNITS)  79.00 

_  STATICA  .  (PEASLKEOI _ _ 1_  2 _ JL _ 4. _ 5 _ 4 _ 

CIAPETER  IPSASLREO)  4t.5CC  32.500  29.CC0  16.900  11.900  9.90C 

hClCt.T  ICM.I _ a, _ Q.2QC0C  C^4CflC0_  0.40Q0C  __0.AQ0QQ  l.OOQOfl _ IJUUftl _ 

CIAPETER  (CP.I  0.42924  0.48944  C.3764S  C.2T872  0.1T326  0.08284  0. 

^  .  ACTUAL  LING  IF  CF  KEFERENCE _ 4.C0CCC _ CP. _ PEASURED  LENGTH  CF  KEFERENCE  ..  269.90C00  ..UNITS 

740  peak  FEICHT  (PEASUREC  LMTSl  79.C0 

STAtICA  (PEASLMEC.I.  _  1  2  3  4  5  4 

CIAPETER  (PEASLKECI  42.9CC  34. COO  26.CC0  19.00C  12.900  7.000 

FEICF.r  ..  JCP.l.  ..  0.  0.2CCCC  0.4CCC0  C.40CO0  O.BCOCO  l.OOCOO  1.11732 

CIAPETER  (CP.I  0.43319  C.9C492  0.3<T34  C. 28305  0. 18422  0.10428  0. 

ACfLALLENGtF  CF  RfFiRgNCP  A.CCCCC  Ci>, _ PtASURCD  LLN6T1j.CP_RER.ERENCE _ 2,68..,9.C.CQ.a _ UNUS. 

741  PEAK  height  (PtASUREO  LNITSI.  .  ..77,00. _ _  _  .  .  _ _ _ _ _ _ 

SfAflOA  (PEASUREOl  123454 

CIAPETER  IPEASUREOJ.,^.  .43.000  34.C00  2T.CC0  21.000  15.0CC  ,8.504 _ 

FCIGFT  (CN.I  0.  C.ECCOe  C.4CCC0  C.40C0C  0.80000  l.OOCCC  :.t9140 

eiAPITAR  ICP.l _ 0-AA349  0-94AX>  0-40374  C-Xi4QX_  Q.2243Q  0.I2T10  Q. _ 

actual  LENGTF  CP  REFERENCE  A.OQOOC  CP.  PEASURED  LENGTH  CF  REFERENCE  267.90000  UNIfS 


742  peak  height  (PEASUREC  UNITSI  E8.C0 

_ JTAT.10N_J.P.£A2UA£&i _ I _ 2 _ 9 - A - J - A - 2 _ 

CIAPETER  (PEASUREOl  44.CC0  37.CCC  30.900  24.000  17.SC0  13.000  7.000 

<CB,I _ (L _ C-7CCC0  O.ACflCQ  C.&OCQQ _ 0,80000 _ l.QOQDC  _1.2C0QQ  J.J1Q99 _ 

CIAPETER  (CP.I  0.49949  0.59121  0.45438  C. 35794  0.26071  0. 19367  C.1C428  0. 

_ ACIUAL_HNCIF-.CF_REF.tlUNC£ _ A.^OO.Oft___CM. _ P.EASUREC.LkNCT.H  .OF  .REFIRINCE _ 268,34400 _ UHUi. 


JkUK _ 

_  FUCL-FXCH.  . 

(CC./SECI 

OXlOANJ.FLOli.. 

(CC./SECI 

iTT.xceee 

..HOLE  FRACTXC1L„. 
1NM18IT0R 
a. 

¥CL0«..flCN _ 

(CC./SECI 

_ CCNE-ARE* _ P.LAU-SP.EED. _ 

ISO.  CP.I  (CP./SECI 

Q.JU£6ft _ P94.T3144 _ 

ECUIVRIENCE. 

RATIO 

1.10404 

736 

91.400CG 

137.3CCC0 

0. 

249.64099 

0.19899 

277.88917 

1.SI499 

M  IH  ■ 

i97.3ccee 

0. _ 

.X6Z,32U5. . 

_ ■0,.9.Z310..  ■■ 

9T4. T I4TI 

_ L..865U. 

798 

I24.0CCC0 

137.3CC00 

0« 

289.22602 

1.09941 

259.43590 

2.19091 

_ i24.ocege 

-  124.RCCeC  _ 

_ 0- 

273.T4458 

..i,10J£0 _ 

247.07991 

.  .  laJlUa. 

740 

124.00CCC 

116.9CCC0 

o« 

262.92146 

1.19208 

227.88642 

2.8I423 

fl. 

1.24648 

201.89848 

741 

124.0CCCC 

lA.oceoo 

0. 

229.22872 

1.43982 

199.2C607 

1.49300 

»Mic  cutifi  tii'u 

ccifFicifKift  rrSSciio  fcuivHfikCf  Mfio  nr*tmm~^ 

_ A1 _ U _ A1 _ _ CEVIAriCN  AJ .$fl«CL.._.FLAHC  SHifi _ 

•i•AfeTf  Cl  *  I.IMIE  Cl  01  4.Uf*f  01  l.iTIO  I.TIIOO  III.O^M 

miii  ^  looimiikCf  Riltticb  ^ocoicfcc  otvuT  ion  ’■“ 

Jtam  sMiC _ 0L4i>i  into _ _  _ : _ I 


04tu  Aca^imi  - 

■"  >uit  •iii'i  '"'  iNt-uiiCNiVusiir 

_ KVCHCCU  -  _ _ _  .  . 

tiuriill _ C _ Ei^l^CUL-hUKIU _ 0 _ lull 

_ CAi<LJOiRCt.j(Mitu _ I -  exmiiKeNUL  ccnciticns.  nurur _ i _ 

_.JlUKI..V(LCC.l.l.t-A.USTClC»<tOI’C1KIC  CCNCCNIAAf ICM  JI«TIG.  ..  lOJ.MfiT _ Cf!.«/Sfl 

rAiiriM  £um  matiiy - - - 2UjJuia _ «■/*« 

.MIfL.£ONCeiltM1JOPlJlT  SIOICNKmeiRIC  CCKfiltfCNS . . 0»72l|20£.lf  MOLECVUSA 

JUtI  Cf.llCtilTtUlOlLJJ.-CCIiOU10WS.  01  UAXIKUr.  UELCCIfJt _ OaOJMJE  20-  .MlfCtIUS/ 


APPENDIX  B 


Fortrar 


rogram  and  Sample  Printouts 
for  Routine  PSR 


ROUflNC  PSA 


MumiTiouisui  im 


III  NOVPLIi  I  •  1|  NOVNII 
ECrOA  (liOVMf  NOVPLI 


TAPE  )i  40i 

AAO  HAlTf  OUTPUT  TAPt  »>«  VECTOA  IWUVAW*  NOVAAI 
CO  TO  6S0 

C  CALCUIATIOM  Of  A1SID4IAL  SUWl  Of  SOUAWES  AWO  CAOil 
tSO  IPlJF^NSri  900t6»l|7AS 


»2  MKITE  OUTPU 
GO  10  410 


Ui 

esso4 

115 

ESS04 

116' 

E$$04 

11/ 

ESS'04 

118 

MO  WAITE  OUIPUT  lAPE  3.  SS 
ftO  UKfTE  OUTPUT  TAPE  At  JS*  (M 
760  WAITE  OUTPUT  TAPE  3»  40,  11  .VECTOAI I  «l*OVAAI . I* t .liQWWH 


__  782  06FK  «  VECIOAINOVPLiWOVPLI  -  l.O _ 

?«*1  IFIVECTOAlIilil  792ir<»4«01O 

y^*2“”wATft  OutPUT  TAPE  J$  - - - 

0  _ I _ 


78J  FOAHAT  lAlM  EAAOA  KESIDUAL  SQUAAE  VAAIAALE  lAillN  IS  NEOATlVEiPAOi 


ESSOA _ iJ« 

"¥5S04  140" 

eSS04  141 
feTSS4  142" 
ESS04  141 


JIT:T:fH8TTTi|  irXWf  !:1 7  A'n  HJl  A« 


810_SICMAUI  •  SUAIF  IVECTOA  IlilU 
tiOO  V^TOKI 


•40  IF  tlFCOCN)  ^Ot  WTOi  1000 
TAPE 


474  N0VM2  • 
•79  OU  ••» 


6V0  WAltC  OUIPUT  TAPE  3i  40i  I 1 • VECTOAl I iHOVAA) » l-l 
«  I 


iTiI<lWTS^4riT.TI - 


1002  NSTPMi  •NOS  TCP  - 


CO  TO  1781 

1020  VNIN  •  0.0 _ 

io7o"vnaa  ■  0.0 


1090  1  •  liNOVMI 

1041  IF  IVECfOA  lltlll  1042U090 


1042  WHITE  OUIPUT  TAPE  !•  1044i  I*  NOSTEP 
1045  CO  TO  010  _  _ 


1040  IFlVECfOAlIill  -  TOU  1050fl0«0il0»0 
1080  VAR  •  VECTQRIltNOVAH)  •  VECTUAINOVAA 


1120  INDEXINOINI  •  I 

1130  COCNINOINI  ■  VECTOAITiNOyAA)  •  SICNAINOVAAl 


180  NOMIN 
00  GO  TO 


1210  VMAX  •  VAN 

1220  NOPAA  •  I _ 

1090  CONTINUE 


_ on  TO  810 

U40  white' UUIPUT  TAPE  7t49i  SK 

1240_W_TO  1390 _ 

1245  IF  IlFCNSri  OUOti290il244 


ESS04  174 
ESS04  175 
ESS04  1/4 


ESS04  100 
ESS04 _ ^101 


ESS04  102 


rno 

404 

IF  (VNINI  1140*11701004 

WHITE  OUIPUT  TAPE  3.  404 

ESS04 

ESS04 

190 

ivi 

00  10  VIM 

ESS04 

192 

1170 

VNIN  •  VAR 

ESS04 

193 

8] 


NUHSA'ifO  kCSIcAllCH  F(.AME  SKfcO  KECHCSSlOM  •  XOUilNE  1922  '  *  M301PlCiri0««  i  •  iItJN 

Tii6iiftViliio1»esufc¥'?6Ri»l»Airuir“TE$rM  = 


NUHrtCR  UF  FMOMIENS  IN  THIS  ftUN 
ACciimitE  uMA^soukcriMMaEii 


HlbHEST  CONTAIIMIOA  CODE  NUNdEA  IN  USE 
'ACCEPrAHE  fXFEAlNENUi  CONDMlbNI  NUMIU* 


ENTAUS  TO  6E  ACAD  FAOH  IH£  IN'hJT'bATA  CUDS 


Lisr  SlbNlFtCANCE 

-  UNACCEPTAALE  data  group  SEAJA^NUNAEAS 
PTAALE  FUEL  CLASS  NUNAEAS 


-  ACCEPTAALE  fuel  CLASS'GAOUP  NUNAEAS 


IMililiilf 


ION  CONIAOC  DATA 
HON  CONTAOL  DATA  


37  4  40  O 


LIST  OF  COUNT  TEST  CAITEAIA 


LIST  UF  ACCEPTCa  DATA  SERIAL  NUNAEAS  ANO  FUEL  NAMES 

_ 2_t.lMAJlt _ 1 _ i  .PROPANE _ juMOUlANfi— _ 

tt  NOCTANE  9  N-OECANE  10  N-UNDECANE 

_ 14^^2HHPRQPAHE  _  ia_2rME-P_fcNTANE _ 

20  234HbPCNrANb  21  224N£PLNrAHE  22  22SNCHeXANE 

_  26  NECtCPENrANt  27  NECVCHbXANE  20  12HEC,VC.HCX?C.. 

3)  ClS*eurE9E-2  J4  rx^QUlENE  2  35  N-PUNIENE*! 

_ Ji  2rMi;dUr^\l£^  _  AQ  2«NEeUlENEt2  41  JkMEP£NXEN-aC  , 

4b  CVCLOHEXCNE  47  4VI^CVCHEX->i  40  CVC0CI0ENE1:» 


_S  NPENTANE  ,  _ 
11  M-IRIUEC4NE 

»ll_UNE.SUJA'i£— 
23  CYCLOPAOPASE 
29  12NECVCHEX-T 
36  NPENTENE>2 
42  4NEPE9rE9*2T 
49  PftbPV9E 


6  N-HEXANE 
12  ISOOUTUE 

..llL.22.NfAM.U9| _ 

24  CVCIOPENTNE 
31  PROPYLENE 
37  ISO’OUTYLENE 
43  244MCPENrNEl 
50  1-BurVNE 


7  NHEPTANE 
13  ISO-PENTANE 
_!!>„  24H{P.tNrANf._.. 

29  CVCLONEXANE 
32  AUTENE-I 

30  IME-AUTENE-l 
4S  2HESUIDENE13 
52  TOLUENE 


•  59  I-LUIFLEFHAN 

_  _  66_NaUCNLUAlQE _ 

72  l-2CLA0rANE 


60  L-4F0UrANE  61  UlFPENIANE 

.6XJ2&y;ilkUS10£ _ .6t-L-j;L«fcN.UL«_ 

r4  2-CL-PAOFENE  75  ALLVLCHLAIOE 


OS  MC-PH-REIONE  66  NE4FPHKET0NE  07  NEOEN/OATE 


97  1-2CL2FPAPAN  99  M-F-TOLUENE  111  N-AU-ANINE 

_116_  HHEPAOeNlNE  117  ISOPACHLiL10E_119^  I SPAPtLANlNE, 

124  OATHOXVLENE  125  NETAXYLENE 

130  OtISOFASLFOX  139  N-OU-SULFIOE  140  ISOO^ACAFTN 


63  NE-CHLOAlOe  64  ET-CHLORlOE  65  NPKCNLOAIOE 

__frl-m2.N£.e!ipAi»B _ 10.XCk2Ni£i£A.H£-  .  71_lr4CLSWIA!lt_  .... 

76  2CL0UrENE-2  77  ETHANOL  70  N-PAOFANOL 

_  orOJ-HE-EmA  ■  03  VIrHE-ElitEA _ 04,  P-F-ANISPLE _ 

00  M-CAESDL  89  l-OCLPAOPANE  90  1-2CLPA0FANC 


112  4NE-PYAlDtNE  113  ETHLNDlANINE  114  2NE-FVAID1NE 
_UP-JiaiaUil£iii£..a2.lJXAJ!E.UNll^lE-.i.T^  SECAUM.l|U!IL_._ 
127  l-F-PENTANE  120  26-LUriOlNE  129  NPACNLOAIOI 

135  1-F-HEKANE  1 3.#LAimJI.ARl.N.E_  WPASUS 

142  FAOPANE 


PONSANtU  iU;5CARCH  CURPIIHA'U  O'!  FLAKE  SPeEb''KECK£SilON'  -  " 'A0UnNfT922  '  “'^MDOIF  iVaTIVn  T  -  *’aUN”  3 


'OEPENUENl  VAHlAdLE  -'i  uVax/s'CHAX'xIO''’*^'' 


L~r.T  UF  INUEPCnOENI  CUNtRtOUrORS  WITH  ThCIR  COAKESPONDING  HEUAESSIUN  INDICIES 

_ 0  X  -2 . > _ _ ^ _ * _ ^ _ •• . 

*37  PRIM-H  30  SeC-M  39  TERI-H 


84 


I  fiiiii  CAiu  run  titso  Ri.uAissiM  mojaam 

Vii.l»WKE  KCiiUSSirM 

7niibi.6M  w  I  “  .  ^ 

"^'o  or  OAU  •  uti  .  . . . 

'no  Of  VA:<1ABIES  • . 4 .  . .  . .  . . 

~n.CHIEO  OC6.'ETs“'uirFkiE6iiN^  '“UArOU  ” 

VVtVEt  lb  iNTCA'VAAIABLC"*”  0.6UOU2  . 

TleVwL*  ib’iufNbW'vA'ituBle ‘'i'*"  biobbol  ' . . 


SfANOAKO  EHi(OA'of~V~«  YfoTiiTWS 

“TfeV  NO,  ‘  i"  . .  ■ 

_ V.AKI*ALE,.|N.rCAlNft _ 1 _ _ _  ... _ _ 

F  LEVEL  )tS.A6A> 

_ SUNOARO  EM^OR  OF  V  •  llA.i.lA9  . 

CONSIA'^r  0. 

_ _ VARIABLE  ..  COEfflCIENf  SID  ERROR  OF  COEf  ...  ..COkF/STA  EHUOR 


JHz.  L- 

^•2AA60E..02 

_ ll.iS2Q0JL.Ql_, 

.Q.1.7.Tjr.Sl;.02..  _ 

jltEF  NO.  .. 

.2. . . 

VARIABLE 

ENfERING 

2 

SfANDAtO 

EMROK  OF  V  • 

6).»0B) 

variable 

CUEFFlClENf 

SVD 

ERROR  OP  COEF 

COEF/SfO  ERROR 

X-  t 

O.ia/SBc  02 
._D.a23LBE.  02 

0.L04B9E  01 
.V.1UKE..UI. 

0.17BA4C  02 

. U...itt4.V*C.  02 _ 

Mn. 

_.J  . .  . . 

VARIABLE 

_ A. LEVEL  . 

SfANOAKO 
.j _ coNSfAvr 

ENrEKlNU 

.  ..  o.soo 

ERROR  OF  V  - 

,0. . 

VARIABLE 

) 

A).72AT 

COEFFlCIENr* 

S  t  o' 

’lRROr  'oF  COEF 

COCP/SfO  ERROR 

.JIOMFLCICO  . 

1 

.  »-  2 

X-  ) 

...  )  SfERS  OF 

O.loBSAE  02 
0.I2241E  02 
-O.BAAOAL  01 
lkCURfcSSlil.1 

O.DSliE  Ul 
0.tlBl)2C  01 
O.llUtk  02 

U.14S40C  02 
O.IOITEr  02 
•O.TOAOBC  00 
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APPENDIX  C 

Fortran  Programs  for  Routines 
FSRTL  and  FSRDM 


87 


RCHINO  4 

.  TIKE  CAM . . 

READ  INAUT  TAAf  2«  1000 
.  WRITE  OUTAWT  TAAE  4«  1000 
NUMOCA  OP  6A0UAS  OW  TARE 
-R£AO  lhA.UT^itft.2>.ia0t»  Ifl 


WRITE  OUTAUT  lAAE  4*100li  lAl  ASRTl. 

MAIN  DATA  fiAOUAS _ ASRK- 

00  20  lAl  ASRTL 

READ  INAUT  TAAE  2*I002*  ISAII».  tm»  AMS*  lAMCi  lAllO*  lAMRt - ASRTL.. 

1  tNOS«  INEC.  I  0IX}«  R«2«4I*  lAAOC* ILAOCCNlKl  t '  AOCIUI  tMlt  ‘  ASRTL 

_2 _ lAAGC-l...  _  A.S<tK. 


ASRTL  9 

-ASRTl _ 4 

ASRTL  T 

-ASRTL _ • 

ASRTL  9 

J>S1UL. 


20  WRITE  OUTAUT  TAAE  4«2002*ISANtANl»  ANSt  lAWC*  1AN0»  IAIM« 
.  1  INOSt  INCC*  1  0IR|«  ll*2«4l.» 

2  IAAGC  I 

_  CONTRIBUTOR  NAMES _ _ 

READ  INAUT  TAAE  2»  100S»ICNLtn»l»l*  2001 
_ WRIJE-.QUTAUl..lAAf  4i 200%, IfeMI  1 1  Will »  2001 


OEAENOCNT  VARlABtES 

READ  INAUT  TAAE  2.  1001*  (ONNI  |  |.•.4•1•4I _ 

WRITE  OUTAUT  TAAE  4,  200S*  tONNIt)tl«lf  41 

ENO  FILE  4  _ _ _ _ 

REWIND  4 

.-WRUi.4WlAUT  TAAE  >>  TOPS _ 


ASRTL 

_ ASRTL- 

ASRTL 
.  ASRTL  . 

ASRTL 

-ESRTt- 


CALL  EXIT 

1000  AORNATIlHiaiX*  AOH. _ 

1  I 

loot  AORMATIIHO*  141. _ 

1002  AORNATIU*  2A4*  014.  2A12.9  /  SAI2*Si  I4»  I4t  A12.S  /Il4i  AI2*4t 


•  ASRTL 
-ASRTL- 
ASRTL 
.ASRTL  .. 
ASRTL 
ASRTl 


_ t4<  fl2>4<  14,  AI2.4»  Ut  AI2.4  1  I 


ASRTL 

-ASRTL.- 

ASRTL 

..ASRTL- 

ASRTL 

-A.SJlJ.l_ 


1003  A0RMAI(12A41 
..2002  AORNATUHOr  t4»  .2A4*-S14t- 
■  1  •  14.  E12.4,  14*  C12«4« 

_2003  AORNATIIHO*.  914XtiUI/JlM.  «.  J 
1  4XA4«  4RA4*  4XA4  I  I 


ASRTl, 


tOUTINI  PSIlON 


C  NOGCOr  -  NIMtU  OP  OAT*  6K0UPS  CUAMNTLV  OH  TAPI 

».C_- _ NOCCHG...  -  —  !-  J«UNIEA.0P.  DATA.6AOUPl.*LTfll£0 _ 

e  I0CCH6UI  -  CODE  NUMOER  OP  ALTERED  OAT*  CROUPS 

.€ _ NAOO  _ >;.JHMilR  OP  DATA  GROUPS  ADDED  .TO -NASTCR..OCCM _ ^ 

e  ICONAM  -  CONTRIRUTOR  NANS  PUNCH  SWITCH  l-PUNCH  EHINIT 

.JC _ lOPNAN _ -  WABlAAtg  »UilfW  AMlTrW  ' 

OINENStON  tOCCHGUODOO)  •LF6CCNI90)*  P6CI.I90>|  GNLIEOOtf  DNNUI 

_ l.t  0(*l  _ _ _ _ _  -  .  - 

REWIND  A 

_C _ READ  INPUT  .CONTROL-DATA  ..  - 

READ  INPUT  TAPE  2.  2000*  NOGPOT*  NDGCHGt  ICONANt  IDPNAM 

_ iPt,NOCCHGt  I0»  2fl>-iD _ _ _ _ 

ID  READ  INPUT  TAPS  2«  2001*  (IDGCHGUli  R«l*  NODCM  I 

„C _ READ  TITLE  AND  NUH2ER  OP.  GROUPS.PROiL.msm..TAPi  AMD  PIllCM 

20  READ  INPUT  TAPE  «•  1000 

_ _  READ  INPUT  TAPE  A* '  lOOiV  NOGCOT  „ .  _ _ _ 

WRITE  OUTPUT  TAPE  IA»  1000 

_ HRtJ1...0UTPiU^lARE _ 2«-LIMIfi _ 

L  •  1 

_ .WRITS  OUTPUT  TAPE  lA.  lOOl,  NDGCOl. _ 

WRITE  OUTPUT  TAPE  2«  lOOl.  NOGCOT 

_ IPl  NOGCOT  NOGPOTI  .100«..l0«.  .20 _ 1 - 

20  00  AS  1«U  NOGCOT  i 

1  INOS*  TNEC»  I  DU)*  K>2*AI*  1PPGC»ILPGCCNU)*  PGCLU)iK«l* 

_ 2  IPPCC  )  .  _ _ 

iPt  I  *  NOGPOT)  22*  22*  2A 

22  IPt  1  -  lOGCHGtL).)  .40*..22»-40 _ L: - 

22  L  •  L  ♦  1 


t  INOS*  INEC*  I  DU)*  K«2f*)*  IPPGC*  ILFGCCNIK)*  PGCLUItR*!* 
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